yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #07661
[Branch ~yade-dev/yade/trunk] Rev 2865: Let spheres of SimpleShear be wiser (reducing GSTS->defaultDt)
------------------------------------------------------------
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$]"))
