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
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
