Re: [Yade-users] [Question #698647]: Specify makeClumpCloud composition
Question #698647 on Yade changed:
https://answers.launchpad.net/yade/+question/698647
Status: Open => Answered
Karol Brzezinski proposed the following answer:
Hi Bence,
the simplest way to do this is to repeat c1 three times and c2 seven
times in the list passed to makeClumpCloud(). However, you won't get
exactly 30%, rather approximately 30%.
Cheers,
Karol
##
from yade import pack
O.bodies.append(geom.facetBox((0,0,0),(0.6,1,0.8),wallMask=55,fixed=True))
sp=pack.SpherePack()
c1=pack.SpherePack([((0,0,0),.05),((0.1,0,0),0.05),((0.2,0,0),.05)])
sp=pack.SpherePack()
c2=pack.SpherePack([((0,0,0),.05),((0.1,0,0),0.05)])
sp.makeClumpCloud((-0.5,-1,-0.8),(0.5,1,0.8),[c1,c1,c1,c2,c2,c2,c2,c2,c2,c2],num=-1)
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
NewtonIntegrator(gravity=(0,-9.81,0),damping=0.2),
]
O.dt=.5*PWaveTimeStep()
O.saveTmp()
c1_count = 0
c2_count = 0
for b in O.bodies:
if isinstance(b.shape,Clump):
if b.state.mass >1.5:
c1_count += 1
else:
c2_count += 1
print("C1/(C1+C2) = {:.2f}".format(c1_count/(c1_count+c2_count)))
--
You received this question notification because your team yade-users is
an answer contact for Yade.
___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #698647]: Specify makeClumpCloud composition
New question #698647 on Yade: https://answers.launchpad.net/yade/+question/698647 Hello Yade mates, I would like to model the soil using makeClumpLoud with different spherePacks at once. My question is is it possible to specify (in this code) the exact percentage of c1 and c2? E.g 30% c1 and 70% c2 Thanks for your advices in advance. Bence My code: from yade import pack O.bodies.append(geom.facetBox((0,0,0),(0.6,1,0.8),wallMask=55,fixed=True)) sp=pack.SpherePack() c1=pack.SpherePack([((0,0,0),.05),((0.1,0,0),0.05),((0.2,0,0),.05)]) sp=pack.SpherePack() c2=pack.SpherePack([((0,0,0),.05),((0.1,0,0),0.05)]) sp.makeClumpCloud((-0.5,-1,-0.8),(0.5,1,0.8),[c1,c2],num=-1) sp.toSimulation() O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), NewtonIntegrator(gravity=(0,-9.81,0),damping=0.2), ] O.dt=.5*PWaveTimeStep() O.saveTmp() -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698572]: Is it possible to constrain the motion of a body along a circle
Question #698572 on Yade changed: https://answers.launchpad.net/yade/+question/698572 Rohit John posted a new comment: Dear Karol and Jan, Thanks for taking time from your work to come up with really creative solutions to my problem. I will test the solutions and let you know what I find. Kind regards, Rohit K. John -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698592]: When using a FlowEngine, is the confining pressure of the triaxial test an effective stress?
Question #698592 on Yade changed: https://answers.launchpad.net/yade/+question/698592 Status: Open => Answered Robert Caulk proposed the following answer: I confirm Jérômes instinct. "Effective stress" depends on perspective (i.e. effective stress on which grain??). Application of an "effective stress" from an external loading with TriaxialStressController would require making generalizations about pore-pressure inside the specimen. Which seems like a futile process for any dynamic flowengine scenario with non homogenous pore pressures. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698487]: Change boundary condition from rigid wall to periodic boundary
Question #698487 on Yade changed: https://answers.launchpad.net/yade/+question/698487 Jérôme Duriez posted a new comment: Hi, The YADE numerical medium is either finite in size (=> non-periodic use of YADE, often with rigid boundaries) or infinite and periodic (=> periodic use of YADE). You can not have it both ways, and a YADE simulation can not switch from one nature to another, I think. (Note that the use of periodic = True in makeCloud, while simulation itself is not periodic, is somewhat suspect to me) -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698592]: When using a FlowEngine, is the confining pressure of the triaxial test an effective stress?
Question #698592 on Yade changed: https://answers.launchpad.net/yade/+question/698592 Jérôme Duriez posted a new comment: I would say it is the total stress (in your [1] example, triax.goal* is constant but strain evolves during the consolidation phase, no ?) but I'm no longer that sure. Maybe the answer is in Catalano PhD manuscript (in a § about sample boundaries treatment ?) / will come from someone else. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698605]: The distribution of contact normal direction and contact normal force are not uniform after isotropic consolidation
Question #698605 on Yade changed: https://answers.launchpad.net/yade/+question/698605 Jan Stránský posted a new comment: > Despite the apparent randomness in variation of contact forces, regular trends emerge when normal and tangential components of interparticle forces are averaged over groups of contacts with similar orientations. ok. Just have in mind the case that in one bin you have one force 9 N, in another bin 3 forces 3 N. Depending on what you really want, plain sum or average gives (of course) totally different results Anyway, in your case, you assign full force to the bins, no matter their relation to the plane. So if the force is perpendicular to the plane, you still assign its full value, even if it sounds logical not to include it at all. Consider assigning just the projection of the force to the plane. Or just select "not too much perpendicular" interactions. Or ... There are really many options of the evaluation. I suggest: - using (trying at least) the same method for both forces and direction - investigate why the 3 forces plot is exactly the same - investigate if wall interactions were used for direction plot (they are used in "standard" utils.plotDirections) - take into account force direction w.r.t. plane - test the method on artificial input with known results For further reasonable discussion, please provide a MWE (final packing and plot evaluation, no simulation). Cheers Jan -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
