yade-users team mailing list archive

Thread
Date

about the stress of particle?

To: yade-users <yade-users@xxxxxxxxxxxxxxxxxxx>
From: 焦利青 <jiaoliqing@xxxxxxx>
Date: Thu, 7 Oct 2010 10:06:38 +0800 (CST)

Hi,
Sorry for disturbing you.
In shop.cpp file, there is the code to describe the stress of particles:
Real minRad=(geom->refR1<=0?geom->refR2:(geom->refR2<=0?geom->refR1:min(geom->refR1,geom->refR2)));
        Real crossSection=Mathr::PI*pow(minRad,2);
       
        Vector3r normalStress=((1./crossSection)*geom->normal.dot(phys->normalForce))*geom->normal;
        Vector3r shearStress;
for(int i=0; i<3; i++){
            int ix1=(i+1)%3,ix2=(i+2)%3;
            shearStress[i]=geom->normal[ix1]*phys->shearForce[ix1]+geom->normal[ix2]*phys->shearForce[ix2];
            shearStress[i]/=crossSection;
        }
question:
1. "refR1" means the radius of particle 1, this is the constant, right?
2. how to get the shear stress of particle? where is  "shearStress[i]=geom->normal[ix1]*phys->shearForce[ix1]+geom->normal[ix2]*phys->shearForce[ix2]" from? could you list the mathematical equation?
Thanks.
Best Regards
Liqing