← Back to team overview

yade-dev team mailing list archive

Re: periodic cell stress tensor

 

By the way, you'll find a few references for Love-Weber definition of the stress tensor in this paper : http://linkinghub.elsevier.com/retrieve/pii/S0020722508001183,
where the older reference, is :
[16] <http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6V32-4XRM86B-3&_user=1531784&_coverDate=12%2F31%2F2009&_rdoc=1&_fmt=high&_orig=browse&_origin=browse&_zone=rslt_list_item&_cdi=5718&_sort=d&_docanchor=&view=c&_ct=1&_refLink=Y&_acct=C000053571&_version=1&_urlVersion=0&_userid=1531784&md5=afb544c838dbe190d98637ddb68e3970&searchtype=a#bbib16> A.E. H Love, A Treatise on the Mathematical Theory of Elasticity, Cambridge University Press, Cambridge (1927).

Bruno

On 27/09/10 16:32, Bruno Chareyre wrote:
Hi Jan,

My two cents about symmetry of stress :

f*branch.transpose()
R*T*n[i]*t[j]

This definition of stress is assuming static equilibrium. In this situation, there is a mathematical proof that the result will be symmetric, even if it there is no obvious symmetry in the formula. IIRC, the proof is based on (1) the decomposition of branch into O1C-O2C, and (2) the fact that the total torque of all forces on one grain is 0. If we are not stricly at static equilibrium, this definition does not apply. DEM is never stricly at equlibrium, so this stress is never exactly symmetric, but it is usually very close.
If you do this :
.5*(f*branch.transpose() + branch*f.transpose())
R*T*.5*(n[i]*t[j]+n[j]*t[i])

which is standard way of tensors symetrization.
you get something symmetric. Ok, it looks better, but you are in fact only hidding the initial mistake : abusing static equlibrium assumption, and the result in itself is not better. I'd keep the non-symmetric definition, and the non symmetry of the result would just remind people to be carefull with that. For dynamic cases, there are other definitions of stress, considering kinetic energy of grains, probably not what you want at the moment.

Cheers.

Bruno

_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp



--
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : +33 4 56 52 86 21
Fax : +33 4 76 82 70 43
________________


References