# yade-users team mailing list archive

## about the stress of particle?

 Thread Previous • Date Previous • Date Next • Thread Next

```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
```

 Thread Previous • Date Previous • Date Next • Thread Next