contact normal force and stress

Submitted by prasantud on Tue, 10/18/2011 - 13:46

pair/gran/local does not save normal force

Hi everybody,
I have two questions
1. How to get the normal force at the contact, as the pair/gran/local only accounts for the shear force as mentioned by Jean in
node/212

2.The value of contact stress
conact stress sigma(alpha,beta)

sigma(alpha,beta)=1/V*[sum(i) sum(i neq j)*[rij(alpha)*fij(beta)/2]]

fij(beta)=fij(beta)_tangential+fij(beta)_normal

If normal force is neglected, will it give valid result for contact stress.

regards,
Samantaray

ckloss's picture

ckloss | Tue, 10/18/2011 - 16:33

>>1. How to get the normal force at the contact
>>2.The value of contact stress
you will have to calculate that with a post-processing step from compute pair/gran/local data, or modify compute pair/gran/local for this purpose

Christoph

prasantud | Fri, 12/02/2011 - 17:32

ok i can substract the tangential force from the total force computed by pair/local to get contact normal forces.
thanks