------------------------------------------------------------ revno: 2865 committer: jduriez <jduriez@c1solimara-l> branch nick: yade timestamp: Wed 2011-05-18 12:54:46 +0200 message: Let spheres of SimpleShear be wiser (reducing GSTS->defaultDt) modified: pkg/dem/SimpleShear.cpp pkg/dem/SimpleShear.hpp
-- lp:yade https://code.launchpad.net/~yade-dev/yade/trunk Your team Yade developers is subscribed to branch lp:yade. To unsubscribe from this branch go to https://code.launchpad.net/~yade-dev/yade/trunk/+edit-subscription
=== modified file 'pkg/dem/SimpleShear.cpp' --- pkg/dem/SimpleShear.cpp 2011-05-12 15:13:50 +0000 +++ pkg/dem/SimpleShear.cpp 2011-05-18 10:54:46 +0000 @@ -119,7 +119,6 @@ body = shared_ptr<Body>(new Body); body->groupMask=1; shared_ptr<NormalInelasticMat> mat(new NormalInelasticMat); shared_ptr<Aabb> aabb(new Aabb); -// shared_ptr<SphereModel> gSphere(new SphereModel); shared_ptr<Sphere> iSphere(new Sphere); body->state->pos =position; @@ -138,18 +137,10 @@ aabb->color = Vector3r(0,1,0); - -/* gSphere->radius = radius; - // de quoi avoir des bandes (huit en largeur) de couleur differentes : - gSphere->color = ((int)(Mathr::Floor(8*position.X()/length)))%2?Vector3r(0.7,0.7,0.7):Vector3r(0.45,0.45,0.45); - gSphere->wire = false; - gSphere->shadowCaster = true;*/ - iSphere->radius = radius; - iSphere->color = Vector3r(0.8,0.3,0.3); + iSphere->color = ((int)(floor(8*position.x()/length)))%2?Vector3r(0.7,0.7,0.7):Vector3r(0.45,0.45,0.45);// so that we have eight different colour bands body->shape = iSphere; -// body->geometricalModel = gSphere; body->bound = aabb; } @@ -212,11 +203,11 @@ shared_ptr<GlobalStiffnessTimeStepper> globalStiffnessTimeStepper(new GlobalStiffnessTimeStepper); globalStiffnessTimeStepper->timeStepUpdateInterval = timeStepUpdateInterval; - globalStiffnessTimeStepper->defaultDt=1e-5; + globalStiffnessTimeStepper->defaultDt=3e-6; shared_ptr<KinemCTDEngine> kinemEngine (new KinemCTDEngine); kinemEngine->compSpeed = 10.0; - kinemEngine->targetSigma=1000.0; + kinemEngine->targetSigma=2000.0; shared_ptr<InteractionLoop> ids(new InteractionLoop); === modified file 'pkg/dem/SimpleShear.hpp' --- pkg/dem/SimpleShear.hpp 2011-05-12 15:13:50 +0000 +++ pkg/dem/SimpleShear.hpp 2011-05-18 10:54:46 +0000 @@ -33,7 +33,7 @@ bool generate(std::string& message); - YADE_CLASS_BASE_DOC_ATTRS(SimpleShear,FileGenerator,"Preprocessor for creating a numerical model of a simple shear box.\n\n\t- Boxes (6) constitute the different sides of the box itself\n \t- Spheres are contained in the box. The sample is generated by default via the same method used in TriaxialTest Preprocesor (=> see in source function GenerateCloud). But import of a list of spheres from a text file can be also performed after few changes in the source code.\n \n \t Launching this preprocessor will carry out an oedometric compression, until a value of normal stress equal to 1 MPa (and stable). But with others Engines :yref:`KinemCNDEngine`, :yref:`KinemCNSEngine` and :yref:`KinemCNLEngine`, respectively constant normal displacement, constant normal rigidity and constant normal stress paths can be carried out for such simple shear boxes.\n \n \tNB about micro-parameters : their default values correspond to those used in [Duriez2009a]_ and [Duriez2011]_ to simulate infilled rock joints.", + YADE_CLASS_BASE_DOC_ATTRS(SimpleShear,FileGenerator,"Preprocessor for creating a numerical model of a simple shear box.\n\n\t- Boxes (6) constitute the different sides of the box itself\n \t- Spheres are contained in the box. The sample is generated by default via the same method used in TriaxialTest Preprocesor (=> see in source function GenerateCloud). But import of a list of spheres from a text file can be also performed after few changes in the source code.\n \n \t Launching this preprocessor will carry out an oedometric compression, until a value of normal stress equal to 2 MPa (and stable). But with others Engines :yref:`KinemCNDEngine`, :yref:`KinemCNSEngine` and :yref:`KinemCNLEngine`, respectively constant normal displacement, constant normal rigidity and constant normal stress paths can be carried out for such simple shear boxes.\n \n \tNB about micro-parameters : their default values correspond to those used in [Duriez2009a]_ and [Duriez2011]_ to simulate infilled rock joints.", /* ((string,filename,"../porosite0_44.txt",,"file with the list of spheres centers and radii"))*/ ((Real,thickness,0.001,,"thickness of the boxes constituting the shear box [$m$]")) ((Real,length,0.1,,"initial length (along x-axis) of the shear box [$m$]"))
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
