[Yade-users] [Question #697615]: displacement vs friction graph plotting
New question #697615 on Yade: https://answers.launchpad.net/yade/+question/697615 We have written a code for a sphere on a horizontal plate with some contact law between them and a normal load is applied on the sphere. sp=sphere((0,0,0),1) plate=box((0,0,-1.1),(10,10,.1)) O.bodies.append(sp) O.bodies.append(plate) plate.state.vel =(0,.002,0) sp.setPermF(5,(0,0,-1)) mat = FrictMat(frictionAngle=0.5) O.materials.append(mat) InteractionLoop( [Ig2_Box_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ) Can you please suggest the code for plotting graph of "Displacement of plate v/s friction at interface" . Thanks in Advance -- 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 #697612]: Velocity distribution in paraview
New question #697612 on Yade: https://answers.launchpad.net/yade/+question/697612 Hi, When I used paraview to look at the distribution of particle velocity, I met a problem that the direction of all the particle velocity (linVelVec) is always the same, like Figure 1 shows (in the attached link[1]). The way I visualise it in paraview is as follow, which can also be seen in attached Figure 2: Glyph Type = Arrow Scalars = linVelLen Vectors = None Scale mode = scalar The vtk files are also attached in the link. The file is exported from a triaxial compression simulation by using VTKRecorder(fileName='fileName',recorders=['all']). Do you have any comments on why it only shows one direction for all the velocity? Thanks Leonard [1]https://we.tl/t-cRsRDHJJ9c -- 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 #697610]: How to get the acceleration of body
New question #697610 on Yade: https://answers.launchpad.net/yade/+question/697610 Hi I'd like to ask that how to get the acceleration of body in yade? The acceleration I ask here is referring to F=ma. It seems that I can gain the acceleration by F/m, where F is the resultant force which can be gained by O.forces.f(b.id), and m is the mass of the body. I just wonder is there any built-in functions that I can get the acceleration? Thanks! Leonard -- 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 #697600]: Contact (dilation) angle between particles
New question #697600 on Yade: https://answers.launchpad.net/yade/+question/697600 Hi all, I am simulating the direct shear of a granular assemblage with different particle size distributions (PSD) (e.g., normal distribution and power-law distribution). I want to output the contact angle between all the particles to see the distribution of it in the entire granular packing for different PSD models. Is there any way to do that? Thanks in advance! Cheers, Chien-Cheng -- 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 #697593]: How to apply normal load and friction
Question #697593 on Yade changed: https://answers.launchpad.net/yade/+question/697593 Status: Open => Answered Jan Stránský proposed the following answer: > how can I write code for Normal Force one option is O.forces.setPermF [2] > how can I write code for Friction mat = FrictMat(...,frictionAngle=0.5) # [3] O.materials.append(mat) > which contact law to use as I am new to this? For start, use the "most basic" Cundall-Strack law [4]. Once you get familiar, you can add different parameters and/or use different contact models Cheers Jan [2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ForceContainer.setPermF [3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.FrictMat [4] https://yade-dem.org/doc/user.html#basic-dem-model -- 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 #697597]: PeriIsoCompressor.sigma different from getStress()
New question #697597 on Yade: https://answers.launchpad.net/yade/+question/697597 Hello, I am trying to compact a cloud in a periodic cell using PeriIsoCompressor. At first it seems to work well : the cloud gets denser until PeriIsoCompressor.sigma.mean() reaches the target mean stress. But when I compare PeriIsoCompressor.sigma to getStress(), the mean stress is completely different (~26 times lower). Please see below a mwe and its output when using yade at git revision 4852d5601. Am I correctly using PeriIsoCompressor ? Thanks in advance, Sacha ## MWE ## O.periodic = True O.cell.setBox((1, 1, 1)) sp=pack.SpherePack() sp.makeCloud(minCorner=(0,0,0), maxCorner=(1, 1, 1), rMean=5e-3, rRelFuzz=5e-4, num=5000, periodic=True, seed=42) sp.toSimulation() iso_comp = PeriIsoCompressor(charLen=5e-3, stresses=[-1e5], doneHook="O.pause()") O.engines = O.engines + [iso_comp] O.run(wait=True) print("\nPIC mean stress: ", iso_comp.sigma.mean()) print("getStress mean stress: ", getStress().trace()/3) ## Output with yade 4852d5601 ## Welcome to Yade 2021-06-08.git-4852d56 Using python version: 3.6.9 (default, Jan 26 2021, 15:33:00) [GCC 8.4.0] TCP python prompt on localhost:9000, auth cookie `yecssk' XMLRPC info provider on http://localhost:21000 Running script mwe.py PIC mean stress: -99503.28987122637 getStress mean stress: -3824.356053577036 -- 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 #697593]: How to apply normal load and friction
Question #697593 on Yade changed: https://answers.launchpad.net/yade/+question/697593 Status: Needs information => Open Pawan gave more information on the question: I want to apply "Normal Force" on sphere. Can you please suggest how can I write code for Normal Force, Friction, and which contact law to use as I am new to this? Thanks in advance. -- 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 #697548]: Simulating the time-dependent behavior of single grain
Question #697548 on Yade changed: https://answers.launchpad.net/yade/+question/697548 Status: Open => Answered Jan Stránský proposed the following answer: > for days or months I am afraid that then DEM (Yade) is not the tool to use.. Yade is suitable for dynamic simulations, which are "short". It uses explicit time integration, which needs "very short" time step to make the simulation stable. You could maybe use some implicit integration scheme, but I am not sure how it would work for "non-continuous" phenomenons Or some event-driven approach, but then there might be a problem with non-linearity.. Or some combination of above? Some "long-time-step" approach for "calm" periods without breakage, skipping to explicit time integration when breaking occurs until it stabilizes, then again "long-time-step" approach... Are new contacts important to your problem? If not, I would consider other methods then DEM. All in all, sounds like challenging topic :-) > I found a class yade.wrapper.CpmMat. It has "damLaw" for damage evolution in uniaxial tension. I do not know if it is can reach my object. Cpm was developed to simulate short-term behavior of concrete. damLaw parameter just switches between exponential or linear softening branch, so it (neither damLaw neither Cpm) is not suitable for your problem. But, as discussed above, I see the biggest problem in "time length" of your problem (w.r.t. DEM). 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
Re: [Yade-users] [Question #697578]: Import VTK
Question #697578 on Yade changed: https://answers.launchpad.net/yade/+question/697578 Status: Open => Answered Jan Stránský proposed the following answer: Hello, there is no predefined function to import vtk file into Yade. But all needed data are in the file, so in principle it is possible to load it. A better approach is to save not only vtk file, but also other more suitable file for loading. E.g. export.text [1] and ymport.text [2] pair. Cheers Jan [1] https://yade-dem.org/doc/yade.export.html#yade.export.text [2] https://yade-dem.org/doc/yade.ymport.html#yade.ymport.text -- 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 #697336]: How to close the function doneHook in PeriTriaxController
Question #697336 on Yade changed: https://answers.launchpad.net/yade/+question/697336 Jan Stránský posted a new comment: Does triax.doneHook="" what you want or not? 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
Re: [Yade-users] [Question #697593]: How to apply normal load and friction
Question #697593 on Yade changed: https://answers.launchpad.net/yade/+question/697593 Status: Open => Needs information Jan Stránský requested more information: Hello, please provide more information [1]. > I want to apply some load on the sphere What is "some load"? > I want to apply ... friction in their interface. friction is not "applied". Standardly, friction is a material/interaction property, which influences tangential force. But it depends on specific contact law, some do not have friction, then the resulting force may differ among the laws, ... Cheers Jan [1] https://www.yade-dem.org/wiki/Howtoask -- 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