iMechanica - Comments for "Derivatives of the invariants of a tensor"
https://imechanica.org/node/1403
Comments for "Derivatives of the invariants of a tensor"enAbout the derivatives done by biswajit
https://imechanica.org/comment/27729#comment-27729
<a id="comment-27729"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18557#comment-18557">Dear Biswajit,
when</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>I am agreeing with Attila. 'p' it self is a fucntion of sigma11 and biswajit has not taken it to account. I have deirved the derivatives of third invariant of deviatoric stress tensor with respect to direct terms and also by indirect terms of deviatoric stress tensor. but those expressions are not similar to the above expressions..since p is also a fucntion of direct stress terms; derivatives of the above expression should be considered as a derivatives of a mutiplication...</p>
<p>Gayan </p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 05 Oct 2015 12:27:09 +0000Gayan Aravindacomment 27729 at https://imechanica.orgDerivative, Gradient, etc
https://imechanica.org/comment/18625#comment-18625
<a id="comment-18625"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18583#comment-18583">derivative, gradient etc</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Dear Prof. Acharya,
</p>
<p>
In my previous post I actually meant derivatives of real valued functions of (deviatoric) symmetric tensors (such as invariants of that tensor as in our discussion) and not of course tensor valued functions. By the way, your comments concerning the gradient and derivatives of real valued functions of (deviatoric) symmetric tensors are very helpfull.
</p>
<p>
Thanks
</p>
<p>
Mohsen
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Sat, 03 Mar 2012 18:45:34 +0000M. Jahanshahicomment 18625 at https://imechanica.orgderivative, gradient etc
https://imechanica.org/comment/18583#comment-18583
<a id="comment-18583"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18576#comment-18576">Derivatives on Constraint Manifolds</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
In everything you say above, I suspect you mean derivatives of real valued functions of (deviatoric) symmetric tensors and not derivatives of the tensors themselves. For the answer to the question you ask, see my comment on Feb. 26 in this post entitled"Re:J3 derivative", item 2).
</p>
<p>
The reason that the gradient of a differentiable real valued function on a finite dimensional vector space with an inner product has to be an element of the space is not because one 'wants' to have it that way - it is a theorem. At this point, it may also be a good idea to realize (and if you are not familiar with these notions, to look up a good book on advanced calculus, e.g. Fleming - Functions of Several Variables) that the derivative of a real valued. fn. on a vector space is a linear transformation on the space to the reals (often called a linear functional). Thus, it is not a member of the space itself. However, there is a very important representation theorem for linear functionals (the proof of which is easy in the finite dimensional case) that shows that every linear functional L on the vector space can be represented by a *unique* element l *in the space itself* such that
</p>
<p>
L(h) = l.h for all h in the vector space. If L was the derivative of the real valued fn. we are thinking of, l is the gradient of the function.
</p>
<p>
In our discussion, det(S) and (1/3)tr(s^3) are equal real-valued functions on the set of deviatoric symmetric tensors. Thus it only makes sense to ask that their derivatives, or their gradeints in the above sense, be equal in this space. But the gradients are then members of the space of deviatoric symmetric tensors.
</p>
<p>
From the above definition you can now see why, if one has the gradient of the 'same' function on the whole space, subtracting off the projection of the gradient on the normal to the subspace will give the right answer for the unique gradient of the restriction of the function on the subspace.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Wed, 29 Feb 2012 15:00:04 +0000Amit Acharyacomment 18583 at https://imechanica.orgDerivatives on Constraint Manifolds
https://imechanica.org/comment/18576#comment-18576
<a id="comment-18576"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18562#comment-18562">derivatives on constraint manifolds</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Dear Prof. Acharya,
</p>
<p>
If I'm correct, this means that we want to keep the derivative onto the constraint manifold and if not so, the constraint would be meaningless (there would be no difference between derivatives of deviatoric symmetric invertible tensors and derivatives of symmetric invertible tensors). But as a consequence of this constraint all the derivatives of deviatoric symmetric invertible tensors should be traceless (such as dJ_3/ds in our discussions). Should it be the case?
</p>
<p>
Mohsen
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Wed, 29 Feb 2012 06:52:23 +0000M. Jahanshahicomment 18576 at https://imechanica.orgre: J3 derivative and dp/dsigma_11
https://imechanica.org/comment/18568#comment-18568
<a id="comment-18568"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18557#comment-18557">Dear Biswajit,
when</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
The point of the above was to show one way in which brute force can be used to get solution 1.
</p>
<p>
We have assumed dJ3/dp = 0. But, as Amit had pointed out somewhere else, dJ3/dp derivative is ill defined because J3 can take an infinite number of values for a given p.
</p>
<p>
Option 4, dJ3/ds = Dev[s^2], gives you the tracefree solution that Amit points to, i.e, the solution that is tangent to the DevSymm manifold.
</p>
<p>
-- Biswajit
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 20:14:00 +0000Biswajit Banerjeecomment 18568 at https://imechanica.orgderivatives on constraint manifolds
https://imechanica.org/comment/18562#comment-18562
<a id="comment-18562"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18555#comment-18555">Derivative of J_3 of a Deviatoric Tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Your remark about the derivative being tangent to the constraint manifold is correct - see my comment entitled "Differentiation" earlier on in this post (from a long time ago).
</p>
<p>
The equation of the manifold is tr(s) = 0 (it is actually a linear subspace of the space of symmetric second order tensors) - for the invertible part, one has to do a little song and dance but things work out as if you did not worry about the invertible part.
</p>
<p>
The way you calculate the normal is as one learns in advanced calculus - you calculate the gradient fo the function that defines the level set. Assuming the Froebenius norm on the space of second order tensors (i.e the trace inner product metric) is involved in this calculation as in forming the unit normal.
</p>
<p>
Look at my comment on "differentiation" above and let me know if that clarifies your question about why the 'derivative' of J_3 on a manifold should be tangent to the constraint manifold.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 15:15:22 +0000Amit Acharyacomment 18562 at https://imechanica.org Re: J3 etc.
https://imechanica.org/comment/18561#comment-18561
<a id="comment-18561"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18558#comment-18558">My question is related to</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
I gave you the unique solution for the gradient of J_3 as a function of symmetric deviatoric invertible tensors, as well as symmetric deviatoric tensors. J2 does not have to be zero for that answer to be valid.
</p>
<p>
I haven't checked what the text books exactly state for the gradient of J2 on the space of symmetric deviatoric invertible tensors. However, where your derivation goes wrong is in interpreting your first equality in (11), <img src="http://imechanica.org/files/images/J311First.png" alt="First equality" width="124" height="48" />, as valid on the whole space of symmetric tensors. This is not true. The second equality in (11), <img src="http://imechanica.org/files/images/J311Second.png" alt="Second equality" width="107" height="51" />, is true in the whole space. if you restrict attention to the subset where J_3 is indeed equal to 1/3tr(s^3), i.e deviatoric tensors, on that set the two 'derivatives' match indeed, as I showed.
</p>
<p>
May be it is easier to see what is going on in a simple setting - take two real valued functions defined on the x-y plane. Let these two functions agree on a straight line in the x-y-plane - so think of two arbitrary surfaces that intersect on a straight line. Would it be reasonable to expect that the gradients of the these two functions should be equal everywhere just because their values match on the line? Even on the line, would it be reasonable to expect the 'derivative' of the functions to be equal in directions not along the line on which they have common values? Of course, if you restricted attention to the functions only along the line then it is eminently reasonable to expect that the 'derivatives' of the restrictions of the two functions on the line should be exactly equal, and indeed this is what happens.
</p>
<p>
Makes sense?
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 15:05:00 +0000Amit Acharyacomment 18561 at https://imechanica.orgMy question is related to
https://imechanica.org/comment/18558#comment-18558
<a id="comment-18558"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18553#comment-18553">derivative of J_3 of a deviatoric tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
My question is related to real-valued symmetric second order invertible tensor <strong>s</strong>.
</p>
<p>
Only one solution should exist for the derivative of<em> J3 </em>with respect to <strong>s</strong>.
</p>
<p>
</p>
<p>
In this general case, <em>J2</em> is not zero. However, Method1 and Method2 are the same only if <em>J2</em>=0.
</p>
<p>
But, I am interested in the case when<em> J2</em> not =0, <em>J3</em> not =0.
</p>
<p>
</p>
<p>
Consequently, I am still surprised why exist two alternatives for this derivative in textbooks.
</p>
<p>
Where is the "tricky" step in the derivations causing two solutions?
</p>
<p>
</p>
<p>Attila</p>
<p>
</p>
<p> </p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 09:28:00 +0000Attila Kossacomment 18558 at https://imechanica.orgDear Biswajit,
when
https://imechanica.org/comment/18557#comment-18557
<a id="comment-18557"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18550#comment-18550">Option 3: J3 derivative</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Dear <a href="http://imechanica.org/user/1095" title="View user profile.">Biswajit</a>,
</p>
<p>
</p>
<p>
when you compute d(J3)/d(sigma11), I think you forgot to compute d(p)/d(sigma11). p itself is obviously depends on sigma11. Therefore a term 1/3 is missing in that line, I think.
</p>
<p>
Attila
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 09:18:22 +0000Attila Kossacomment 18557 at https://imechanica.orgDerivative of J_3 of a Deviatoric Tensor
https://imechanica.org/comment/18555#comment-18555
<a id="comment-18555"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18553#comment-18553">derivative of J_3 of a deviatoric tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Dear Prof. Acharya,
</p>
<p>
The expression
</p>
<p>
dJ_3/ds = s^2 - (1/3)(s^2:I)I
</p>
<p>
actually means that we are computing the tensor tangent to the constraint manifold. This is the key difference between the derivations in previous posts. But the questions is that how we can compute the normal to this manifold (or what is the equation for this manifold) and why the derivative of J_3 should be tangent to the constraint manifold (in some sense a covariant derivative on constraint manifold)?
</p>
<p>
Regards
</p>
<p>
Mohsen
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 08:45:07 +0000M. Jahanshahicomment 18555 at https://imechanica.orgderivative of J_3 of a deviatoric tensor
https://imechanica.org/comment/18553#comment-18553
<a id="comment-18553"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18539#comment-18539">details</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
What's going on here is interesting indeed.
</p>
<p>
Without loss of generality, let's consider your (1) to be valid for symmetric second order tensors and (2) to be valid for symmetric, invertible second order tensors.
</p>
<p>
Consider your eqn. 11.
</p>
<p>
First equality: <img src="http://imechanica.org/files/images/J311First.png" alt="First equality" width="124" height="48" />
</p>
<p>
The first equality there holds if you view J_3 and tr(s^3) to be defined as functions *ONLY* on the space of symmetric deviatoric tensors. If you go off the manifold defiend by tr(s) = 0, then that relation does not hold.
</p>
<p>
Second equality: <img src="http://imechanica.org/files/images/J311Second.png" alt="Second equality" width="107" height="51" />
</p>
<p>
However, the second equality is written considering tr(.) as a function on the space of symmetric tensors, and in that space that is the correct definition of the derivative of tr(s^3) following (1). It can be shown by being careful about the meanings of the derivative and the gradient of a real valued function on a vector space that the gradient of tr(s^3) on the manifold defined by tr(s) = 0 can be written down by subtracting from the gradient of the same function on the whole space of symmetric second order tensors its projection on the normal to the constraint manifold. All this amounts to simply subtracting off the projection fo s^2 on I/(sqrt 3), so that the correct result 12 should be
</p>
<p>
<br />
<img src="http://imechanica.org/files/images/J3DerivGen.png" alt="J3 derivative in symm dev space" width="294" height="38" />---- (**)
</p>
<p>
This result holds for all symmetric deviatoric tensors, whether invertible or not.
</p>
<p>
Now let's go to your Method 1: The definition 5 and formulae 2 and 7 require J_3 to be defined only on the space of invertible symmetric tensors (no deviatoric requirement).The gradient of J_3 on the space of invertible symmetric tensors is indeed given by (7) - the domain is not a vector space but an open set of the space of symmetric tensors and this suffices. If you now ask for what the gradient of the J_3 should be on the space of invertible, deviatoric symmetric tensors, then again you take the formula in 7 and subtract off its component along I/(sqrt 3) and you get the eqn (**) above, by using your (8) and (9) which is valid for symmetric deviatoric tensors.
</p>
<p>
So d J_3/ds on the space of symmetric, invertible deviatoric tensors is the same no matter what formula you use. Of course, your corrected method 2 is more general because it gives the formula on the whole space of symmetric deviatoric tensors.
</p>
<p>
So I would summarize by saying that what your are looking for when you write dJ_3/ds is the gradient of the function det(s) on the (linear) manifold/subspace of the vector space of second order tensors defined by tr(s) = 0. This requires a little care, but it all works out consistently
</p>
<p>
I hope this helps.
</p>
<p>
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 04:13:00 +0000Amit Acharyacomment 18553 at https://imechanica.orgOption 4: J3 derivative
https://imechanica.org/comment/18551#comment-18551
<a id="comment-18551"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18539#comment-18539">details</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Another possibility :)
</p>
<p>
<img src="//imechanica.org/files/images/J3DerivOpt4.png" alt="J3 deriavtive option 4" width="663" height="470" />
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 03:48:04 +0000Biswajit Banerjeecomment 18551 at https://imechanica.orgOption 3: J3 derivative
https://imechanica.org/comment/18550#comment-18550
<a id="comment-18550"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18539#comment-18539">details</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p><img src="http://imechanica.org/files/images/J3Deriv.png" alt="j3" width="650" height="773" /></p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 28 Feb 2012 02:48:00 +0000Biswajit Banerjeecomment 18550 at https://imechanica.orgdetails
https://imechanica.org/comment/18539#comment-18539
<a id="comment-18539"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18532#comment-18532">Re: J3 derivative</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
I made a detailed derivation steps regarding to my question.
</p>
<p>
I think method 2 is the correct solution for the derivative of J3.
</p>
<p>
But what about Method 1?
</p>
<p>
</p>
<p>
Attila
</p>
<p>
<img src="http://www.mm.bme.hu/~kossa/j3der.png" alt=" " width="735" height="1310" />
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 27 Feb 2012 07:29:10 +0000Attila Kossacomment 18539 at https://imechanica.orgRe: J3 derivative
https://imechanica.org/comment/18532#comment-18532
<a id="comment-18532"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18531#comment-18531">Re: J3 derivative</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Amit
</p>
<p>
Agree with your caveats (1). It's been a long time since I thought about these things and I can't recall where I first saw that result. But it did look strange and hence the contortions to get to it.
</p>
<p>
(2) is quite a bit deeper and I hadn't thought in those terms. That's worth more examination from others in this forum.
</p>
<p>
</p>
<p>
-- Biswajit
</p>
<p>
</p>
<p>
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Sun, 26 Feb 2012 22:04:35 +0000Biswajit Banerjeecomment 18532 at https://imechanica.orgRe: J3 derivative
https://imechanica.org/comment/18531#comment-18531
<a id="comment-18531"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18520#comment-18520">J3 derivative</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
I admit that I have not checked the solution you provide and neither the references you give. But here are a couple of thoughts you may want to comment on:
</p>
<p>
1) Biswajit's formula above is valid only when the symmetric deviatoric tensor where the derivative is being done is invertible. A general traceless tensor need not satisfy this condition.
</p>
<p>
2) The other question that just occurred to me (and I have not thought about it much at all) is, if we think of J3 in your case as defined on the space of symmetric deviatoric tensors, which is a vector space (as far as I can think of), then calculating the gradient of a differentiable real valued function on this space should yield an element of the space itself - i.e. a symmetric deviatoric tensor. So as you say, s^2 need not be deviatoric, so what are you thinking of as the domain of the function J3 you consider - symmetric-deviatoric tensors, or only symmetric tensors, or invertible-symmetric-deviatoric tensors?
</p>
<p>
I note here that if we were to consider the fn J3 to be defined on the domain of invertible, symmetric deviatoric tensors the formulae that both you and Biswajit write will have to be corrected to remove their hydrostatic parts to obtain the gradient and then they would yield identical formulae.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Sun, 26 Feb 2012 18:41:40 +0000Amit Acharyacomment 18531 at https://imechanica.orgIt is obvious.
I didn't
https://imechanica.org/comment/18524#comment-18524
<a id="comment-18524"></a>
<p><em>In reply to <a href="https://imechanica.org/node/1403">Derivatives of the invariants of a tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
It is obvious.
</p>
<p>
I didn't write Tr[s].
</p>
<p>
Tr[s^2] and Tr[s^3] are not 0.
</p>
<p>
Attila
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Sun, 26 Feb 2012 07:48:33 +0000Attila Kossacomment 18524 at https://imechanica.orgRe: J3 derivative
https://imechanica.org/comment/18521#comment-18521
<a id="comment-18521"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/18520#comment-18520">J3 derivative</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
s is the deviatoric part of sigma and Tr(s) = 0.</p>
<p>
P.S. This comment is irrelevant. Please disregard.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Sun, 26 Feb 2012 04:10:00 +0000Biswajit Banerjeecomment 18521 at https://imechanica.orgJ3 derivative
https://imechanica.org/comment/18520#comment-18520
<a id="comment-18520"></a>
<p><em>In reply to <a href="https://imechanica.org/node/1403">Derivatives of the invariants of a tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Dear All!
</p>
<p>
</p>
<p>
There is a tricky question: why do we get different result for the derivative of J3 with respect to s if we start from the definition J3=1/3*Tr[s^3] ?
</p>
<p>
<strong>Solution 1:</strong>
</p>
<p>
given above in this blog entry in eq (32) as: s^2 - 1/2 Tr[s^2] * 1
</p>
<p>
<strong>Solution 2:</strong>
</p>
<p>
derivative of 1/3*Tr[s^3] with respect to s is simply s^2
</p>
<p>
[See: Jirásek, M. and Bazant (2000): Inelastic analysis of structures, p. 653 eq. (D.55); or see<br />
Zienkiewicz & Taylor (2000): The nite element method Volume 2, p. 435 eq. (A.25)]
</p>
<p>
[See: Itskov (2009): Tensor algebra and tensor analysis for engineers 2nd ed., p. 124 eq. (6.52)]
</p>
<p>
</p>
<p>
Best Regards,
</p>
<p>
Attila
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Sat, 25 Feb 2012 23:52:08 +0000Attila Kossacomment 18520 at https://imechanica.orgx + h off the shell surface
https://imechanica.org/comment/2854#comment-2854
<a id="comment-2854"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/2853#comment-2853">does the point leave the surface?</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Nachiket,
</p>
<p>
Yes, you are correct.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Wed, 16 May 2007 01:04:23 +0000Amit Acharyacomment 2854 at https://imechanica.orgdoes the point leave the surface?
https://imechanica.org/comment/2853#comment-2853
<a id="comment-2853"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/2848#comment-2848">Differentiation</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
You wrote:
</p>
<p>
For geometric understanding, consider a scalar function φ defined on a<br />
two dimensional shell surface in 3-d space. Let h be a tangent vector<br />
to the shell at x and let the shell not be flat at x. Then φ(x + h) is<br />
not strictly defined for every non-zero tangent vector h, however small.
</p>
<p>
I think this is because the point (x+h) leaves the surface of the shell. Am I correct?
</p>
<p>
-Nachiket
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Wed, 16 May 2007 00:51:55 +0000Nachiket Gokhalecomment 2853 at https://imechanica.orgDifferentiation
https://imechanica.org/comment/2848#comment-2848
<a id="comment-2848"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/2841#comment-2841">Re: Tensor derivatives</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Biswajit,
</p>
<p>
Consider the invertible tensors I and -I. Their sum is not invertible. Hence. not a vector space.
</p>
<p>
Differentiation of a function at x requires the value of the function to be defined at x + h for h sufficiently small. If the domain of the function involved is a vector space this is not a problem, and it is not a problem even if it is not the entire space but only an open set of a vector space.
</p>
<p>
For geometric understanding, consider a scalar function φ defined on a two dimensional shell surface in 3-d space. Let h be a tangent vector to the shell at x and let the shell not be flat at x. Then φ(x + h) is not strictly defined for every non-zero tangent vector h, however small. For the same reason the definition of the directional derivative also changes to
</p>
<p>
d/ds φ(f(s))|s=0, where f is a curve on the manifold with f(0) = x.
</p>
<p>
If (y^i) is a local parametrization of the shell at x, and identifying the derivative at x with the gradient vector
</p>
<p>
grad φ = φ_,i e^i where (e^i) is the dual basis corresponding to the parametrization. Thus the gradient vector is tangent to the shell whereas if φ was defined everywhere then the gradient vector would have a component normal to the shell.
</p>
<p>
Practically, this comes up all the time when dealing with finite rotations in FE implementations of shells, while deriving Jacobians.
</p>
<p>
Hope this helps.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 15 May 2007 03:36:12 +0000Amit Acharyacomment 2848 at https://imechanica.orgRe: Tensor derivatives
https://imechanica.org/comment/2841#comment-2841
<a id="comment-2841"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/2824#comment-2824">a slightly easier way</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Amit,
</p>
<p>
Thanks for the pointer. I'll try to write out the derivation when I get the time and post it here.
</p>
<p>
Could you elaborate on why "the set of all invertible second order tensors is not a vector space" and how that affects differentiation? An example of a situation where an excursion leaves the manifold will really be helpful.
</p>
<p>
Biswajit
</p>
<p>
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 14 May 2007 18:56:39 +0000Biswajit Banerjeecomment 2841 at https://imechanica.orgderivative of the inverse
https://imechanica.org/comment/2826#comment-2826
<a id="comment-2826"></a>
<p><em>In reply to <a href="https://imechanica.org/node/1403">Derivatives of the invariants of a tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Of course, the easiest way for the deriv. of the inverse is to work from SS^-1 = I.
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 14 May 2007 11:39:02 +0000Amit Acharyacomment 2826 at https://imechanica.orga slightly easier way
https://imechanica.org/comment/2824#comment-2824
<a id="comment-2824"></a>
<p><em>In reply to <a href="https://imechanica.org/node/1403">Derivatives of the invariants of a tensor</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>
Biswajit,
</p>
<p>
Indeed, algebra and calculus of the type that never loses its charm.
</p>
<p>
Another way: The derivatives of the first and second invariants are easily done by realizing that the trace of a second order tensor is its inner product with the identity. For the second invariant, couple this observation with the chain rule.
</p>
<p>
The third invariant, as well as the derivative of the inverse of a (invertible) tensor is more interesting. I learnt the following from my advisor Don Carlson.
</p>
<p>
Once you know how to do the derivatives of the first and second invariants as above, and then realize that a tensor satisfies its own characteristic equation (Cayley-Hamilton theorem), then taking a derivative of the tensor characteristic equation gives the derivative of the third invariant in terms of the derivatives of the first two invariants, and the derivatives of the tensor itself and its square.
</p>
<p>
Multiply the tensor characteristic equation by the inverse of the tensor and then take a derivative; then, if one knows the derivatives of the invariants, then the derivative of the inverse falls out.
</p>
<p>
The next interesting question is how to do all this if the domain of these functions was a nontrivial manifold, e.g. suppose your s was incompressible (i.e. det s = 1). Of course, even with the inverse there is a bit of an issue to wade through as the set of all invertible second order tensors is not a vector space, but one gets away because it is at least an open set and diferentiation works out as sufficiently small excursions from a point in the domain remains in the domain (a requirement for, say, your definition of the directional derivative to make sense).
</p>
<p>
best,
</p>
<p>
Amit
</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/1403%23comment-form">Log in</a> or <a href="/user/register?destination=node/1403%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 14 May 2007 03:09:29 +0000Amit Acharyacomment 2824 at https://imechanica.org