Re: [Yade-users] [Question #292911]: How to count number of particles in a fragment, connected by cohesive bonds
Question #292911 on Yade changed: https://answers.launchpad.net/yade/+question/292911 Jan Stránský proposed the following answer: Hi Varun, sorry, I was just stupid, after O.step() the interactions are there :-) adding following functions at the end of the script (after O.step() ) should do what you want. You can use it also during simulation. Note that this approach is probably not very efficient, but should work. In case the inefficiency is too much, a more sophisticated graph theory should give better answer. In my case, I got several "aggregates" composed of 4 particles, but always only 3 of them were connected by cohesively and one was standalone, which corresponded to the result of aggregs() function. cheers Jan PS: used approach is a nice exercise of recursive function call :-) # def addBodyToAggreg(body,aggreg): # auxiliary function, add body [yade.Body instance] and all its neighbors into aggreg (python set instance) if body.id in aggreg: # do nothing if b is already in aggreg ... return aggreg.add(body.id) # ... otherwise add it to aggreg intrs = body.intrs() for i in intrs: # and add also all its neighbors ... if not isinstance(i.phys,CohFrictPhys): # ... but only that connected with CohFrictPhys interactions continue if i.phys.cohesionBroken: # ... but only not yet broken continue i2 = i.id1 if i.id2==body.id else i.id2 # choose the other body of interaction b2 = O.bodies[i2] addBodyToAggreg(b2,aggreg) # and add it and all its neighbors to aggreg def aggregs(): # actual function to detect standalone aggregates ids = set(b.id for b in O.bodies if isinstance(b.shape,Sphere)) # first make a set of all spheres ids ret = [] while len(ids)>0: # while there are still some particles not assigned to any aggregate ... i = ids.pop() # ... choose one random ... b = O.bodies[i] a = set() # ... create new aggregate (set of sphere ids) addBodyToAggreg(b,a) # ... and add the sphere together with all its neigbors to aggregate for bid in a: # delete all used ids from ids ids.discard(bid) ret.append(a) return ret aggs = aggregs() for a in aggs: print a # 2016-05-02 18:27 GMT+02:00 VG <question292...@answers.launchpad.net>: > Question #292911 on Yade changed: > https://answers.launchpad.net/yade/+question/292911 > > VG posted a new comment: > Hello Jan, > > I have defined sample_material as CohFrictMat and also using > > Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,label="cohesiveIp") > in the interaction loop. I am not sure what am I missing here. Could you > please help me identify the problem in my model set up. > > Thanks > Varun > > -- > 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 > -- 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 #292911]: How to count number of particles in a fragment, connected by cohesive bonds
Question #292911 on Yade changed: https://answers.launchpad.net/yade/+question/292911 Jan Stránský proposed the following answer: Hi Varun, the script works fine from the point of view that there are no errors. However, in your question you describe cohesive bonds. If I got your question correctly, your fragment is a set of particles mutually connected with cohesive bonds. After running your script, there are no cohesive bonds. This was my point. Thanks Jan 2016-05-02 1:22 GMT+02:00 VG <question292...@answers.launchpad.net>: > Question #292911 on Yade changed: > https://answers.launchpad.net/yade/+question/292911 > > VG posted a new comment: > Hello Jan, > > Thanks for your response. The script seems to be running fine. Is there > a problem with the interactions in my model set up ? > > Regards > Varun > > -- > 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 > -- 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 #292846]: Particle packing - How to avoid overlap
Question #292846 on Yade changed: https://answers.launchpad.net/yade/+question/292846 Jan Stránský proposed the following answer: In this case, after update of yadedaily, the fix should be there soon cheers Jan 2016-05-02 1:27 GMT+02:00 VG <question292...@answers.launchpad.net>: > Question #292846 on Yade changed: > https://answers.launchpad.net/yade/+question/292846 > > VG posted a new comment: > Hello Jan, > > I installed Yade using the below instructions for pre-built package on > website: https://yade-dem.org/doc/installation.html > > sudo bash -c 'echo "deb http://www.yade-dem.org/packages/ trusty/" >> > /etc/apt/sources.list' > wget -O - http://www.yade-dem.org/packages/yadedev_pub.gpg | sudo apt-key > add - > sudo apt-get update > sudo apt-get install yadedaily > > Thanks > Varun > > -- > 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 > -- 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 #292846]: Particle packing - How to avoid overlap
Question #292846 on Yade changed: https://answers.launchpad.net/yade/+question/292846 Status: Open => Answered Jan Stránský proposed the following answer: Hi Varun, > VG is still having a problem: > Thinking about this problem more, this might not actually be a bug. As you > outlined, the difference between the previous version and "fixed" version > is: > - Previous version: Pack all the smaller spheres within the predicate > (bigger sphere), such that the centers of smaller spheres are within the > boundary of predicate. > - Fixed version: Only the smaller spheres which are entirely contained > within the predicate are considered. > > What the fixed version does can also possibly be attained by > pack.inSphere((s[0][0],s[0][1],s[0][2]),(s[1]-particle_radius)). > Correct me if I am wrong. The description of previous and fixed versions is correct in pack.inSphere, it was a real bug. Of course, you can play with settings like using s[1] or s[1]-particle_radius, but this is different point. > > That still doesn't solve the physical problem I am trying to set up. I > should give a better description of the problem I am trying to solve, if > there is a better set up possible for the problem: > > -- I have a granular sample with given input particle size distribution. > Under mechanical loading, these particles will break and I want to get > the resulting particle size distribution. > > -- I am representing the input particle size distribution with spherical > particles of diameters D1, D2, D3 etc. All of such big particles are > composed of smaller particles of diameter 'd', since d is the smallest > particle size of interest. > > -- Initially, within each big spherical particle, all the smaller > particles should be densely packed and connected with cohesive bonds. As > the loading is applied, the initially defined cohesive bonds within each > big sphere will break progressively leading to smaller fragments, giving > a new particle size distribution. > thanks for clarifying, but your description was OK before :-) how you want to create the initial packing is the matter of your preferences. randomDensePack within a big agregate is a good choice. For the big aggregates, you can use: - makeCloud like now - randomDensePack to get denser packing of bigger aggrgates - use something like randomDensePack on already "meshed" aggregates, see Klaus's answer #3 in [1] > > > The major problems I am facing in the current script are: > > -- Smaller particles from two different predicates (big spheres) tend to > overlap. > this should not happen with newer version > > -- I am noticing that sometimes cohesive bond gets formed between two > smaller particles of different predicates. I would like to avoid this, > since I only want non-cohesive contact between the two different big > spheres. > > -- Initially, within a single big sphere itself, sometimes there are > particles which are shown not to have any cohesive interaction, leaving > that small single particle free to fly. > this is material model dependent. I am not very familiar with CohFrictMat, so somebody else have to help. Probably open a new question on this topic, since it is different from the original problem. cheers Jan [1] https://answers.launchpad.net/yade/+question/292672 -- 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 #293247]: Physic interaction betweem particles and polyhedra
Question #293247 on Yade changed: https://answers.launchpad.net/yade/+question/293247 Status: Open => Answered Jan Stránský proposed the following answer: Hi Jonathan, currently you can have only one instance (thus one set of parameters) of p2_FrictMat_FrictMat_MindlinPhys. If needed, it could be modified to distinguish between sphere-sphere and sphere-box interactions, but it would need some coding and probably would be a bit hacky.. Do you have more boxes in your simulation? If no, I would use box instead of polyhedra. cheers Jan 2016-05-04 16:52 GMT+02:00 Jonathan Pergoli < question293...@answers.launchpad.net>: > New question #293247 on Yade: > https://answers.launchpad.net/yade/+question/293247 > > Hi everybody! > > I'm trying to simulate a box falling on a cylindrical pack of spheres. The > spheres have coefficient of restitution = .1 and the box has a COR = .6. > For the spheres I have used Ip2_FrictMat_FrictMat_MindlinPhys(en=.1,es=.1) > . As regard the box, I would like to use the same but if I write: > ... > > [p2_FrictMat_FrictMat_MindlinPhys(en=.1,es=.1),p2_FrictMat_FrictMat_MindlinPhys(en=.6,es=.6)] > the former properties are overwritten by the latter. > > The box is made using utils.polyhedron(). > > I have also used Ip2_FrictMat_PolyhedraMat_CundallStrack(), but I cannot > set the COR. > > Can someone help me,please? > > Thank you :) > > -- > 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 > -- 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 #293216]: Creating Ellipse or Ellipsoid
Question #293216 on Yade changed: https://answers.launchpad.net/yade/+question/293216 Status: Open => Answered Jan Stránský proposed the following answer: Hi Felipe, > 1-> Creating one ellipse or ellpsoid; > 2-> Simulating bouncing.py changing the sphere to ellipse; > ellipsoids are currently not implemented in Yade. A workaround is to create a clump composed of overlaping spheres with surface approximating an ellipsoid. > > 3-> How is it work the function "pack.inEllipsoid"; > inEllipsoid should work in the same way like other pack.predicates, like pack.inCylinder, pack.inAlignedBox, see [1] > > 4-> Is it possible to create one ellipsoid like one pack.cylinder. I have > created a pack cylinder and right now i would like to change this pack for > one pack in ellipsoid. > did you mean pack.inCylinder? if yes, changing inCylinder to inEllipsoid (and probably the arguments) should be enough.. cheers Jan [1] https://yade-dem.org/doc/user.html#constructive-solid-geometry-csg -- 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 #293247]: Physic interaction betweem particles and polyhedra
Question #293247 on Yade changed: https://answers.launchpad.net/yade/+question/293247 Jan Stránský proposed the following answer: instead of using utils.polyhedron() using utils.box() :-) + changing some Ig2 and Law2 fomr Polyhedra to Box cheers Jan 2016-05-04 18:13 GMT+02:00 Jonathan Pergoli < question293...@answers.launchpad.net>: > Question #293247 on Yade changed: > https://answers.launchpad.net/yade/+question/293247 > > Jonathan Pergoli posted a new comment: > Hi Jan, > > How can be modified? > > I have only one box. > > -- > 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 > -- 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 #292846]: Particle packing - How to avoid overlap
Question #292846 on Yade changed: https://answers.launchpad.net/yade/+question/292846 Jan Stránský proposed the following answer: Hi Varun, I meant to use clumps only for the stage of defining initial packing for your actual simulation. The steps I meant are: 1) use makeCloud to create aggregates 2) use randomDensePack to divide each aggregate into a set of spherical particles 3) apply some compression to make the packing of aggregates denser 1) and 2) you do already, the rest is just to make the packing overall denser cheers Jan 2016-05-04 20:26 GMT+02:00 VG <question292...@answers.launchpad.net>: > Question #292846 on Yade changed: > https://answers.launchpad.net/yade/+question/292846 > > VG posted a new comment: > Thanks for your response, Jan! > For the generation of big aggregates, one of the procedures you suggested > is not very clear to me: > > >"- use something like randomDensePack on already "meshed" aggregates, see > > Klaus's answer #3 in [1]" > > This procedure talks about creating clumps of spherical particles. As > far as I understand, they behave in a rigid manner and probably I can't > have breakable cohesive bonds within a clump. Also, I dint quite get > what you mean by already "meshed" aggregates ? > > I will post another question about the issue with CohFrictMat. > > -- > 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 > -- 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 #293471]: Circular RandomDensePack with CpmMat
Question #293471 on Yade changed: https://answers.launchpad.net/yade/+question/293471 Jan Stránský proposed the following answer: Hi Mohsen, > >Not to create any cohesive interactions, use value less than 1 for > intRadius. > I formerly have check it; i reduced the value of intRadius to 0.1, but > again there were cohesive bonds between particles and the particles did not > FLOW freely. > if you reset the value to 1 after first step, the cohesive interactions are still created until 10th step, which is the default value for cohesiveThresholdIter, see below. > >Furthermore, check [1] and make it something like 1 in Ip2. > Jan there is something confusing! when is changed my code to use > 'cohesiveThresholdIter' i received an strange error: > ... > [Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter(-1))], > NameError: name 'cohesiveThresholdIter' is not defined > the correct syntax is Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter=1) Yade [9]: O.engines[2].physDispatcher.functors[0] > -> [9]: > _ > although my material is Cpm, Yade [9] tells my functor is FrictMat this is probably because of the error during O.engines=..., after that the new values are not set and the default values remain (Ip2_FrictMat_FrictMat_FrictPhys among them). to conlude, after setting intRadius=0.5 and Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter=1), everything should work :-) cheers Jan PS: reading the answer once more, cohesiveThresholdIter=0 should work even with intRadius=1 -- 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 #290947]: Memory leak during updating position of tetrapoly elements
Question #290947 on Yade changed: https://answers.launchpad.net/yade/+question/290947 Status: Open => Answered Jan Stránský proposed the following answer: Hi Mike, correct, please report a bug (such that we do not forget about it). Just now I am not able to say where it comes from.. cheers Jan 2016-04-14 19:02 GMT+02:00 MikeZhao <question290...@answers.launchpad.net>: > New question #290947 on Yade: > https://answers.launchpad.net/yade/+question/290947 > > There exists sever memory leaking problem when using intermediate > variables to update positions of tetraPoly in a loop: > > running the code, the memory will keep leaking: > # > a = [[0,0,0],[0,0,1],[0,1,0],[1,0,0]] > O.bodies.append(yade.utils.tetraPoly(a)) > for i in range(100): > b = [[0,0,0],[0,0,1],[0,1,0],[i+2,0,0]] > O.bodies[0].shape.setVertices(b) > O.bodies[0].state.ori = O.bodies[0].shape.GetOri() > O.bodies[0].state.pos = O.bodies[0].shape.GetCentroid() > O.bodies[0].state.mass = O.bodies[0].mat.density * > O.bodies[0].shape.GetVolume() > O.bodies[0].state.inertia = O.bodies[0].mat.density * > O.bodies[0].shape.GetInertia() > > > -- > 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 > -- 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 #290946]: Several Errors occured while using tetraPoly elements
Question #290946 on Yade changed: https://answers.launchpad.net/yade/+question/290946 Status: Open => Answered Jan Stránský proposed the following answer: Hi Mike, at least some of the errors comes from CGAL, an external library used for polyhedra intersection computation. Could you please provide a script(s) reproducing the problems? thanks Jan 2016-04-14 19:08 GMT+02:00 MikeZhao <question290...@answers.launchpad.net>: > Question #290946 on Yade changed: > https://answers.launchpad.net/yade/+question/290946 > > MikeZhao gave more information on the question: > It seems that if you are using four points which are on the same plane > to update the position of a tetraPoly, you will encounter the first > error. Just found it out coincidently. > > -- > 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 > -- 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 #290950]: Particles crossing the wall
Question #290950 on Yade changed: https://answers.launchpad.net/yade/+question/290950 Status: Open => Answered Jan Stránský proposed the following answer: Hello, please read [1], especially the MWE part (there might be other problems than stiffness) > Can I solve this problem by increasing the Young's modulus for both the > materials and what values should I use here ? > Also, wouldn't increasing the Young's Modulus for particles arbitrarily > lead to unphysical response ? > > Probably you can solve the problem by increasing the modulus of spheres. Depending on your problem and expected result, it might be different material, which would not be desired. Increasing modulus of wall has limited effect. The interaction stiffness is computed as a weighted harmonic mean of the two stiffnesses [2]. In the case of soft sphere (E1) and inifinitely stiff wall it would result in interaction stiffness 2*E1*r1, so anyway given by the lower value. > In addition, I am using O.dt=0.5*PWaveTimeStep() in my simulation. Now, > increasing the Young's modulus would also reduce the required time step. > > correct cheers Jan [1] https://yade-dem.org/wiki/Howtoask [2] https://yade-dem.org/doc/formulation.html#normal-stiffness -- 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 #291326]: Packing and predicate dimension
Question #291326 on Yade changed:
https://answers.launchpad.net/yade/+question/291326
Status: Open => Answered
Jan Stránský proposed the following answer:
Hi Luis,
could you please provide a python script?
thanks
Jan
2016-04-18 3:13 GMT+02:00 Luis Barbosa <question291...@answers.launchpad.net
>:
> New question #291326 on Yade:
> https://answers.launchpad.net/yade/+question/291326
>
> Hello everyone,
>
> In my script I'm using a gts surface to create a random dense pack.
>
> However, the dimensions of the pack are different from gts surface, as
> described by these two warnings:
>
> UserWarning: Packing's dimension
> (Vector3(1.7155834296905184,0.9247992240716244,1.7113872666787984)) doesn't
> fully contain dimension of the predicate (Vector3(2,1,2)).
> if dimP[0]>dimS[0] or dimP[1]>dimS[1] or dimP[2]>dimS[2]:
> warnings.warn("Packing's dimension (%s) doesn't fully contain dimension of
> the predicate (%s)."%(dimS,dimP))
>
> Follows the image of the pack. The blue line is gts surface.
>
> https://drive.google.com/open?id=0B9LhrwAf_vasWTd5cG5MektKOXc
>
> How can I create a pack of spheres with the same predicate of gts?
>
> Thanks,
> Luis
>
> --
> 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
>
--
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 #291114]: virtual time and time step
Question #291114 on Yade changed: https://answers.launchpad.net/yade/+question/291114 Status: Open => Answered Jan Stránský proposed the following answer: Hello, for simple cases it works for me. Could you please post your script? thanks Jan PS: it could be related to https://bugs.launchpad.net/yade/+bug/1544022 2016-04-15 2:27 GMT+02:00 VG <question291...@answers.launchpad.net>: > New question #291114 on Yade: > https://answers.launchpad.net/yade/+question/291114 > > I am new to Yade and this question might sound a little silly. I am > running a Yade simulation and from the controller, when I look at the > virtual time at any given instant, it is not equal to Number of iterations > X delta_t (fixed). What's the reason for this ? > > For example, at a particular instant, here are the three values I see from > the controller: > > virt 001s898m453mu002n > iter #49708131, 1528.0/s > delta_t (fixed) = 4.22594404226e-08 > > > -- > 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 > -- 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 #291326]: Packing and predicate dimension
Question #291326 on Yade changed:
https://answers.launchpad.net/yade/+question/291326
Status: Open => Answered
Jan Stránský proposed the following answer:
Hi Luis,
thanks for the files. For some reason, the "smooth" surface, created by a
wall in triaxial compression, is placed wrongly. There are several
solutions, all parameters of randomDensePack, see docs [1]:
spheresInCell=500 (or another number). It does internally periodic triaxial
compression, overcoming the problem with fixed wall position
cropLayers=3, it simply enlarge the compressed volume
dim, according to the documentation, it could also help
my suggestion is 1), since it also reduces time needed for preparation.
cheers
Jan
[1] https://yade-dem.org/doc/yade.pack.html#yade.pack.randomDensePack
2016-04-19 2:17 GMT+02:00 Luis Barbosa <question291...@answers.launchpad.net
>:
> Question #291326 on Yade changed:
> https://answers.launchpad.net/yade/+question/291326
>
> Status: Answered => Open
>
> Luis Barbosa is still having a problem:
> Hi Jan,
>
> Sure I can.
>
> #!/usr/bin/python
> # -*- coding: utf-8 -*-
> #==
> from yade import pack
> import gts
> from yade import utils
> surf=gts.read(open('tetrahedron.gts'))
> pred=pack.inGtsSurface(surf)
> vol = surf.volume()
> #==
> # Spheres
>
> O.materials.append(JCFpmMat(type=1,young=1e8,poisson=0.3,frictionAngle=radians(30),density=3000,tensileStrength=1e6,cohesion=1e6,jointNormalStiffness=1e7,jointShearStiffness=1e7,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=0.0,label='sphereso'))
> #==
> O.bodies.append(pack.gtsSurface2Facets(surf,color=(0,0,1)))
> sp = pack.randomDensePack(pred,radius=0.08,rRelFuzz=0.01, material =
> 'sphereo',useOBB=False,returnSpherePack=True)
> sp.toSimulation(color=(0.9,0.8,0.6))
> #==
> O.engines=[
> ForceResetter(),
> InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
> InteractionLoop(
> [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
>
> [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
>
> [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=False)]
> ),
>
> GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.2),
> NewtonIntegrator(damping=0.9,gravity=[0,0,-9.81]),
> ]
> ve=utils.getSpheresVolume()
> print vol,ve
> #===
> from yade import qt
> qt.View()
> qt.Controller()
> #
>
> And here follows the gts surface...
> https://drive.google.com/open?id=0B9LhrwAf_vasanZ4bVc4b0FHVUE
>
> Thanks,
> Luis
>
> --
> 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
>
--
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 #291320]: Volume of Complex Geometries
Question #291320 on Yade changed: https://answers.launchpad.net/yade/+question/291320 Status: Open => Answered Jan Stránský proposed the following answer: Hi Luis, gts.Surface (argument of pack.inGtsSurface) has a method volume. Once you have volume of all bodies (the part you already know I guess) and volume of predicate (see below) you can calculate whatever from these values. pred = pack.inGtsSurface(surf) vol = surf.volume() # the same surf cheers Jan 2016-04-17 21:47 GMT+02:00 Luis Barbosa < question291...@answers.launchpad.net>: > New question #291320 on Yade: > https://answers.launchpad.net/yade/+question/291320 > > Hello everyone, > > I'm working with crushable aggregates, which have complex geometries > (non-symmetric), by using pack.inGtsSurface. > > What I'm facing now is how to calculate the porosity of these aggregates, > since utils.porosity and utils.voxelPorosity don't calc the total volume of > gts geometry. > > Is there another way to get this volume and calc porosity? > > Thank you! > Luis > > -- > 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 > -- 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 #291529]: Polyhedra longest edge to vtk
Question #291529 on Yade changed: https://answers.launchpad.net/yade/+question/291529 Status: Open => Answered Jan Stránský proposed the following answer: Hello, I don't understand the question 1 at all, sorry :-) Could you specify it bit? What do you mean by "prescribe", "parameter" and how would you use it in simulation? Why do you want something to prescribe before creation if you import the polyhedra? Thanks for answers. Otherwise it should be easy. cheers Jan 2016-04-20 11:52 GMT+02:00 velimier <question291...@answers.launchpad.net>: > New question #291529 on Yade: > https://answers.launchpad.net/yade/+question/291529 > > Hi, > I am creating simulation by importing voronoi cells to yade as polyhedra. > Along with the computation of the voronoi cells I can determine the > longest edge of each cell/polhedra. > My goal is to export the polhedrons as .vtk and to be able to postprocess > them using their longest edge as parameter. > Questions are: > 1. How to prescribe the longest edge to each polyhedron prior the creation > of polyhedrons as "custom" parameter and use it during the simulation? > 2. How to export the longest edge of each polyhedron to .vtk to be able > to postprocess it in paraview? > > > -- > 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 > -- 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 #291529]: Polyhedra longest edge to vtk
Question #291529 on Yade changed: https://answers.launchpad.net/yade/+question/291529 Jan Stránský proposed the following answer: > > > If I understood you well, this aproach can used for Sphere (or any shape) > as well to create attribute that I need for certain simulation? > exactly :-) 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 #291529]: Polyhedra longest edge to vtk
Question #291529 on Yade changed:
https://answers.launchpad.net/yade/+question/291529
Jan Stránský proposed the following answer:
Hello,
thanks for explanation, now it makes more sense :-) just another
terminology point, in this case you would probably assign attribute or
member, parameter is what you give to a function as input.
to your question:
from yade import export,polyhedra_utils
polyhedrons = [polyhedra_utils.polyhedra(defaultMaterial(),size=(.5,.5,.5))
for i in xrange(3)]
for i,b in enumerate(polyhedrons): b.state.pos = (i,0,0)
O.bodies.append(polyhedrons)
maxEdgeLengths = (1,2,3)
# polyhedrons and maxEdgeLengths would be obtained differently in your case
for poly,maxLen in zip(polyhedrons,maxEdgeLengths):
poly.maxEdgeLength = maxLen # you can assign whatever like this, you only
need not to conflict existing names. E.g. poly.state=1 would give error
vtk = export.VTKExporter('/tmp/p.vtk')
vtk.exportPolyhedra(what=[('maxLen','b.maxEdgeLength')]) # maxLen is
paraview name, b.maxEdgeLength is ours assigned value
cheers
Jan
2016-04-20 12:58 GMT+02:00 velimier <question291...@answers.launchpad.net>:
> Question #291529 on Yade changed:
> https://answers.launchpad.net/yade/+question/291529
>
> velimier posted a new comment:
> Hi Jan,
> maybe "assign" is the better term. I meant to create each polyhedron and
> to asign parameter to it so I can access the parameter during the
> simulation. Or is there option to get max edge length of each polyhedron
> in yade?
>
> --
> 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
>
--
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 #290946]: Several Errors occured while using tetraPoly elements
Question #290946 on Yade changed: https://answers.launchpad.net/yade/+question/290946 Jan Stránský proposed the following answer: See [1] for a reference Jan [1] https://bugs.launchpad.net/yade/+bug/1380103 2016-04-14 23:52 GMT+02:00 Jan Stránský < question290...@answers.launchpad.net>: > Question #290946 on Yade changed: > https://answers.launchpad.net/yade/+question/290946 > > Status: Open => Answered > > Jan Stránský proposed the following answer: > Hi Mike, > at least some of the errors comes from CGAL, an external library used for > polyhedra intersection computation. Could you please provide a script(s) > reproducing the problems? > thanks > Jan > > > 2016-04-14 19:08 GMT+02:00 MikeZhao <question290...@answers.launchpad.net > >: > > > Question #290946 on Yade changed: > > https://answers.launchpad.net/yade/+question/290946 > > > > MikeZhao gave more information on the question: > > It seems that if you are using four points which are on the same plane > > to update the position of a tetraPoly, you will encounter the first > > error. Just found it out coincidently. > > > > -- > > 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 > > > > -- > 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 > -- 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 #290947]: Memory leak during updating position of tetrapoly elements
Question #290947 on Yade changed: https://answers.launchpad.net/yade/+question/290947 Jan Stránský proposed the following answer: @Francois: it would be nice also for other state variables. Anyway, the most significant leaks takes place in O.bodies[0].shape.setVertices(b). Look at mentioned bug, #8 [1], it was solved by implementing setVertices for facet. There you pass 3 instances of Vector3, doing no automatic conversion from tuple/list to std::vector. I will try the same fin this case. Jan [1] https://bugs.launchpad.net/yade/+bug/1041084/comments/8 2016-04-15 8:13 GMT+02:00 Anton Gladky <question290...@answers.launchpad.net >: > Question #290947 on Yade changed: > https://answers.launchpad.net/yade/+question/290947 > > Anton Gladky proposed the following answer: > Indeed, those two bugs should be merged. > > Anton > > > 2016-04-15 8:03 GMT+02:00 Kneib François < > question290...@answers.launchpad.net>: > > Question #290947 on Yade changed: > > https://answers.launchpad.net/yade/+question/290947 > > > > Kneib François proposed the following answer: > > Hi Mike, > > This is a known bug: https://bugs.launchpad.net/yade/+bug/1041084 > > I wrote recently a couple of C++ functions as a temporary workaround for > pos, vel, ori, angVel and color : > https://github.com/yade/trunk/commit/71f93009ee9a1dd348c08c269256020eacb2e46a > > > > If you compile Yade from sources, I can write similar functions for mass > > and inertia if you wish. > > > > François > > > > -- > > 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 > > -- > 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 > -- 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 #290947]: Memory leak during updating position of tetrapoly elements
Question #290947 on Yade changed: https://answers.launchpad.net/yade/+question/290947 Jan Stránský proposed the following answer: In the newest trunk version b.shape.setVertices(v1,v2,v3,v4) should work without memory leaks cheers Jan 2016-04-15 10:23 GMT+02:00 Jan Stránský < question290...@answers.launchpad.net>: > Question #290947 on Yade changed: > https://answers.launchpad.net/yade/+question/290947 > > Jan Stránský proposed the following answer: > @Francois: it would be nice also for other state variables. > > Anyway, the most significant leaks takes place > in O.bodies[0].shape.setVertices(b). > > Look at mentioned bug, #8 [1], it was solved by implementing setVertices > for facet. There you pass 3 instances of Vector3, doing no automatic > conversion from tuple/list to std::vector. I will try the same fin this > case. > > Jan > > [1] https://bugs.launchpad.net/yade/+bug/1041084/comments/8 > > > 2016-04-15 8:13 GMT+02:00 Anton Gladky < > question290...@answers.launchpad.net > >: > > > Question #290947 on Yade changed: > > https://answers.launchpad.net/yade/+question/290947 > > > > Anton Gladky proposed the following answer: > > Indeed, those two bugs should be merged. > > > > Anton > > > > > > 2016-04-15 8:03 GMT+02:00 Kneib François < > > question290...@answers.launchpad.net>: > > > Question #290947 on Yade changed: > > > https://answers.launchpad.net/yade/+question/290947 > > > > > > Kneib François proposed the following answer: > > > Hi Mike, > > > This is a known bug: https://bugs.launchpad.net/yade/+bug/1041084 > > > I wrote recently a couple of C++ functions as a temporary workaround > for > > pos, vel, ori, angVel and color : > > > https://github.com/yade/trunk/commit/71f93009ee9a1dd348c08c269256020eacb2e46a > > > > > > If you compile Yade from sources, I can write similar functions for > mass > > > and inertia if you wish. > > > > > > François > > > > > > -- > > > 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 > > > > -- > > 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 > > > > -- > 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 > -- 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 #291326]: Packing and predicate dimension
Question #291326 on Yade changed: https://answers.launchpad.net/yade/+question/291326 Status: Open => Answered Jan Stránský proposed the following answer: Hi Luis, > "No suitable packing in database found, running PERIODIC compression" > > this is just information for the user, it is OK. > And the pack is not dense, but is similar to a "cloud". > this is strange. Could you please post a final picture? In my case it worked as expected. 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 #291326]: Packing and predicate dimension
Question #291326 on Yade changed: https://answers.launchpad.net/yade/+question/291326 Status: Open => Answered Jan Stránský proposed the following answer: Hi Luis, you can also improve the result with setting noPad=True to pack.inGtsSurface. Computing the volume fraction of spheres, I got 57%. randomDensePack in periodic (optimal) conditions has usually around 62% of volume fraction (with monodisperse packing), so the packing is OK. Could you try it and (dis)confirm that it works or not? Thanks Jan 2016-04-21 14:47 GMT+02:00 Luis Barbosa < question291...@answers.launchpad.net>: > Question #291326 on Yade changed: > https://answers.launchpad.net/yade/+question/291326 > > Status: Answered => Open > > Luis Barbosa is still having a problem: > Hi Jan, > > https://drive.google.com/open?id=0B9LhrwAf_vasN1RvNW9Qa3liZUE > > As you can see the spheres are not overlaping, in some case not even > touch. > > Thanks, > Luis > > -- > 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 > -- 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 #293445]: Polyhedra Elements, Infinite Loop.
Question #293445 on Yade changed: https://answers.launchpad.net/yade/+question/293445 Status: Open => Answered Jan Stránský proposed the following answer: Hi Mike, according to PolyhedraVolumetricLaw (which was renamed to Law2_PolyhedraGeom_PolyhedraPhys_Volumetric), it seems you have an older version of Yade. Since then, there were some changes in polyhedra code. Is it possible for you to try the script with some newer version? cheers Jan 2016-05-08 8:57 GMT+02:00 MikeZhao <question293...@answers.launchpad.net>: > New question #293445 on Yade: > https://answers.launchpad.net/yade/+question/293445 > > I have been working on building finite-discrete element model to simulate > fracture generation, and polyhedra elements in yade are used for DEM. > > After several steps, it seems that an infinite loop occurs when > yade.O.step() is called. I want to trace back to the source of this bug, > but it seems O.step() has called several functions in boost library, and it > is hard for me to understand what codes in Yade are involved for O.step(). > > Could you please give me some advice on locating this bug? It will be very > helpful if someone can tell my what codes in Yade are called for O.step(). > > By the way, the engine I used is as following: > O.engines = [ > ForceResetter(), > InsertionSortCollider([Bo1_Polyhedra_Aabb()]), > InteractionLoop( > [Ig2_Polyhedra_Polyhedra_PolyhedraGeom()], > [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], > [PolyhedraVolumetricLaw()] > ), > NewtonIntegrator(damping=0), > ] > > -- > 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 > -- 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 #293295]: Want to connect all the particles in an aggregate using cohesive bond
Question #293295 on Yade changed: https://answers.launchpad.net/yade/+question/293295 Status: Open => Answered Jan Stránský proposed the following answer: Hi Varun, I see the main problem in how to correctly create cohesive interactions between particles with CohFrictMat. Unfortunately, I am not familiar with this model. You have to wait for somebody experienced with the model :-) or you can search CohFrictMat in previous questions [1].. cheers Jan [1] https://answers.launchpad.net/yade 2016-05-04 23:42 GMT+02:00 VG <question293...@answers.launchpad.net>: > New question #293295 on Yade: > https://answers.launchpad.net/yade/+question/293295 > > Hello everyone, > > Here is a description of the problem I am trying to solve: > > -- I have a granular sample with spherical particles of diameters D1, D2, > D3 etc. All of such big particles (lets call them agglomerates) are > composed of smaller particles of diameter 'd', since d is the smallest > particle size of interest. > > -- Initially, within each big spherical particle, all the smaller > particles should be densely packed and connected with cohesive bonds > (CohFrictMat). As the loading is applied, the initially defined cohesive > bonds within each big sphere will break progressively leading to smaller > fragments. > > Now, the problem I am facing: Initially, within a single big sphere > itself, sometimes there are particles which are shown not to have any > cohesive interaction, leaving the small individual particle free to fly. > How can I alleviate this problem ? > > This is a follow-up question to > https://answers.launchpad.net/yade/+question/292846 > > Here is the script: > > > from yade import pack > > > # > # Set up run > > # > run_name="PeriodicCohesive_MWE" > data_root_dir="." > > > > > # > # Materials > > # > plate_material=FrictMat( > young=200e9 > ,poisson=0.3 > ,density=8000. > ,frictionAngle=radians(30) > ,label='plate_mat') > > O.materials.append(plate_material) > > > sample_material=CohFrictMat( > young=4e9 > ,poisson=0.25 > ,density=1400 > ,frictionAngle=radians(30) > ,normalCohesion=1e8*1.2 > ,shearCohesion=.4e8*1.2 > ,momentRotationLaw=True > ,label='sample_mat') > O.materials.append(sample_material) > > > > > # > # Component dimensions and operating condition > > # > # Granular material dimension > sample_diameter=2e-4 > #sample_diameter=126e-6 > sample_radius=sample_diameter/2.0 > # Sub-particle dimension > particle_diameter=74e-6 > particle_radius=particle_diameter/2. > > > # > # grinding plate dimension > > # > > rotvel=2./3.*pi*(1.5+0.5)*.254 > > > # > # Periodic Geometry > > # > > > # Set up periodic boundary conditions > O.periodic=True > xExt=4*sample_diameter > yExt=3.*sample_diameter*2 #to block the periodicity in y direction > zExt=xExt > xLim=xExt > yLim=yExt/4 > zLim=zExt > O.cell.hSize=Matrix3( > xExt, 0, 0, > 0, yExt, 0, > 0, 0, zExt) > > > length=xExt > height=yExt > width=zExt > > # Top and bottom plate thickness > thickness=0.1*height > > > > > > bottomBoxes = [] > for ix in (0,1,2): > for iz in (0,1,2): > bottomBoxes.append(box( # create 3x3 boxes with 1/3 cell > size > center=(xExt/6.*(1+2*ix),yLim - > thickness/2.0,zExt/6.*(1+2*iz)) > ,extents=(xExt/6.,thickness/2.0,zExt/6.) > ,wire=False > ,material='plate_mat' > )) > > bottom_id,bottom_ids = O.bodies.appendClumped(bottomBoxes) # bottom_id is > the clump id, > > > O.bodies[bottom_id].state.blockedDOFs='xyzXYZ' > > > > > #
Re: [Yade-users] [Question #293445]: Polyhedra Elements, Infinite Loop.
Question #293445 on Yade changed: https://answers.launchpad.net/yade/+question/293445 Status: Open => Answered Jan Stránský proposed the following answer: Hi Mike, could you please provide a complete script which does this problem? Thanks Jan 2016-05-09 7:52 GMT+02:00 MikeZhao <question293...@answers.launchpad.net>: > Question #293445 on Yade changed: > https://answers.launchpad.net/yade/+question/293445 > > MikeZhao gave more information on the question: > The new engine is as following: > O.engines = [ > ForceResetter(), > InsertionSortCollider([Bo1_Polyhedra_Aabb()]), > InteractionLoop( > [Ig2_Polyhedra_Polyhedra_PolyhedraGeom()], > [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], > [Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()] > ), > NewtonIntegrator(damping=0), > ] > > According to my test, this bug happens during Yade.O.step(). > > I have added several cout codes at the begining and the end of several > places. It turns out that this bug doesn't locate in > Ig2_Polyhedra_Polyhedra_PolyhedraGeom::go, > Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys::go or > Law2_PolyhedraGeom_PolyhedraPhys_Volumetric::go. > > Could anybody tell me where this bug possibly locate? Thanks a lot. > > -- > 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 > -- 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 #293468]: Force controller
Question #293468 on Yade changed: https://answers.launchpad.net/yade/+question/293468 Status: Open => Answered Jan Stránský proposed the following answer: Hello Thang, one option is to apply force on one or more particles. To apply force on one body use O.forces.setPermF(bodyId,force) To apply distributed load, you can identify a boundary layer of particles and to each apply certain force such that in total the forces approximates the distributed load. cheers Jan [1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ForceContainer.setPermF 2016-05-09 4:17 GMT+02:00 ntthangvn <question293...@answers.launchpad.net>: > Question #293468 on Yade changed: > https://answers.launchpad.net/yade/+question/293468 > > Description changed to: > hi every YADE users, > > My name is Thang, a PhD student in Monash Uni. I have tried to model a > cantiliver beam support to a load P=1kN > Normally, I use Displacement controller by applying velocity on one > particle or addtional facets, but I dont know how to use force controller > by applying additional load on this beam. Could anyone help me to solve the > problem? > > Thank you so much for your kind help! > > -- > 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 > -- 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 #293471]: Circular RandomDensePack with CpmMat
Question #293471 on Yade changed: https://answers.launchpad.net/yade/+question/293471 Status: Open => Answered Jan Stránský proposed the following answer: Hi Mohsen, randomDensePack creates packing with minimal but existing overlaps. you can check it by for i in O.interactions: print i.geom.penetrationDepth Not to create any cohesive interactions, use value less than 1 for intRadius. Furthermore, check [1] and make it something like 1 in Ip2. cheers Jan [1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Ip2_CpmMat_CpmMat_CpmPhys.cohesiveThresholdIter 2016-05-09 6:17 GMT+02:00 mohsen <question293...@answers.launchpad.net>: > New question #293471 on Yade: > https://answers.launchpad.net/yade/+question/293471 > > In the name of god > Hello all > I want to simulate a gravity deposition. The particles are generated in a > circular predicate. If FrictMat is used, there is no problem and the > particles flow freely on the horizontal wall; But in the case of CpmMat, it > seems that the bonds between particles is made. > With regarding to this point that i have setted intRadius = 1 and there is > no initial overlap, why bonds between particles are made? > > > from yade import pack > > mat=CpmMat(young=1e9,neverDamage=True,epsCrackOnset=1e-2) > O.materials.append(mat) > > ###generating circular predicate and filling with spheres > > pred=pack.inCylinder(centerBottom=(0.5,0.5,0),centerTop=(0.5,0.5,1),radius=.25) > spheres=pack.randomDensePack(pred,spheresInCell=500,radius=.55e-1) > O.bodies.append(spheres) > > O.bodies.append(wall(-0.1,axis=2)) > > intRadius = 1 > > O.engines = [ > ForceResetter(), > > InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius,label='is2aabb'),Bo1_Wall_Aabb()]), > InteractionLoop( > > [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intRadius,label='ss2d3dg'),Ig2_Wall_Sphere_ScGeom()], > [Ip2_CpmMat_CpmMat_CpmPhys()], > [Law2_ScGeom_CpmPhys_Cpm()] > ), > NewtonIntegrator(gravity=(0,0,-10),damping=0.5) > ] > O.dt = 0 > O.step() > is2aabb.aabbEnlargeFactor = 1. > ss2d3dg.interactionDetectionFactor = 1. > > O.dt=.5*PWaveTimeStep() > > > -- > 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 > -- 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 #293468]: Force controller
Question #293468 on Yade changed: https://answers.launchpad.net/yade/+question/293468 Status: Open => Answered Jan Stránský proposed the following answer: Hi Thang, it will be constant, but it is possible to change this value each time step to create linear increase Jan 2016-05-09 14:56 GMT+02:00 ntthangvn <question293...@answers.launchpad.net>: > Question #293468 on Yade changed: > https://answers.launchpad.net/yade/+question/293468 > > Status: Answered => Open > > ntthangvn is still having a problem: > Dear Jan, > > Thank you very much for your quick response. > > As reading from your link O.forces.setPermF(bodyId,force) means "set the > value of permanent force on body", is that this force will increase to > the designed load each time step or it will be constant, Jan? > > many thanks > Thang > > -- > 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 > -- 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 #293235]: Update algorithm of shear displacement
Question #293235 on Yade changed:
https://answers.launchpad.net/yade/+question/293235
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello Jibril,
thanks for checking the documentation and code so much in detail. After
checking the code, the quick answer is that you are right and documentation
and actual implementation do differ.
Another discussion is which one is correct and how much it influence the
result or if the difference is negligible (Bruno?).
Also interesting note, the parts of code and documentation you mentioned
are unchanged since 2010 :-)
cheers
Jan
2016-05-04 11:13 GMT+02:00 Jibril Coulibaly <
question293...@answers.launchpad.net>:
> New question #293235 on Yade:
> https://answers.launchpad.net/yade/+question/293235
>
> Hi everyone,
>
> My questions concern the algorithm used for the shear dispalcement update
> and its implementation is the class ScGeom. I would like to fully
> understand the update procedure but it seems like there are discrepancies
> between the documentation and the code.
>
> In the documentation (
> https://yade-dem.org/doc/formulation.html#fig-shear-2d ), the update is
> said to be "done by perpendicular projection to the plane first (which
> might decrease |\uT|) and adding what corresponds to spatial rotation of
> the interaction instead". The formulae for the corresponding increments
> don't indicate a projection and look more like infinitesimal rotations.
> These 2 rotations, acting of the previous shear dispalcement vector, are
> eventually linearly added, not composed, to form the current shear
> dispalcement vector.
>
> Getting into the corresponding code from the class ScGeom (
> https://github.com/yade/trunk/blob/master/pkg/dem/ScGeom.cpp ), the
> ScGeom::rotate function clearly performs the composition of the two
> approximate rotations, not the addition, about what is defined as the
> orthonormal axis first (line 17) and the twist axis second (line 18).
> The determination of these axes is performed in the ScGeom::precompute
> function. While the definition of the orthonormal axis (line 27)
> corresponds to what the documentation says, the definition of the twist
> angle and axis (lines 28,29) differs from the formula in the documentation
> in that is uses the variable "normal" for the current normal vector n° and
> not "currentNormal".
>
> Regarding the documentation, shouldn't the second rotation of the shear
> dispalcement (around the twist axis) be applied to the shear displacement
> already updated by the first rotation (around the orthonormal axis), as it
> is in the function ScGeom::rotate ? That is :
>
> (\Delta \uT)_2&=-(\prevuT + (\Delta \uT)_1) \times\left(\frac{\Delta t}{2}
> \currn \cdot (\pprev{\vec{\omega}}_1+\pprev{\vec{\omega}}_2)\right) \currn
>
> instead of
>
> (\Delta \uT)_2&=-\prevuT\times\left(\frac{\Delta t}{2} \currn \cdot
> (\pprev{\vec{\omega}}_1+\pprev{\vec{\omega}}_2)\right) \currn
>
> And keeping the summation of the 3 increments as it is.
>
> Regarding the ScGeom::precompute function, shouldn't the twist axis be
> around the current normal vector n° as it is in the documentation ? That is
> replacing the variable "normal" by "currentNormal" in lines (28,29) ? This
> way the rotations would be properly composed :
> - rotation around the orthonormal axis first brings the previous normal
> vector together with the current normal vector.
> - rotation around the twist axis (that is colinear to the current normal)
> then performs the remaining rotation.
>
> I would really appreciate your help and hope I have been understandable
> and clear.
>
> --
> 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
>
--
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 #305894]: track coordinates and size of bond fractures ?
Question #305894 on Yade changed:
https://answers.launchpad.net/yade/+question/305894
Jan Stránský proposed the following answer:
Hi James,
according to docs [1,2], one option is ti use
Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(...,
*recordCracks=True,Key='test'*)
It should not be difficult to add some more information in the source code.
another option (with some modifications also applicable to different
materials) is to use a python function and call it with PyRunner. For this
purpose, set (Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(...,
*neverDamage=True*) [3] to make broken interactions still exist and make
possible to get information form them ("crack width" etc.)
###
def processBrokenInteractions():
print
print 'iter',O.iter
for i in O.interactions:
if i.phys.isBroken:
b1 = O.bodies[i.id1]
b2 = O.bodies[i.id2]
w = i.phys.initD-i.geom.penetrationDepth # not 100% sure
print 'cp',i.geom.contactPoint,'mat',b1.mat,'normal',i.geom.normal,'width',w
# write to file instead of print possible
# extract whatever you want from the interaction i
###
cheers
Jan
[1]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM.recordCracks
[2]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM.Key
[3]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM.neverErase
2016-07-28 0:43 GMT+02:00 James X <question305...@answers.launchpad.net>:
> Question #305894 on Yade changed:
> https://answers.launchpad.net/yade/+question/305894
>
> James X posted a new comment:
> Hi Jan,
>
> Thank you very much. I have included a simplified version of what I am
> trying to do. It is currently a pseudo-2D simulation, though I would
> ideally like to be able to handle both 3D and 2D. What I would hope for
> is to be able to export, at any user-defined timestep, the fracture
> event coordinates, the coordinates of the element pair that broke, and
> the material of the two elements.
>
> However, having just the coordinates of the fracture event themselves
> itself would be great.
>
> Just as a FYI, the first half is just a 2D gravity-settling script to
> generate a 2D mesh. While I would like to do 3D simulations, the
> limitations of my experimental data does not allow me to go 3D at the
> moment.
>
>
> import numpy as np
> import time
>
> from yade import pack
> from yade import geom
> from yade import polyhedra_utils
> from yade import utils
> from yade import linterpolation
> from yade import export
> from yade import ymport
>
> #other random parameters
> start = time.time()
>
> #this function just makes a 2D sphere mesh and exports it as
> .sphere
> def Make_Packed_file(filename, xSize, ySize, radius, rRelFuzz_value,
> nbIter=3000):
> #filename = filename #output a msh file
> sp=pack.SpherePack()
>
> sp.makeCloud(minCorner=(0,0,0),maxCorner=(xSize,ySize*2.5,0),rMean=radius,
> rRelFuzz = 0)
> sp.toSimulation()
> for b in O.bodies:
> b.state.blockedDOFs = 'zXY'
>
> box=O.bodies.append(utils.geom.facetBox((xSize/2,1.25*ySize,0),(xSize/2,1.25*ySize,radius),wallMask=63))
> # the same box of makeCloud
>
> O.engines=[
> ForceResetter(),
>
> InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]),
> InteractionLoop(
>
> [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
> [Ip2_FrictMat_FrictMat_FrictPhys()],
> [Law2_ScGeom_FrictPhys_CundallStrack()]
> ),
> NewtonIntegrator(gravity=(0,-2.0,0),damping=.1),
> PyRunner(iterPeriod=500,initRun=False,command='print("elapsed time
> is: " + str(time.time() - start) ) '),
>
> DomainLimiter(lo=(-0.01,-0.01,-0.01),hi=(xSize,ySize,1.1*radius),iterPeriod=nbIter-1)
> ]
> O.dt=utils.PWaveTimeStep()
> O.stopAtIter=nbIter
> O.run()
> O.wait()
> export.text(filename+'.spheres')
>
>
> ###random parameters
> radius = 3
> smoothContact=True
> jointFrict=radians(20)
> jointDil=radians(0)
> new_scale = 1.0
> xSize = 250;
> ySize = 150;
> rRelFuzz_value = 0;
>
> #actually make this packed mesh
> Make_Packed_file("simplified_packed_2D_mesh", xSize, ySize, radius,
> rRelFuzz_value)
>
> O.reset() #reset everything, just in case
>
>
> #defin the material
> def mat_1(): return
> JCFpmMat(type=1,young=30.0e9,poisson=0.3,frictionAngle=radians(30),density=3000,tensileStrength=1.23e8,cohesion=1.23e8,jointNormalStiffness=1e7,jointShearStiffness=1e7,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=
Re: [Yade-users] [Question #305894]: track coordinates and size of bond fractures ?
Question #305894 on Yade changed: https://answers.launchpad.net/yade/+question/305894 Status: Open => Answered Jan Stránský proposed the following answer: Hi James, I think it should be possible with current version of Yade. For sure it can be make possible with additional implementation. Could you please send a script of current version of such simulation? Thanks Jan 2016-07-27 23:27 GMT+02:00 James X <question305...@answers.launchpad.net>: > New question #305894 on Yade: > https://answers.launchpad.net/yade/+question/305894 > > I'd like to find the coordinates of the center when a bond breaks in a > simulation. For example, in a uniaxial simulation, as the material is > stretched and the cohesive bonds between two spheres break, I'd like to be > able to export the coordinates for a flow simulation in a different > program, ideally with some additional idea of how 'wide' this channel is. > > Does anyone know if this is possible? > > -- > 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 > -- 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 #311636]: USA 1.8.00-4.45-2.7.9.0 BitDEFEnder antivirus S.u.p.p.o.r.t p.h.o.n.e N.u.m.b.e.r BitDEFEnder p.h.o.n.e N.u.m.b.e.r Call.for BitDEFEnder .....UsA 1.8.00-4.45-2.7.9
Question #311636 on Yade changed: https://answers.launchpad.net/yade/+question/311636 Status: Open => Invalid Jan Stránský rejected the question: spam -- 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 #311658]: Talk:SUPPORT+I.8.O.O.4.4.5.2.7.9.O.++++ BitDEFEnder customer service number :De..s.c.r..i.b.e. .I.8.O.O.4.4.5.2.7.9.O. BitDEFEnder .A.N.T..I.V.I.R.U.S. ..T.e.c.h.
Question #311658 on Yade changed: https://answers.launchpad.net/yade/+question/311658 Status: Open => Invalid Jan Stránský rejected the question: spam -- 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 #312121]: aViRA+Tech+SUPPORT+UsA I.8.O.O.4.4.5.2.7.9.O aViRA HeLpLinE phone Number aViRA SUppOrt PhOnE NuMbEr aViRA HeLplIne pHoNe NuMbEr.....Usa
Question #312121 on Yade changed: https://answers.launchpad.net/yade/+question/312121 Status: Open => Invalid Jan Stránský rejected the question: spam -- 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 #312187]: aViRA™ I.8.O.O.4.4.5.2.7.9.O aViRA Antivirus tech support phone number aViRA Antivirus technical Support phone Number ,aViRA Antivirus support phone number
Question #312187 on Yade changed: https://answers.launchpad.net/yade/+question/312187 Status: Open => Invalid Jan Stránský rejected the question: spam -- 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 #303013]: The question about Peri3dController
Question #303013 on Yade changed:
https://answers.launchpad.net/yade/+question/303013
Jan Stránský proposed the following answer:
Hello,
sorry for late answer.
The stresses are pretty similar (the nonzero stresses differ 1 and 2 %).
After its creation, Peri3dController uses some rough estimations at the
beginning (to prescribe strain such that stress values are the prescribed
one), so especially if the material is not elastic any more (like after
4000 steps in your case), it behaves slightly differently than ongoing 8000
simulation, so you get slightly different final stresses.
If you prescribe all deformations (stressMask=0), you get almost the same
results in both cases (O.iter is 8000 in one and 8002 in the other example,
so the simulations are not exactly the same anyway)
If some stresses are prescibed, you get some differences
cheers
Jan
2016-07-24 9:52 GMT+02:00 liukeqi <question303...@answers.launchpad.net>:
> New question #303013 on Yade:
> https://answers.launchpad.net/yade/+question/303013
>
> Hi, everyone
> I want to use the engine Peri3dController to reach the goal state, but I
> suffered some question. I used the identical initial condition, and used
> two ways to reach the same goal of strain, but I got different results. The
> first way is that I used the stressMask=0b011100 and goal=(0,-8e-4,0,0,0,0)
> and nSteps = 8000 to compute once. The second way is that I used the
> stressMask=0b011100 and goal=(0,-4e-4,0,0,0,0) and nSteps = 4000 to compute
> twice by using strain=p3d.strain(imitating the example
> peri3dController_triaxialCompression.py ). I thought that the final goal is
> same, and the rate of applied strain is also same (goal/nSteps), so the
> final stress should be same, too, should not it? But the final stress is
> different! The following is my code, can you give me some suggestions?
> Thank you.
>
> # The initial condition
> #!/usr/bin/python # This is server.py file
> import string
> from yade import plot,qt
> from yade.pack import *
> from yade import pack, plot
>
>
> O.materials.append(CpmMat(young=150e9,frictionAngle=atan(0.8),poisson=.2,sigmaT=500e8,epsCrackOnset=1e-6,relDuctility=8))
>
> initSize=1.2
>
>
> sp=pack.randomPeriPack(radius=.05,initSize=Vector3(initSize,initSize,initSize))
>
> sp.toSimulation()
>
> O.save('/tmp/tt.gz')
>
> # The first way, using goal=-8e-4 and
> nSteps=8000 to compute once
> #!/usr/bin/python # This is server.py file
> import string
> from yade import plot,qt
> from yade.pack import *
> from yade import pack, plot
>
> loadFile='/tmp/tt.gz'
> O.load(loadFile)
>
> O.dt=PWaveTimeStep()/2
>
> plot.plots={'ex':('sx',)}
> def plotAddData():
> plot.addData(
> sx=p3d.stress[1],
> ex=p3d.strain[1],
> )
>
> EnlargeFactor=1.5
> EnlargeFactor=1.0
> O.engines=[
> ForceResetter(),
>
> InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=EnlargeFactor,label='bo1s')]),
> InteractionLoop(
>
> [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=EnlargeFactor,label='ig2ss')],
> [Ip2_CpmMat_CpmMat_CpmPhys()],[Law2_ScGeom_CpmPhys_Cpm()]),
> NewtonIntegrator(),
> Peri3dController(
>
> nSteps=8000,
> # how many time steps the simulation will last
> # after reaching
> nSteps do doneHook action
> doneHook='print
> "Simulation with Peri3dController finished."; O.pause()',
>
> # the prescribed
> path (step,value of stress/strain) can be defined in absolute values
>
> # or in relative
> values
>
> # if the goal
> value is 0, the absolute stress/strain values are always considered (step
> values remain relative)plot.plots={'ex':('sx',)}
>
> # if ##Path is not
> explicitly defined, it is considered as linear function between (0,0) and
> (nSteps,goal)
> # as in yzPath and
> xyPath
> # the relative
> values are really relative (zxPath gives the same - except of the sign from
> goal value - result as yyPath)
>
>
> # variables used
> in the first step
>
Re: [Yade-users] [Question #320069]: Box interacton with gravity
Question #320069 on Yade changed: https://answers.launchpad.net/yade/+question/320069 Status: Open => Answered Jan Stránský proposed the following answer: Hi Bernardo, there might be several reasons. Please provide a MWE (see [1], section 3) thanks Jan [1] https://yade-dem.org/wiki/Howtoask Dne 3. 8. 2016 21:52 napsal uživatel "Bernardo Albuquerque" < question320...@answers.launchpad.net>: > > New question #320069 on Yade: > https://answers.launchpad.net/yade/+question/320069 <https://answers.launchpad.net/yade/+question/320069> > > Hello all, > > I am trying to make an oedometric test but for some reasons I need the upper plane to be a body that interacts with the spheres inside the box. I attempted to create a box object and define its density according to the stress I want to aply in the test. For that i used the simple equation > > density = GoalStress * OedometerCrossSection / (gravity * BoxVolume) > > However, I cant achieve the stress I am seeking. I am calculating this stress with O.forces.f(box#)/OedometerCrossSection. The interesting fact is that even if I change the box density the calculated stress is always the same. I first thought that the explanation to my problem was that box objects are not affected by gravity, but I performed a simple script and eliminated that hypothesis. Anybody knows what it could be? > > -- > 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 -- 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 #313648]: randome dense packing in a cylinder with a wire on its axis-howto?
Question #313648 on Yade changed: https://answers.launchpad.net/yade/+question/313648 Status: Open => Answered Jan Stránský proposed the following answer: Hi Jean-François, as a reference for other users who might have the same problem, could you please write how did yo solve it or what was misinterpreted? Thanks Jan 2016-08-04 12:27 GMT+02:00 Jeff <question313...@answers.launchpad.net>: > Question #313648 on Yade changed: > https://answers.launchpad.net/yade/+question/313648 > > Status: Answered => Open > > Jeff is still having a problem: > Hi Jan > Sorry for the late answer but I was away from my desk > > I explored my problem a bit further and I am happy to report that there is > no trace of bug iassociated to my problem > just a poor use of the code by a beginner > All is just working fine now ! > > thanks for your patience and willingness to answer though even if the > question was naive > > JF > > > Jean-François Leon > Galtenco > > Tel [Main-Global]: +1 (415) 322-0376 > Cell [France]: +33 (651) 168-014 > > On Mon, Aug 1, 2016 at 6:17 PM, Jan Stránský < > question313...@answers.launchpad.net> wrote: > > > Your question #313648 on Yade changed: > > https://answers.launchpad.net/yade/+question/313648 > > > > Status: Open => Answered > > > > Jan Stránský proposed the following answer: > > Hi Jean-François, > > > > In fact I initially tried the difference between 2 cylinders you suggest > > and > > > I then moved to the revolution surface because I found that the > > > inside cylinder was fully filled as soon as its diameter was close to > > the > > > sphere diameter. > > > > > > I am not sure if I got it, could you please send a cript reproducing this > > problem? There might be some bug in the code.. > > > > I thought that the revolution surface approach I posted will give me > > an explicit > > > control of the angular discretization. I found with a few experiment > that > > > it was indeed playing a role on the performance of the algorithm but I > > > could never come close to the wire I need to model... so my post. > > > > > > The problem with inGtsSurface is that it evaluates the sphere > in/ousideness > > only approximately in 6 axis aligned directions [1]. It should not be > > difficult to make it more exact. > > > > I assume when you mention other aproach that you are suggesting gravity > > filling > > > with sphere factory ? something else? > > > > > > gravity deposition is one option, some kind of compression within desired > > volume is another similar option. > > > > Random dense pack first create a packing an then it is "cropped" to the > > geometry, so near boundary the packing is not optimal. You can also use > > randomDensePack(...,cropLayers=3), from certain point of view it improves > > the results [2,3]. Using this parameter and difference of two cylinders, > I > > got (from my point of view :-) OK results. > > > > cheers > > Jan > > > > [1] > > > > > https://yade-dem.org/doc/yade.pack.html#yade._packPredicates.inGtsSurface.__init__ > > [2] https://answers.launchpad.net/yade/+question/291326 > > [3] https://yade-dem.org/doc/yade.pack.html#yade.pack.randomDensePack > > > > > > 2016-08-01 23:42 GMT+02:00 Jeff <question313...@answers.launchpad.net>: > > > > > Question #313648 on Yade changed: > > > https://answers.launchpad.net/yade/+question/313648 > > > > > > Status: Answered => Open > > > > > > Jeff is still having a problem: > > > Hi Jan, > > > Thank you for your answer. > > > > > > in fact I initially tried the difference between 2 cylinders you > suggest > > > and I then moved to the revolution surface because I found that the > > inside > > > cylinder was fully filled as soon as its diameter was close to the > > sphere > > > diameter. > > > I thought that the revolution surface approach I posted will give me > an > > > explicit control of the angular discretization. I found with a few > > > experiment that it was indeed playing a role on the performance of the > > > algorithm but I could never come close to the wire I need to model... > so > > my > > > post. > > > > > > I assume when you mention other aproach that you are suggesting gravity > > > filling with sphere factory ? something else? > > > > > > Thanks anyway > > > JF &
Re: [Yade-users] [Question #313648]: randome dense packing in a cylinder with a wire on its axis-howto?
Hi Jean-François, as a reference for other users who might have the same problem, could you please write how did yo solve it or what was misinterpreted? Thanks Jan 2016-08-04 12:27 GMT+02:00 Jeff <question313...@answers.launchpad.net>: > Question #313648 on Yade changed: > https://answers.launchpad.net/yade/+question/313648 > > Status: Answered => Open > > Jeff is still having a problem: > Hi Jan > Sorry for the late answer but I was away from my desk > > I explored my problem a bit further and I am happy to report that there is > no trace of bug iassociated to my problem > just a poor use of the code by a beginner > All is just working fine now ! > > thanks for your patience and willingness to answer though even if the > question was naive > > JF > > > Jean-François Leon > Galtenco > > Tel [Main-Global]: +1 (415) 322-0376 > Cell [France]: +33 (651) 168-014 > > On Mon, Aug 1, 2016 at 6:17 PM, Jan Stránský < > question313...@answers.launchpad.net> wrote: > > > Your question #313648 on Yade changed: > > https://answers.launchpad.net/yade/+question/313648 > > > > Status: Open => Answered > > > > Jan Stránský proposed the following answer: > > Hi Jean-François, > > > > In fact I initially tried the difference between 2 cylinders you suggest > > and > > > I then moved to the revolution surface because I found that the > > > inside cylinder was fully filled as soon as its diameter was close to > > the > > > sphere diameter. > > > > > > I am not sure if I got it, could you please send a cript reproducing this > > problem? There might be some bug in the code.. > > > > I thought that the revolution surface approach I posted will give me > > an explicit > > > control of the angular discretization. I found with a few experiment > that > > > it was indeed playing a role on the performance of the algorithm but I > > > could never come close to the wire I need to model... so my post. > > > > > > The problem with inGtsSurface is that it evaluates the sphere > in/ousideness > > only approximately in 6 axis aligned directions [1]. It should not be > > difficult to make it more exact. > > > > I assume when you mention other aproach that you are suggesting gravity > > filling > > > with sphere factory ? something else? > > > > > > gravity deposition is one option, some kind of compression within desired > > volume is another similar option. > > > > Random dense pack first create a packing an then it is "cropped" to the > > geometry, so near boundary the packing is not optimal. You can also use > > randomDensePack(...,cropLayers=3), from certain point of view it improves > > the results [2,3]. Using this parameter and difference of two cylinders, > I > > got (from my point of view :-) OK results. > > > > cheers > > Jan > > > > [1] > > > > > https://yade-dem.org/doc/yade.pack.html#yade._packPredicates.inGtsSurface.__init__ > > [2] https://answers.launchpad.net/yade/+question/291326 > > [3] https://yade-dem.org/doc/yade.pack.html#yade.pack.randomDensePack > > > > > > 2016-08-01 23:42 GMT+02:00 Jeff <question313...@answers.launchpad.net>: > > > > > Question #313648 on Yade changed: > > > https://answers.launchpad.net/yade/+question/313648 > > > > > > Status: Answered => Open > > > > > > Jeff is still having a problem: > > > Hi Jan, > > > Thank you for your answer. > > > > > > in fact I initially tried the difference between 2 cylinders you > suggest > > > and I then moved to the revolution surface because I found that the > > inside > > > cylinder was fully filled as soon as its diameter was close to the > > sphere > > > diameter. > > > I thought that the revolution surface approach I posted will give me > an > > > explicit control of the angular discretization. I found with a few > > > experiment that it was indeed playing a role on the performance of the > > > algorithm but I could never come close to the wire I need to model... > so > > my > > > post. > > > > > > I assume when you mention other aproach that you are suggesting gravity > > > filling with sphere factory ? something else? > > > > > > Thanks anyway > > > JF > > > > > > Jean-François Leon > > > Galtenco > > > > > > Tel [Main-Global]: +1 (415) 322-0376 > > >
Re: [Yade-users] [Question #315620]: The parameter calibration of CpmMat
Question #315620 on Yade changed: https://answers.launchpad.net/yade/+question/315620 Status: Open => Answered Jan Stránský proposed the following answer: Hello, But I do not know the reasonable range of the per microscope parameter. Can > you give me suggestions? For a set of reasonable parameters, see e.g. yade/examples/concrete, for instance uniax.py [1]. What is "reasonable range" is difficult to say and depends very much on your material and situation where you want to use the model :-) And I also want to know if any of these parameters (young, sigmaT, > relDuctility, poisson, frictionAngle, epsCrackOnset) have conversion > formula to experiment parameter. unfortunatelly no, you have to use some kind of optimization to find values such that the macroscopic response match your requirements. > For example, does the parameter ”young” equal to macroscopic Young’s > moduli in experiment or not. I do not understand the statement “In our > case, however, we simply run elastic simulation to determine the actual > E/kN ratio (3.40). To obtain desired macroscopic modulus of E∗, the value > of kN is scaled by E∗/E.” young parameter equal normal modulus E of bonds represented as rods (considering their crossection is A=pi*r*r) in formula F = k*u = (EA/L)*u, L being distance of particles and u their mutual displacement (penetrationDepth). The resulting macroscopic modulus depends on how "dense" is your interaction network. However, for the same network, the micro and macro moduli are proportional, it is the meaning of the second part of the paragraph. If you double micro modulus, you will get double macro modulus. So you can set the correct young value from one simulation (knowing the ratio of resulting and desired macro modulus). All the other input parameters does not behave so nicely :-) A standard approach would be firstly determine elastic parameters (young, poisson), which are independent on inelasticity. Then determine the inelastic parameters (using fixed young and poisson). For calibration, also consider using the dimensional analysis mentioned in [2] from your previous message. cheers Jan [1] https://github.com/yade/trunk/blob/master/examples/concrete/uniax.py 2016-08-02 17:57 GMT+02:00 liukeqi <question315...@answers.launchpad.net>: > New question #315620 on Yade: > https://answers.launchpad.net/yade/+question/315620 > > Hello, > > I want to calibrate microscope parameters to simulate the > stress-strain curve of concrete obtained by experiment in the reference[1]. > But I do not know the reasonable range of the per microscope parameter. Can > you give me suggestions? > > And I also want to know if any of these parameters (young, sigmaT, > relDuctility, poisson, frictionAngle, epsCrackOnset) have conversion > formula to experiment parameter. For example, does the parameter ”young” > equal to macroscopic Young’s moduli in experiment or not. I do not > understand the statement “In our case, however, we simply run elastic > simulation to determine the actual E/kN ratio (3.40). To obtain desired > macroscopic modulus of E∗, the value of kN is scaled by E∗/E.” in the page > 52 of reference[2]. If it means that the young parameter in CpmMat do not > equal to the Young in experiment. > > Thank you. > > [1] http://link.springer.com/article/10.1007/BF02476284 > [2] https://tel.archives-ouvertes.fr/tel-00502402/ > > > > -- > 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 > -- 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 #340018]: Simulation of the shearbox
Question #340018 on Yade changed: https://answers.launchpad.net/yade/+question/340018 Status: Open => Answered Jan Stránský proposed the following answer: Hi Felipe, there might be several reasons, please provide a MWE [1]. Additionally please specify: - more in detail what you would like to achieve (e.g. I **assume** constant motion=motion at constant velocity) - how the requirements differs from actual simulation (not important with MWE as anybody can try the provided script). thanks Jan [1] https://yade-dem.org/wiki/Howtoask PS: from the videos made from just a few pictures, it is not possible to see anything :-) 2016-08-12 2:12 GMT+02:00 Felipe <question340...@answers.launchpad.net>: > New question #340018 on Yade: > https://answers.launchpad.net/yade/+question/340018 > > Hello Everybody, > > > I am having some issues about simulation of the shear box. When the > surface was moving down, the motion is not constant. It's sounds like leaps > in motion. > > I have tried to change all my script to build the same program with > different ways. However, the result was the same. I was wondering if anyone > could help me to understand what is happening. > > I am sharing one video which was produced at the ParaView (version 5.0.1) > were simulated on Ubuntu 16.04 LTS. > > https://www.dropbox.com/sh/yu6crs6d4bpu996/AACeT41RsvgpBPmJkcX2gr7ra?dl=0 > > > The best motion will be moving down constantly and when the surface gets > either the correct position or apply stress, the box is going to move > right, like shear motion. > > Do Anyone have any idea to solve this problem? > > Hope getting a solution about this problem. > > Best Regards, > > Felipe > > -- > 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 > -- 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 #344674]: find source code
Question #344674 on Yade changed: https://answers.launchpad.net/yade/+question/344674 Status: Open => Answered Jan Stránský proposed the following answer: Hello, L3Geom is old and unmaintained code, I recommend not to spend time on it Ip2_FrictMat_FrictMat_FrictPhys is located in yade/pkg/dem/FrictPhys.*pp files cheers Jan 2016-08-14 15:22 GMT+02:00 张文卿 <question344...@answers.launchpad.net>: > New question #344674 on Yade: > https://answers.launchpad.net/yade/+question/344674 > > I want to understand how yade works，so I‘m ready to read souce code .But I > cann't find the code of"Ig2_Sphere_Sphere_L3Geom"or"Ip2_FrictMat_FrictMat_ > FrictPhys"or"Law2_L3Geom_FrictPhys_ElPerfPl".Can someone tell me where > they are? > > -- > 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 > -- 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 #341304]: Landing of a satellite
Question #341304 on Yade changed:
https://answers.launchpad.net/yade/+question/341304
Status: Open => Answered
Jan Stránský proposed the following answer:
Hi Jonathan,
not studied your script very much in detail, but you may try:
- decrease number of cores (contrary to what one expects :-). Sometimes
using more cores actually make the simulation take longer time, see [1],
page 470.
- not using polyhedrons, as from its nature it takes much longer time to
evaluate them than spheres
- increase time step (instead of 0.5 use higher factor)
- increase particle size (increasing time step)
- decrease number of particle number (using some tricks with boundary
conditions etc.)
- decrease stiffness or increase mass to increase time step (but this
changes also the physics), perhaps increase mass of just a few smallest
particles if you use some PSD (cannot be determined from your script)
- use timing stats to check if some of PyRunner functions does not use too
much time (not likely since you call them not very often)
Anyway, none of the options I gave decrease the simulation time
drastically..
cheers
Jan
2016-08-12 18:17 GMT+02:00 Jonathan Pergoli <
question341...@answers.launchpad.net>:
> New question #341304 on Yade:
> https://answers.launchpad.net/yade/+question/341304
>
> Hi guys,
>
> I am performing simulation in which a satellite lands on the soil. The
> simulation deal with more than 70.000 particles and take a week to make 10
> seconds of simulations using 10 cores. I'd like to know if it is possible
> to speed it up.
> This is the code:
>
> from yade import utils,ymport,export,plot
> import math as m
>
>
> # Material
> E1=1e+8
> E2=5e+7
> mli=FrictMat(density=643,frictionAngle=0.1489,label="MLI",young=E2)
> MLI=O.materials.append(mli)
> gravel=FrictMat(density=1700,frictionAngle=0.7853,label="gravel",young=E1)
> GRAVEL=O.materials.append(gravel)
>
> # Ground
> s=ymport.textExt('1cm_0g2.txt',format='x_y_z_r')
> sphere=O.bodies.append(s)
>
>
> # Create a vector o spheres to eliminate those that have COG above the
> container
> for i in O.bodies:
> if isinstance(i.shape,Sphere):
> if i.state.pos[2]>.30:
> O.bodies.erase(i.id)
> print i.state.mass
> for i in O.bodies:
> if isinstance(i.shape,Sphere):
> if i.state.pos[2]<0:
> O.bodies.erase(i.id)
> for i in O.bodies:
> if isinstance(i.shape,Sphere):
> x=i.state.pos[0]
> y=i.state.pos[1]
> r=m.sqrt(x**2+y**2)
> if r>.75:
> O.bodies.erase(i.id)
>
> aa=[]
> for i in O.bodies:
> if isinstance(i.shape,Sphere):
> aa.append(i.id)
>
> print len(aa)
>
> # Cylinder
> hc=.30
> c=geom.facetCylinder((0,0,.15),radius=.75,height=hc,
> segmentsNumber=100,wallMask=6,material="gravel")
> O.bodies.append(c)
>
> # SAT
> a=.2774
> b=.2922
> c=.1973
> aa=a/2
> bb=b/2
> cc=c/2
> h=.65
> dist=0
> theta=0
> thetav=0
> v1=(aa,bb,c)
> v2=(aa,-bb,c)
> v3=(-aa,-bb,c)
> v4=(-aa,bb,c)
> v5=(aa,bb,0)
> v6=(aa,-bb,0)
> v7=(-aa,-bb,0)
> v8=(-aa,bb,0)
> V=[v1,v2,v3,v4,v5,v6,v7,v8]
> vz=.19
> R=[[m.cos(theta),0,m.sin(theta)],[0,1,0],[-m.sin(theta),0,m.cos(theta)]]
> v1=(R[0][0]*V[0][0]+R[0][1]*V[0][1]+R[0][2]*V[0][2],R[1][0]*
> V[0][0]+R[1][1]*V[0][1]+R[1][2]*V[0][2],R[2][0]*V[0][0]+R[
> 2][1]*V[0][1]+R[2][2]*V[0][2])
> v2=(R[0][0]*V[1][0]+R[0][1]*V[1][1]+R[0][2]*V[1][2],R[1][0]*
> V[1][0]+R[1][1]*V[1][1]+R[1][2]*V[1][2],R[2][0]*V[1][0]+R[
> 2][1]*V[1][1]+R[2][2]*V[1][2])
> v3=(R[0][0]*V[2][0]+R[0][1]*V[2][1]+R[0][2]*V[2][2],R[1][0]*
> V[2][0]+R[1][1]*V[2][1]+R[1][2]*V[2][2],R[2][0]*V[2][0]+R[
> 2][1]*V[2][1]+R[2][2]*V[2][2])
> v4=(R[0][0]*V[3][0]+R[0][1]*V[3][1]+R[0][2]*V[3][2],R[1][0]*
> V[3][0]+R[1][1]*V[3][1]+R[1][2]*V[3][2],R[2][0]*V[3][0]+R[
> 2][1]*V[3][1]+R[2][2]*V[3][2])
> v5=(R[0][0]*V[4][0]+R[0][1]*V[4][1]+R[0][2]*V[4][2],R[1][0]*
> V[4][0]+R[1][1]*V[4][1]+R[1][2]*V[4][2],R[2][0]*V[4][0]+R[
> 2][1]*V[4][1]+R[2][2]*V[4][2])
> v6=(R[0][0]*V[5][0]+R[0][1]*V[5][1]+R[0][2]*V[5][2],R[1][0]*
> V[5][0]+R[1][1]*V[5][1]+R[1][2]*V[5][2],R[2][0]*V[5][0]+R[
> 2][1]*V[5][1]+R[2][2]*V[5][2])
> v7=(R[0][0]*V[6][0]+R[0][1]*V[6][1]+R[0][2]*V[6][2],R[1][0]*
> V[6][0]+R[1][1]*V[6][1]+R[1][2]*V[6][2],R[2][0]*V[6][0]+R[
> 2][1]*V[6][1]+R[2][2]*V[6][2])
> v8=(R[0][0]*V[7][0]+R[0][1]*V[7][1]+R[0][2]*V[7][2],R[1][0]*
> V[7][0]+R[1][1]*V[7][1]+R[1][2]*V[7][2],R[2][0]*V[7][0]+R[
> 2][1]*V[7][1]+R[2][2]*V[7][2])
> p=utils.polyhedron((v1,v2,v3,v4,v5,v6,v7,v8),fixed=False,
> color=(.6,.45,0),material=&qu
Re: [Yade-users] [Question #326664]: when"make" an error occures
Question #326664 on Yade changed: https://answers.launchpad.net/yade/+question/326664 Status: Open => Answered Jan Stránský proposed the following answer: Hello, please read [1] and add at least the most basic info to your question (version of Yade, operating system etc.). With compiling problems, cmake output might be very useful (showing version of VTK for instance). even if you gave some info in [2], it would not hurt to add at least the link (e.g. for those who read only unanswered questions). We can just guess that the problem is on the same computer. Quite important information, what has changed since [2] or why is the output different? cheers Jan [1] https://yade-dem.org/wiki/Howtoask [2] https://answers.launchpad.net/yade/+question/326649 2016-08-06 12:17 GMT+02:00 sandedadi <question326...@answers.launchpad.net>: > New question #326664 on Yade: > https://answers.launchpad.net/yade/+question/326664 > > [ 72%] Building CXX object CMakeFiles/yade.dir/pkg/dem/VTKRecorder.cpp.o > /home/zwq/yade/trunk/pkg/dem/VTKRecorder.cpp: In member function ‘virtual > void VTKRecorder::action()’: > /home/zwq/yade/trunk/pkg/dem/VTKRecorder.cpp:876:13: error: ‘class > vtkXMLUnstructuredGridWriter’ has no member named ‘SetInput’ > writer->SetInput(spheresUg); > > > -- > 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 > -- 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 #326649]: when "make install" , an error appears
Question #326649 on Yade changed: https://answers.launchpad.net/yade/+question/326649 Status: Open => Answered Jan Stránský proposed the following answer: Hello, what version of Yade do you use? The error seems to come from VTK, do you have to use VTKRecorder? An alternative is to try to compile yade without VTK (cmake ... -DENABLE_VTK=OFF) and use export.VTKExporter instead. cheers Jan 2016-08-06 9:32 GMT+02:00 sandedadi <question326...@answers.launchpad.net>: > New question #326649 on Yade: > https://answers.launchpad.net/yade/+question/326649 > > I have cofigure done myfile but when I input "make install" ,an error > occured.my "cmake file" is here: > -- The C compiler identification is GNU 4.8.4 > -- The CXX compiler identification is GNU 4.8.4 > -- Check for working C compiler: /usr/bin/cc > -- Check for working C compiler: /usr/bin/cc -- works > -- Detecting C compiler ABI info > -- Detecting C compiler ABI info - done > -- Check for working CXX compiler: /usr/bin/c++ > -- Check for working CXX compiler: /usr/bin/c++ -- works > -- Detecting CXX compiler ABI info > -- Detecting CXX compiler ABI info - done > -- Found PythonInterp: /usr/bin/python (found version "2.7.6") > -- Try OpenMP C flag = [-fopenmp] > -- Performing Test OpenMP_FLAG_DETECTED > -- Performing Test OpenMP_FLAG_DETECTED - Success > -- Try OpenMP CXX flag = [-fopenmp] > -- Performing Test OpenMP_FLAG_DETECTED > -- Performing Test OpenMP_FLAG_DETECTED - Success > -- Found OpenMP: -fopenmp > -- Looking for Q_WS_X11 > -- Looking for Q_WS_X11 - found > -- Looking for Q_WS_WIN > -- Looking for Q_WS_WIN - not found > -- Looking for Q_WS_QWS > -- Looking for Q_WS_QWS - not found > -- Looking for Q_WS_MAC > -- Looking for Q_WS_MAC - not found > -- Found Qt4: /usr/bin/qmake (found version "4.8.6") > -- Found PkgConfig: /usr/bin/pkg-config (found version "0.26") > -- Version is set to 1.20.0 > -- Looking for XOpenDisplay in /usr/lib/x86_64-linux-gnu/ > libX11.so;/usr/lib/x86_64-linux-gnu/libXext.so > -- Looking for XOpenDisplay in /usr/lib/x86_64-linux-gnu/ > libX11.so;/usr/lib/x86_64-linux-gnu/libXext.so - found > -- Looking for gethostbyname > -- Looking for gethostbyname - found > -- Looking for connect > -- Looking for connect - found > -- Looking for remove > -- Looking for remove - found > -- Looking for shmat > -- Looking for shmat - found > -- Looking for IceConnectionNumber in ICE > -- Looking for IceConnectionNumber in ICE - found > -- Found X11: /usr/lib/x86_64-linux-gnu/libX11.so > -- Found OpenGL: /usr/lib/x86_64-linux-gnu/libGL.so > -- GTS using gts-config /usr/bin/gts-config > -- Using GTS from /usr > -- Found GL2PS: /usr/lib/libgl2ps.so > -- Found CGAL: /usr/include/CGAL, /usr/lib/libCGAL.so > -- Found NumPy: version "1.8.2" /usr/lib/python2.7/dist- > packages/numpy/core/include > -- Found Loki: /usr/include > -- GCC Version >= 4.8. Adding -ftrack-macro-expansion=0 and -save-temps > -- Boost version: 1.54.0 > -- Found the following Boost libraries: > -- python > -- thread > -- filesystem > -- iostreams > -- regex > -- serialization > -- system > -- date_time > -- Boost_VERSION: 105400 > -- Boost_LIB_VERSION: 1_54 > -- Boost_INCLUDE_DIRS: /usr/include > -- Boost_LIBRARIES: /usr/lib/x86_64-linux-gnu/libb > oost_python.so/usr/lib/x86_64-linux-gnu/libboost_thread.so/ > usr/lib/x86_64-linux-gnu/libboost_filesystem.so/usr/ > lib/x86_64-linux-gnu/libboost_iostreams.so/usr/lib/x86_64- > linux-gnu/libboost_regex.so/usr/lib/x86_64-linux-gnu/ > libboost_serialization.so/usr/lib/x86_64-linux-gnu/libboost_ > system.so/usr/lib/x86_64-linux-gnu/libboost_date_time. > so/usr/lib/x86_64-linux-gnu/libpthread.so > -- Found Eigen3: /usr/local/include/eigen3 (Required is at least version > "2.91.0") > -- Found BZip2: /usr/lib/x86_64-linux-gnu/libbz2.so (found version > "1.0.6") > -- Looking for BZ2_bzCompressInit in /usr/lib/x86_64-linux-gnu/libbz2.so > -- Looking for BZ2_bzCompressInit in /usr/lib/x86_64-linux-gnu/libbz2.so > - found > -- Found ZLIB: /usr/lib/x86_64-linux-gnu/libz.so (found version "1.2.8") > -- Found PythonLibs: /usr/lib/x86_64-linux-gnu/libpython2.7.so (found > version "2.7.6") > -- Found Eigen3, version: 3.2.9 > -- Disable vectorization > -- Found VTK > -- Found OpenMP > -- GTS using gts-config /usr/bin/gts-config > -- Using GTS from /usr > -- checking for one of the modules 'glib-2.0' > -- Found GLib2: glib-2.0 /usr/include/glib-2.0;/usr/ > lib/x86_64-linux-gnu/glib-2.0/include > -- Looking for include file glib/gregex.h > -- Looking for include file glib/gregex.h - not found > -- Found G
Re: [Yade-users] [Question #326654]: make error
Question #326654 on Yade changed: https://answers.launchpad.net/yade/+question/326654 Status: Open => Invalid Jan Stránský changed the question status: Already asked here https://answers.launchpad.net/yade/+question/326649 -- 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 #449020]: Quaternions might not be represented properly in YADE
Question #449020 on Yade changed: https://answers.launchpad.net/yade/+question/449020 Status: Open => Answered Jan Stránský proposed the following answer: Hi Nolan, I would recommend to use a builtin method. If you want to reimplement it to see if you understand it or know how it is computed internally, do it afterwards. So use (I got expected results): # q = O.bodies[..].state.ori axis,angle = q.toAxisAngle() print ... # 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 #452513]: overlapping particle/particle and particle/wall
Question #452513 on Yade changed: https://answers.launchpad.net/yade/+question/452513 Status: Open => Answered Jan Stránský proposed the following answer: Hi Michael, I'm new with Yade Welcome :-) the basic idea of DEM is that inter-particle forces depends on mutual overlap (linearly in the most basic case). So if you apply gravity and the sample is supposed to be in static equilibrium, overlaps are intrinsic part of the situation. One solution is to enlarge particles slightly before simulation and shrink them back to original size afterwards. In this case there will be no overlaps, but the particles would not be touching at the same time (both mutually and also between particle-wall). To make particles just touching the walls, you can shrink the particles individually to fulfill exactly the condition of touching. 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 #443228]: segmentation fault using two kinds of cohefricmats under cohesive triaxial test
Question #443228 on Yade changed: https://answers.launchpad.net/yade/+question/443228 Status: Open => Answered Jan Stránský proposed the following answer: Hi Han, I have tried the your script on Ubuntu 16.04 and got error concerning interactions, but not segmentation fault.. 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 #431222]: change position and orientation of clipping planesin terminal
Question #431222 on Yade changed: https://answers.launchpad.net/yade/+question/431222 Jan Stránský posted a new comment: Hi Alexander, have you tried to set the whole object, renderer.clipPlaneSe3=[(pos,ori),(pos,ori),(pos,ori)] ? 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 #446239]: Yade Installation Process
Question #446239 on Yade changed: https://answers.launchpad.net/yade/+question/446239 Status: Open => Answered Jan Stránský proposed the following answer: Hi Felipe, could you please put the .stl files somewhere to the internet and send us the link so we can test the script completely? Thanks 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 #446239]: Yade Installation Process
Question #446239 on Yade changed: https://answers.launchpad.net/yade/+question/446239 Status: Open => Answered Jan Stránský proposed the following answer: Hi Felipe, thanks for the source codes, I have tried them and had no problem. Do you have any problem with them, or was is a different script? 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 #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Status: Open => Answered Jan Stránský proposed the following answer: Hello, does [1] help? cheers Jan [1] https://yade-dem.org/doc/user.html#rendering-interactions-as-force- chain -- 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 #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Status: Open => Answered Jan Stránský proposed the following answer: Hello, [1] normally works without problems achieving width according to the force.. to help you more, we need more info: - what version of Paraview do you use? - did you actually try the last two points: Set “Vary Radius” to “By Scalar” and Set “Radius” and “Radius Factor” such that the result looks OK ? - also what you save (force magnitude or force vector?) might influence the results, like [2]. Could you please provide a MWE or saved VTK data? cheers Jan [2] https://answers.launchpad.net/yade/+question/441859 -- 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 #473518]: Random dense packing gives low density
Question #473518 on Yade changed: https://answers.launchpad.net/yade/+question/473518 Status: Open => Answered Jan Stránský proposed the following answer: Hi Oskari, In the periodic case, randomDensePack results in packing fraction approx. 0.62, which is pretty good (see [1]). I think in this case the boundary effects are very strong. If you use spheresInCell parameter, first a periodic sample is created and then cropped. For 20x10x10 box, I got 0.52. The limit should be around 0.62 For **just postprocessing** reasons, you can add a box to see actual considered box: O.bodies.append(box((1,.5,.5),(1,.5,.5),wire=True)) For the spheresInCell case, you immediately see that near boundary there are no spheres.. A solution should be not to use spheresInCell, but then I got very strange results.. Depending on your simulation requests, a solution could be to use larger predicate with the risk, that some particles will overlap outside the box. cheers Jan [1] https://arxiv.org/pdf/cond-mat/0003416.pdf -- 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 #478212]: How to visualize displacement feild
Question #478212 on Yade changed:
https://answers.launchpad.net/yade/+question/478212
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello,
e.g. using export.VTKExporter + Paraview:
vtk = export.VTKExporter(...)
vtk.exportSpheres(useRef=True,what=[('dspl','b.state.displ()')])
In paraview:
Glyph, Glyph Type = Arrow, Vectors = dspl, Scale Mode = vector, Scale
Factor = 1
colored by "GlyphVector"
I got arrows from original to current positions, colored according to arrow
length.
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 #482400]: index out of range
Question #482400 on Yade changed: https://answers.launchpad.net/yade/+question/482400 Jan Stránský proposed the following answer: Please put here complete error, system you use and Yade version. Thanks Jan 2017-02-25 8:43 GMT+01:00 sunwei <question482...@answers.launchpad.net>: > Question #482400 on Yade changed: > https://answers.launchpad.net/yade/+question/482400 > > sunwei posted a new comment: > in fact, i had already installed it. so it doesn't work. it still has > the" index out of range" problem. > > -- > 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 > -- 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 #482687]: make doc error
Question #482687 on Yade changed: https://answers.launchpad.net/yade/+question/482687 Status: Open => Answered Jan Stránský proposed the following answer: Hello, have you read yade/doc/sphinx/README ? I create documentation with command(s) cd /path/to/yade/doc/sphinx PYTHONPATH=. yade-trunk -x yadeSphinx.py 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 #478212]: How to visualize displacement feild
Question #478212 on Yade changed: https://answers.launchpad.net/yade/+question/478212 Jan Stránský posted a new comment: refPos is set when the body is created (refPos=pos) in case, you can assign it arbitrary value if wanted/needed: b.state.refPos=whatever cheers Jan 2017-02-24 14:48 GMT+01:00 weijy <question478...@answers.launchpad.net>: > Question #478212 on Yade changed: > https://answers.launchpad.net/yade/+question/478212 > > Status: Open => Solved > > weijy confirmed that the question is solved: > Thanks Jan Stránský, that solved my question. > > -- > 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 > -- 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 #482400]: index out of range
Question #482400 on Yade changed:
https://answers.launchpad.net/yade/+question/482400
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello,
"except ImportError: # argparse not present, print error message
> raise RuntimeError("\n\nPlease install 'python-argparse' package.\n")
exactly as the error says, package python-argparse is missing. On ubuntu run
sudo apt-get install python-argparse
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 #455796]: Local energy dissipation
Question #455796 on Yade changed: https://answers.launchpad.net/yade/+question/455796 Jan Stránský proposed the following answer: https://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
Re: [Yade-users] [Question #455796]: Local energy dissipation
Question #455796 on Yade changed:
https://answers.launchpad.net/yade/+question/455796
Status: Open => Answered
Jan Stránský proposed the following answer:
Hi Jabrane,
Is it possible to have the information about particles and contact (force,
> momentum, velocity, stiffenesses etc...) in a specific area of a numerical
> sample ?
of course:
forcesInCertainArea = [O.forces.f(b.id) for b in O.bodies if b.state.pos[2]
> 1.5] # to get forces from bodies with z coordinate > 1.5
avgForce = sum(forcesInCertainArea,Vector3.Zero) / len(forcesInCertainArea
# to average them
To get better answer about energies, please provide a MWE or specify more
in detail how you compute energy ("globally" does not really says much. Do
you use your own functions in Python? O.energy[...]?
Law2.someKindOfEnergy()? ..).
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 #455796]: Local energy dissipation
Question #455796 on Yade changed: https://answers.launchpad.net/yade/+question/455796 Jan Stránský proposed the following answer: Hi Jabrane, you can (just some brainstorming): 1) hard code some values to the code :-) not very good solution, but if you know it in advance, it would make its job 2) create another member of the Law2, probably being a predicate and there would be another computed value (apart from the global value, also "local" value satisfying the predicate condition would be computed 3) The values of dissipated energy would not be only added to some global value, but stored individually for each interaction. 4) From python, you can somehow store values from previous step and with the values from the current step, you can perform computation according to the same equations, but only at certain location. Reading it again, I think 3) is the most clean and general approach. 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 #452868]: Rotate the axis of the figure generated by post2d.plot
Question #452868 on Yade changed: https://answers.launchpad.net/yade/+question/452868 Status: Open => Answered Jan Stránský proposed the following answer: Hi Weimin, after a quick look, I did not find any solution (although I am not very familiar with matplotlib). A dirty workaround is to create a new simulation with bodies with swapped coordinates and swapped data and plot them: ## from yade import post2d import pylab import random r = random.random sphs = ([sphere((-.5+r(),-.5+r(),-1+2*r()),.05) for i in xrange(20)]) # some random particles O.bodies.append(sphs) flattener = post2d.AxisFlatten(axis=1) extractSomething = lambda b: 2*b.state.pos # extract something extractSomething2 = lambda b: b.data # b.data is just temporary variable to save the swapped data # normal process pylab.figure() post2d.plot(post2d.data(extractSomething,flattener)) # get and swap positions of particles poss = [b.state.pos for b in O.bodies] poss = [Vector3(p[2],p[1],p[0]) for p in poss] # get and swap data data = [extractSomething(b) for b in O.bodies] data = [Vector3(d[2],d[1],d[0]) for d in data] O.switchScene() # not to indluence current simulation # add particles with swapped coordinates O.bodies.append([sphere(p,.05) for p in poss]) # assign them swapped data for b in O.bodies: b.data = data[b.id] pylab.figure() post2d.plot(post2d.data(extractSomething2,flattener)) O.switchScene() # back to the original simulation pylab.show() ## 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 #441859]: force chain width (interaction between polyhedras)
Question #441859 on Yade changed: https://answers.launchpad.net/yade/+question/441859 Status: Open => Answered Jan Stránský proposed the following answer: You need a scalar value for this. If you have saved a full vector, use a Calculator filter with "mag(forceN)" between CellDataToPointData and Tube filters. 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 #441859]: force chain width (interaction between polyhedras)
Question #441859 on Yade changed: https://answers.launchpad.net/yade/+question/441859 Status: Open => Answered Jan Stránský proposed the following answer: Hi Weimin, see [1] cheers Jan [1] https://yade-dem.org/doc/user.html#rendering-interactions-as-force- chain -- 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 #444311]: Gts error while installing yade in ubuntu 14.04
Question #444311 on Yade changed: https://answers.launchpad.net/yade/+question/444311 Status: Open => Answered Jan Stránský proposed the following answer: Hello, what version of Yade do you try to compile? thanks 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 #444311]: Gts error while installing yade in ubuntu 14.04
Question #444311 on Yade changed: https://answers.launchpad.net/yade/+question/444311 Jan Stránský proposed the following answer: Hmm, so 3 years old Yade and almost 3 years old system.. Would it be possible for you to use a newer version, or do you have any special reason to use this version? I am afraid nobody will want to find problems in this case.. alternatively you can install ubuntu 14.04 package, which is also 1.07 version by "sudo apt-get install yade" 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 #355381]: ID Bodies Plot
Question #355381 on Yade changed: https://answers.launchpad.net/yade/+question/355381 Status: Open => Answered Jan Stránský proposed the following answer: Hi Luis, geom.facetBox creates several facets, so it is not a single body and therefore you cannot plot its position.. you can sum the forces from all bodies and define position somehow based on position of individual facets (e.g. average or weighted average o positions of individual facets) f = sum(O.forces.f(id)[2] for id in p2) z1 = sum(O.bodies[id].pos[2] for id in p2) / len(p2) is that what you meant? cheers Jan 2016-08-21 4:32 GMT+02:00 Luis Barbosa <question355...@answers.launchpad.net >: > New question #355381 on Yade: > https://answers.launchpad.net/yade/+question/355381 > > Hi, > > I am defining a plane with name p2: > p2=O.bodies.append(utils.geom.facetBox((0,0,r),(0.01,-0.01, > r),wallMask=32)) > > I am using the id to plot, but it is not convenient to identify in each > simulation. > plot.addData(z1=O.bodies[id].state.pos[2],F=O.forces.f(id)[2]) > > Is there a way to plot the position and forces of this body p2 using a > ids=p2? > > Tks, > Luis > > -- > 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 > -- 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 #497350]: Export a polyhedron into a *.stl or * .gts file
Question #497350 on Yade changed: https://answers.launchpad.net/yade/+question/497350 Status: Open => Answered Jan Stránský proposed the following answer: Hi Tina, you can do it "manually" if you know how to write .stl or .gts format. It is not (yet) implemented inside Yade. 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 #478785]: Export VTK files to txt files
Question #478785 on Yade changed: https://answers.launchpad.net/yade/+question/478785 Jan Stránský proposed the following answer: Hi Jabrane, using PyRunner: ## def saveMyDate(): ... # code creating a file and writing your desired data into it O.engines = [ ... PyRunner(iterPeriod=1,command='saveMyData()'), ] ## 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 #482400]: index out of range
Question #482400 on Yade changed: https://answers.launchpad.net/yade/+question/482400 Jan Stránský proposed the following answer: Thanks for info, empty sys.argv is not standard. A workaround should be argv = sys.argv sys.argv = ['whatever'] from yadeimport import * sys.argv = argv Anyway, why you import Yade into Yade session? :-) (based on #4) 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 #497350]: Export a polyhedron into a *.stl or * .gts file
Question #497350 on Yade changed: https://answers.launchpad.net/yade/+question/497350 Status: Open => Answered Jan Stránský proposed the following answer: Hi Tina, I have tried your code and did not face any problem.. please be more specific on "the *.stl only contains some triangles". The code is sopposed to save surface of the polyhedron as a set of triangles. Are some triangles missing? Or did you expect some other output? I don't use stl myself, but according to [1], "In practice, however, all facets are simple triangles.". Probably it would be possible to save each face as a polygon, not only triangles.. 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 #364275]: sintering
Question #364275 on Yade changed: https://answers.launchpad.net/yade/+question/364275 Status: Open => Answered Jan Stránský proposed the following answer: Hello, please read [1] and be a bit more specific what you expect form the simulation. One approach could be to create a clump from sintered particles. cheers Jan [1] https://yade-dem.org/wiki/Howtoask 2016-08-25 15:52 GMT+02:00 张文卿 <question364...@answers.launchpad.net>: > New question #364275 on Yade: > https://answers.launchpad.net/yade/+question/364275 > > I‘d like to simulate the processs of some spheres sintering。but I don’t > konw how to begin。Is there any examples？thank you very much。 > > -- > 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 > -- 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 #378395]: sleeper on ballast
Question #378395 on Yade changed: https://answers.launchpad.net/yade/+question/378395 Status: Open => Answered Jan Stránský proposed the following answer: Hello, I just tried your script and it worked as expected, the sleeper stopped after contact with ballast and also a force jump was recorded.. what version of Yade do you use? I now used 2016-04-21 (not the newest one..) cheers Jan 2016-09-02 0:47 GMT+02:00 Weimin Song <question378...@answers.launchpad.net> : > New question #378395 on Yade: > https://answers.launchpad.net/yade/+question/378395 > > Dear all, > > I am running a similar simulation to [1], which is about the > sleeper on ballast layer. I have two questions: > > (a) As said in [1], box can be used to model the sleeper, > allowing the behavior as a dynamic body (dynamic=True,fixed=False). I put > the sleeper above the ballast and gave sleeper velocity and set one max > load, however, after the contact between sleeper and ballast, the sleeper > does not move and there is no increase of the load (fz). How should I > modify the code? > > (b) Besides the box, what can be used to model the sleeper > better? > > > #compress.txt > #format x_y_z_r > 0.9470860.8314750.2840210.0477039 > 0.4500340.9523870.1200530.0476128 > 0.6463730.6185150.0409984 0.0409984 > 0.4983260.5631530.0749890.0478509 > 0.3232880.495 0.1649160.0477942 > 0.93579 0.1471240.3980740.053411 > 0.83032 0.5750830.0367456 0.0367457 > 0.4800930.45683 0.2252520.0421577 > 0.4493750.9032460.3974590.0513119 > 0.1389880.1984820.1682990.058746 > 0.3067890.3187810.0649523 0.0649523 > 0.6930170.1574850.0644752 0.0644753 > 0.6496560.87162 0.2650320.0480424 > 0.7491240.2259690.1413520.0408891 > 0.68096 0.2219570.3116080.0474974 > 0.1762160.3861360.3315960.0354314 > 0.4685040.6537220.1212460.0502423 > 0.3124650.1913160.0626370.0626371 > 0.3465620.8213410.0981689 0.0522149 > 0.9515080.90208 0.1483180.0449059 > 0.6179550.6022180.2280090.0496856 > 0.1604210.0357093 0.2737750.0350106 > 0.39686 0.1458470.3035830.0423115 > 0.0546970.8608390.3528740.0443734 > 0.3408060.4568220.3591420.0620779 > 0.8847940.2862950.2276740.0523843 > 0.1293030.12767 0.2418120.0420248 > 0.9297430.4462370.3919470.0619074 > 0.5819250.6323650.3268910.0414924 > 0.6904710.7555760.1522380.047335 > 0.8953280.1870660.1635440.064954 > 0.7947050.4997890.0364907 0.0364908 > 0.0619091 0.8713960.0938363 0.058947 > 0.3722460.3649420.4203680.0350741 > 0.8948590.9507340.3664890.0467399 > 0.74198 0.8895480.0377409 0.0377409 > 0.1746890.4051930.2470110.0508009 > 0.1880150.96039 0.1469720.0396097 > 0.7873660.9396110.3891440.0603188 > 0.3476550.1221370.3945730.0580475 > 0.1588520.2030660.3434610.0357254 > 0.1237260.8142160.1796820.0482266 > 0.6581750.6966860.2268380.0406917 > 0.7657490.3582350.1488250.0404131 > 0.9392010.4662950.2598590.0607987 > 0.9570380.0424501 0.2633380.0421011 > 0.1500030.7862540.3881050.0557299 > 0.7842850.2514940.2886190.0618017 > 0.4163940.5530350.1484 0.0626229 > 0.9635090.6114420.4058870.0356438 > 0.3240060.6445220.3944290.0510247 > 0.5145230.7238830.1546890.0396758 > 0.7922490.3001440.2015330.0375461 > 0.5524670.8614770.2650640.0486583 > 0.6176930.3040280.3143020.055041 > 0.6714660.4300640.1895150.0558252 > 0.1280940.7039440.14021 0.0535805 > 0.3428760.6611060.1445460.0380088 > 0.2227930.4627850.32049
Re: [Yade-users] [Question #376480]: remove the wall and append a box (load-unload-stopunloading)
Question #376480 on Yade changed:
https://answers.launchpad.net/yade/+question/376480
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello,
I don't understand the problem at all, please make it more clear.
Especially the point "However, this step can not be realized".
- does it give some error?
- if no, why it cannot be realized?
- what is the purpose to use box instead of wall?
There is no box in you script and otherwise works ok, so the problem cannot
be tested..
thanks
Jan
2016-09-01 21:02 GMT+02:00 Weimin Song <question376...@answers.launchpad.net
>:
> Question #376480 on Yade changed:
> https://answers.launchpad.net/yade/+question/376480
>
> Weimin Song gave more information on the question:
> Code above has bugs. Below is the modified code:
>
> readParamsFromTable(rMean=.05,rRelFuzz=.3,maxLoad=2e8,minLoad=10)
>
>
> from yade.params.table import *
>
> from yade import pack, plot
>
> aggregate=O.materials.append(FrictMat(density=3e3, young=30e9,poisson=.3,
> frictionAngle=.5,label="aggregate"))
> O.bodies.append(geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=31))
>
>
> sp=pack.SpherePack()
> sp.makeCloud((0,0,0),(1,1,1),rMean=rMean,rRelFuzz=rRelFuzz)
> sp.toSimulation(material=aggregate)
> #sp.toSimulation()
>
> wallMat=O.materials.append(FrictMat(density=4e3, young=500e9,
> poisson=0.1, frictionAngle=0.5, label="wallMat"))
> wall=O.bodies.append(wall((.5,.5,1),axis=2,sense=-1, color=(0,0,1),
> material=wallMat))
>
> O.engines=[
>ForceResetter(),
>InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()
> ,Bo1_Wall_Aabb()]),
>InteractionLoop(
> [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom(),Ig2_
> Wall_Sphere_L3Geom()],
> [Ip2_FrictMat_FrictMat_FrictPhys()],
> [Law2_L3Geom_FrictPhys_ElPerfPl()]
>),
>NewtonIntegrator(gravity=(0,0,-9.81),damping=0.5),
>PyRunner(command='checkUnbalanced()',realPeriod=2,label='checker'),
> ]
> O.dt=.5*PWaveTimeStep()
>
> def checkUnbalanced():
>if O.iter<1: return
>if unbalancedForce()>0.2: return
> # O.bodies.append(wall(max([b.state.pos[2]+b.shape.radius for b in
> O.bodies if isinstance(b.shape,Sphere)]),axis=2,sense=-1,
> material=wallMat))
>global plate# without this line, the plate variable would only
> exist inside this function
>plate=O.bodies[-1] # the last particles is the plate
>plate.state.vel=(0,0,-.4)
>O.engines=O.engines+[PyRunner(command='addPlotData()',iterPeriod=50)]
>checker.command='unloadPlate()'
>
> def unloadPlate():
>if abs(O.forces.f(plate.id)[2])>maxLoad:
> plate.state.vel*=-1
> checker.command='stopUnloading()'
>
> def stopUnloading():
>if abs(O.forces.f(plate.id)[2]) plot.saveDataTxt('1.txt')
> O.pause()
> ##remove the facet and append the box
> O.bodies.erase(wall)
> boxMat=O.materials.append(FrictMat(density=4e3, young=500e9,
> poisson=0.1, frictionAngle=0.5, label="boxMat"))
>
> O.bodies.append(box((.5,.5,1.05),(0.1,0.2,0.05),orientation=Quaternion((1,
> 0, 0), 0), dynamic=True, fixed=False, wire=False, color=(0.5,0.5,0.5),
> highlight=False, material=boxMat))
> O.engines=[
> ForceResetter(),
> # sphere, facet, box
> InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()
> ,Bo1_Box_Aabb()]),
> InteractionLoop(
> # the loading plate is a wall, we need to handle sphere+sphere,
> sphere+facet, sphere+wall
> [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom(),Ig2_
> Box_Sphere_ScGeom()],
> [Ip2_FrictMat_FrictMat_FrictPhys()],
> [Law2_ScGeom_FrictPhys_CundallStrack()]
> ),
> NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4),
> PyRunner(command='checkUnbalanced1()',
> realPeriod=2,label='checker1'),
> ]
>
> O.dt=.05*PWaveTimeStep()
>
> def checkUnbalanced1():
> global box# without this line, the box variable would
> only exist inside this function
> box=O.bodies[-1] # the last particles is the box
> box.state.vel=(0,0,-1)
> O.engines=O.engines+[PyRunner(command='addPlotData()',
> iterPeriod=50)]
> checker.command='unloadBox()'
>
> def unloadBox():
> if abs(O.forces.f(box.id)[2])>maxLoad:
> box.state.vel*=-1
> checker.command='stopUnloading1()'
>
>
> def stopUnloading1():
> if abs(O.forces.f(plate.id)[2]) plot.saveDataTxt('2.txt')
> O.pause()
>
> def addPlotData():
>if not isinstance(O.bodies[-1].shape,Wal
Re: [Yade-users] [Question #381111]: where is penetrationDepth？
Question #38 on Yade changed: https://answers.launchpad.net/yade/+question/38 Status: Open => Answered Jan Stránský proposed the following answer: Hello, penetration depth i computed (if this is what you mean by "defining") in Ig2::go methods. For example, for Ig2_Sphere_Sphere_ScGeom see [1], lines 37-39 cheers Jan [1] https://github.com/yade/trunk/blob/master/pkg/dem/Ig2_Sphere_Sphere_ScGeom.cpp 2016-09-04 14:13 GMT+02:00 张文卿 <question381...@answers.launchpad.net>: > New question #38 on Yade: > https://answers.launchpad.net/yade/+question/38 > > When I read c++ code of ScGeom.hpp,I find this line > YADE_CLASS_BASE_DOC_ATTRS_INIT_COTR_PY(ScGeom, > GenericsSpheresContact,((Real,penetrationDepth，NaN，（Attr:: > noSave|Attr::readonly),...))...). > Where is the code of defining penetrationDepth? > > -- > 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 > -- 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 #383099]: The tutorial Script "3d postprocessing" is giving error
Question #383099 on Yade changed: https://answers.launchpad.net/yade/+question/383099 Status: Open => Answered Jan Stránský proposed the following answer: Hello, I could not reproduced your problem (but found different :-) with newer versions of Yade (1.12 is almost 2 years old). Would it be possible to try a newer version? cheers Jan 2016-09-06 13:32 GMT+02:00 feda <question383...@answers.launchpad.net>: > New question #383099 on Yade: > https://answers.launchpad.net/yade/+question/383099 > > Hello, I was learning how to record a video trying the tutorial's script > (3d postprocessing) but the script was returning the error below; Can > someone please help? I'm using Yade 1.12.0-133-a5eec8b~trusty version. > > Qt has caught an exception thrown from an event handler. Throwing > exceptions from an event handler is not supported in Qt. You must > reimplement QApplication::notify() and catch all exceptions there. > > terminate called after throwing an instance of 'std::runtime_error' > what(): Secondary views not supported > Abandon (core dumped) > > Thanks ! > > > -- > 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 > -- 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 #293468]: Force controller
Question #293468 on Yade changed: https://answers.launchpad.net/yade/+question/293468 Jan Stránský proposed the following answer: Hello, just put the call "O.forces.setPermF(bodyId,force)" inside PyRunner and instead of force use O.time*df, where df is force increment per second 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 #402592]: Plot Over Line with Paraview
Question #402592 on Yade changed: https://answers.launchpad.net/yade/+question/402592 Status: Open => Answered Jan Stránský proposed the following answer: Hi Jabrane, please provide a vtk file for testing. Also consider to ask on Paraview forum where is IMO better chance to get correct answer :-) cheers Jan 2016-09-30 12:08 GMT+02:00 Yor1 <question402...@answers.launchpad.net>: > New question #402592 on Yade: > https://answers.launchpad.net/yade/+question/402592 > > Hello, > > I do a test of excavation and i obtain the stress distribution on the > medium. > I want to trace a curve stress/abscissa of particles (X) on a line. > In the paraview's manual : http://www.paraview.org/Wiki/images/d/d5/ > ParaViewTutorial398.pdf (page 34) > they give these instructions to obtain the curve with filter PlotOverLine : > 1. Open the file disk out ref.ex2, load all variables, (see Exercise > 2.4). > 2. Add the Clip filter to disk out ref.ex2, Uncheck the Show Plane > checkbox in the properties panel, and click (like in > Exercise 2.8). This will make it easier to see and manipulate the line > we are plotting over. > 3. Click on disk out ref.ex2 in the pipeline browser to make that the > active > object. > 4. From the menu bar, select Filters → Data Analysis → Plot Over Line > or apply the Plot Over Line filter using the quick launch (ctrl+space > Win/Linux, alt+space Mac). > 5. Adjust the line so that it goes from the base of the disk straight up > to the top of the mesh using the 3D widget manipulators, the p key > shortcut, or the properties panel parameters. > 6. Once you have your line satisfactorily located, click the button. > > I follow these instructions but i don't obtain the curve. > I don't know if i present well the problem. If you want more informations, > i'm glad to give them. > > I tried also the filter PlotData. The problem is that it takes in account > all the particles. > > Best regards. > Jabrane. > > > -- > 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 > -- 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 #402592]: Plot Over Line with Paraview
Question #402592 on Yade changed: https://answers.launchpad.net/yade/+question/402592 Jan Stránský proposed the following answer: Hi Jabrane, using Dropbox, you have to click "Share" and then "Create a link" and send that link to let people other than you view the files Jan 2016-09-30 15:52 GMT+02:00 Yor1 <question402...@answers.launchpad.net>: > Question #402592 on Yade changed: > https://answers.launchpad.net/yade/+question/402592 > > Yor1 posted a new comment: > Hello Jan, > > You find the vtk file in this link: > https://www.dropbox.com/home?preview=20k_sx60sy11_ > excavation-spheres.55302.vtu > > Best regards. > Jabrane. > > -- > 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 > -- 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 #402592]: Plot Over Line with Paraview
Question #402592 on Yade changed: https://answers.launchpad.net/yade/+question/402592 Jan Stránský proposed the following answer: Hi Jabrane, place the files(s) you are postprocessing somewhere on the internet and send here a link (it is not possible to directly attach files here). Thanks Jan 2016-09-30 15:14 GMT+02:00 Yor1 <question402...@answers.launchpad.net>: > Question #402592 on Yade changed: > https://answers.launchpad.net/yade/+question/402592 > > Yor1 posted a new comment: > Hi Jan ! > > Thank you for the response. > I want to know how can i provide the vtk file. > > Best regards. > Jabrane. > > -- > 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 > -- 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 #402779]: issue about irregular shape generation
Question #402779 on Yade changed: https://answers.launchpad.net/yade/+question/402779 Status: Open => Answered Jan Stránský proposed the following answer: Hello, there are polyhedrons and contact between them implemented in Yade. have a look at examples/polyhedra [1]. Also have a look at publications using this implementation [2,3]. >From its nature, the contact detection and evaluation is more time consuming than sphere-sphere contacts. I am not sure if I got last part of your question, but since Yade is open source, you can redefine or write whatever you want. cheers Jan [1] https://github.com/yade/trunk/tree/master/examples/polyhedra [2] http://www.sciencedirect.com/science/article/pii/S0032591014005191 [3] http://congress.cimne.com/particles2013/proceedings/full/p103.pdf 2016-10-05 11:27 GMT+02:00 张德 <question402...@answers.launchpad.net>: > New question #402779 on Yade: > https://answers.launchpad.net/yade/+question/402779 > > Sir, I am working on a project about granulars with irregular > shapesballasts which used on train track ballast bed. The shape has > great influence on the mechanical behavior of ballast bed. So I want to > simulate the irregular shape of ballast precisely. As common methods using > finite number of spheres bonded up to generate irregular shape granular > like PFC3D, but it still can't accurate describe the contact relation > between ballasts as it still relies on the Hertz contact theory between > ball to ball. > So I have a question that does YADE can generate the polyhedrons instead > of cluster of balls bonded up? > As different from spheres, the contact detection and calculation between > polyhedrons have to redifined. So does YADE have the contact algorithm > between polyhedrons? or new contact algorithm could the be writen into the > procedure of YADE(the algorithm may take the non smooth contact method)? > > > -- > 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 > -- 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 #399298]: Machmaker in CPM model
Question #399298 on Yade changed: https://answers.launchpad.net/yade/+question/399298 Status: Open => Answered Jan Stránský proposed the following answer: Ok, in yadedaily you cannot modify the source code.. could you use yade compiled from source? Then you could implement this by yourself :-) cheers Jan 2016-09-19 12:02 GMT+02:00 Seti <question399...@answers.launchpad.net>: > Question #399298 on Yade changed: > https://answers.launchpad.net/yade/+question/399298 > > Status: Answered => Open > > Seti is still having a problem: > Hi Jan, > > > Thank so much for your reply. I am using "yadedaily" for this code. > > Yade: > Yade 1.20.0 > > Yadedaily : > Yade 2016.06a-20-d4840e5~xenial > > -- > 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 > -- 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 #399298]: Machmaker in CPM model
Question #399298 on Yade changed: https://answers.launchpad.net/yade/+question/399298 Status: Open => Answered Jan Stránský proposed the following answer: Hello, multiple CPM materials was not the issue at the time of its implementation, so currently there is MatchMaker implementation only for 'young' parameter. But it is not difficult to implement it for other parameters. What Yade (from source or packages or yadedaily) do you use? cheers Jan 2016-09-18 11:57 GMT+02:00 Seti <question399...@answers.launchpad.net>: > Question #399298 on Yade changed: > https://answers.launchpad.net/yade/+question/399298 > > Seti gave more information on the question: > I have tried to add below lines to CPM Model, and I faced with error as > per below: > > [Ip2_CpmMat_CpmMat_CpmPhys(sigmaT=MatchMaker(matches=(( > concreteId1,concreteId2,3.5e3),(concreteId2,concreteId2,3. > 5e6))),relDuctility=MatchMaker(matches=((concreteId1,concreteId2,30),( > concreteId2,concreteId2,40], > > Error: > > line 101, in > > AttributeError: No such attribute: sigmaT. > > > Would you please advise? > > Regards, > Seti > > -- > 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 > -- 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 #401682]: Problem with angular velocity
Question #401682 on Yade changed: https://answers.launchpad.net/yade/+question/401682 Status: Open => Answered Jan Stránský proposed the following answer: Hi Jonathan, please provide a working script, i.e. with definitions of v1-v8, MLI etc. such that we can try what does not work. do you use the same versions of yadedaily on both machines? do you use the same script on both machines (not like polyhedron on one and sphere on the other)? meanwhile try to set angMom instead of angVel for non-spherical shapes. cheers Jan 2016-09-21 17:53 GMT+02:00 Jonathan Pergoli < question401...@answers.launchpad.net>: > New question #401682 on Yade: > https://answers.launchpad.net/yade/+question/401682 > > Hi everybody! > > I'm trying to set an initial value of the angular velocity for a body. I > set the angular velocity: > body.state.angVel=(0,0,0.1) > > Then when I start yadedaily at the beginning of the simulation the angular > velocity is correct but when I write O.run() I don't know why the angular > velocity is set to (0,0,0). This happen when I use yadedaily for ubuntu > server 14.04 LTS whereas when I run the simulation with ubuntu 14.04 on my > laptop it works well and this problem is not present. Any solutions? > > Thank you in advance, > > Jonathan > > -- > 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 > -- 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 #399298]: Machmaker in CPM model
Question #399298 on Yade changed:
https://answers.launchpad.net/yade/+question/399298
Jan Stránský proposed the following answer:
Hello,
could you please send me an email what parameters you want to make
MatchMaker? I will modify the files for you, but I will not apply them in
the official version (since it would need some testing that the default
behavior does not change etc.)
cheers
Jan
2016-09-19 18:17 GMT+02:00 Seti <question399...@answers.launchpad.net>:
> Question #399298 on Yade changed:
> https://answers.launchpad.net/yade/+question/399298
>
> Status: Answered => Open
>
> Seti is still having a problem:
> Thanks Jan for kind support. I have downloaded yade from
> "https://github.com/yade/trunk;
>
> is that what you have recommended?
>
> if yes,
> 1- How I can install it?
> installing package was really easy just copy/past the mentioned commend
> lines from "https://yade-dem.org/doc/installation.html; however I am not
> sure how should I install the source codes. would you please assist me ?
>
> 2- I have found the source code of CPM on truck and I suppose I should
> change below lines to add matchmaker to control the cohesion. so dose it
> mean if I do changes in " trunk/pkg/dem/ConcretePM.hpp" and save it
> then run YADE it will work?
>
> ###
> class Ip2_CpmMat_CpmMat_CpmPhys: public IPhysFunctor{
> public:
> virtual void go(const shared_ptr& pp1, const
> shared_ptr& pp2, const shared_ptr& interaction);
> FUNCTOR2D(CpmMat,CpmMat);
> DECLARE_LOGGER;
>
> YADE_CLASS_BASE_DOC_ATTRS(Ip2_CpmMat_CpmMat_CpmPhys,IPhysFunctor,"Convert
> 2 :yref:`CpmMat` instances to :yref:`CpmPhys` with corresponding
> parameters. Uses simple (arithmetic) averages if material are different.
> Simple copy of parameters is performed if the :yref:`material` is
> shared between both particles. See :yref:`cpm-model` for detals.",
> ((long,cohesiveThresholdIter,10,,"Should new
> contacts be cohesive? They will before this iter#, they will not be
> afterwards. If 0, they will never be. If negative, they will always be
> created as cohesive (10 by default)."))
> ((shared_ptr,E,,,"Instance of
> :yref:`MatchMaker` determining how to compute interaction's normal modulus.
> If ``None``, average value is used."))
> );
> };
>
> ##
>
> Thanks
>
> --
> 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
>
--
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 #401682]: Problem with angular velocity
Question #401682 on Yade changed: https://answers.launchpad.net/yade/+question/401682 Jan Stránský proposed the following answer: Hi Jonathan, next time please try the script that it really works ;-) and please doulble check if both scripts are exactly the same for the two versions (especially O.bodies.append, see below) To get some effect, you should use O.bodies.append(body) so that the engines know about the particle. Do you use it in your scripts or not? If you don't add body to O.bodies, than angVel does not change (but the body does not move at the same time, because Yade engines do not know about the particle) With O.bodies.append(body), angVel is reset after O.step(), which is expected (see e.g. [1] and links from there). Setting angMom is the solution. The relation is angMom = inertiaTensor*angVel (similar to momentum=mass*velocity). If this works for you, it should not be a problem to correctly define inertiaTensor. cheers Jan [1] https://answers.launchpad.net/yade/+question/274324 2016-09-22 12:57 GMT+02:00 Jonathan Pergoli < question401...@answers.launchpad.net>: > Question #401682 on Yade changed: > https://answers.launchpad.net/yade/+question/401682 > > Jonathan Pergoli posted a new comment: > Hi Jan, > > sorry I will provide the working script: > > E1=1e+8 > E2=5e+7 > mat=FrictMat(density=643,frictionAngle=0.1489,label="MLI",young=E2) > MLI=O.materials.append(mat) > a=.2754 > b=.2822 > c=.1963 > aa=a/2 > bb=b/2 > cc=c/2 > h=.65 > dist=0 > theta=0 > thetav=0 > v1=(aa,bb,c) > v2=(aa,-bb,c) > v3=(-aa,-bb,c) > v4=(-aa,bb,c) > v5=(aa,bb,0) > v6=(aa,-bb,0) > v7=(-aa,-bb,0) > v8=(-aa,bb,0) > V=[v1,v2,v3,v4,v5,v6,v7,v8] > vz=.19 > R=[[m.cos(theta),0,m.sin(theta)],[0,1,0],[-m.sin(theta),0,m.cos(theta)]] > v1=(R[0][0]*V[0][0]+R[0][1]*V[0][1]+R[0][2]*V[0][2],R[1][0]* > V[0][0]+R[1][1]*V[0][1]+R[1][2]*V[0][2],R[2][0]*V[0][0]+R[ > 2][1]*V[0][1]+R[2][2]*V[0][2]) > v2=(R[0][0]*V[1][0]+R[0][1]*V[1][1]+R[0][2]*V[1][2],R[1][0]* > V[1][0]+R[1][1]*V[1][1]+R[1][2]*V[1][2],R[2][0]*V[1][0]+R[ > 2][1]*V[1][1]+R[2][2]*V[1][2]) > v3=(R[0][0]*V[2][0]+R[0][1]*V[2][1]+R[0][2]*V[2][2],R[1][0]* > V[2][0]+R[1][1]*V[2][1]+R[1][2]*V[2][2],R[2][0]*V[2][0]+R[ > 2][1]*V[2][1]+R[2][2]*V[2][2]) > v4=(R[0][0]*V[3][0]+R[0][1]*V[3][1]+R[0][2]*V[3][2],R[1][0]* > V[3][0]+R[1][1]*V[3][1]+R[1][2]*V[3][2],R[2][0]*V[3][0]+R[ > 2][1]*V[3][1]+R[2][2]*V[3][2]) > v5=(R[0][0]*V[4][0]+R[0][1]*V[4][1]+R[0][2]*V[4][2],R[1][0]* > V[4][0]+R[1][1]*V[4][1]+R[1][2]*V[4][2],R[2][0]*V[4][0]+R[ > 2][1]*V[4][1]+R[2][2]*V[4][2]) > v6=(R[0][0]*V[5][0]+R[0][1]*V[5][1]+R[0][2]*V[5][2],R[1][0]* > V[5][0]+R[1][1]*V[5][1]+R[1][2]*V[5][2],R[2][0]*V[5][0]+R[ > 2][1]*V[5][1]+R[2][2]*V[5][2]) > v7=(R[0][0]*V[6][0]+R[0][1]*V[6][1]+R[0][2]*V[6][2],R[1][0]* > V[6][0]+R[1][1]*V[6][1]+R[1][2]*V[6][2],R[2][0]*V[6][0]+R[ > 2][1]*V[6][1]+R[2][2]*V[6][2]) > v8=(R[0][0]*V[7][0]+R[0][1]*V[7][1]+R[0][2]*V[7][2],R[1][0]* > V[7][0]+R[1][1]*V[7][1]+R[1][2]*V[7][2],R[2][0]*V[7][0]+R[ > 2][1]*V[7][1]+R[2][2]*V[7][2]) > body=utils.polyhedron((v1,v2,v3,v4,v5,v6,v7,v8),fixed= > False,color=(.6,.45,0),material="MLI",wire=False) > body.state.angVel=(0,0,0.1) > r=m.sqrt(aa**2+bb**2) > Rj=m.sqrt(r**2+cc**2) > Ri=0.05 > Rr=Rj*Ri/(Rj+Ri) > mu_rM=0.05 > mu_rG=1.08 > KN=E1*2*Ri*E1*2*Ri/(E1*2*Ri+E1*2*Ri) > KR=3*Ri**2*mu_rG**2*KN/4 > > O.engines=[ > ForceResetter(), > InsertionSortCollider([Bo1_Polyhedra_Aabb()]), > InteractionLoop( > [Ip2_FrictMat_FrictMat_MindlinPhys(en=.45,es=.45, > krot=KR,frictAngle=.7853)], > [Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True)] > ), > NewtonIntegrator(gravity=(0,0,-2.5e-4),damping=0), > ] > > The script is the same and also the version of yadedaily. > > I'll try with angMom. > > Thank you, > > Jonathan > > -- > 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 > -- 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 #335600]: Parallel simulation
Question #335600 on Yade changed: https://answers.launchpad.net/yade/+question/335600 Status: Expired => Answered Jan Stránský changed the question status: Sorry, probably everybody (based on my experience) missed the question -- 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 #335600]: Parallel simulation
Question #335600 on Yade changed: https://answers.launchpad.net/yade/+question/335600 Jan Stránský posted a new comment: Hi Mike, to run Yade in parallel, do everytihn the same as normally and just use -j option at startup. So to run Yade on 4 cores, use yade -j 4 or yade -j 4 script.py 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 #360288]: two-layer compaction
Question #360288 on Yade changed: https://answers.launchpad.net/yade/+question/360288 Jan Stránský proposed the following answer: Hello, for b in underlyingLayer: b.state.blockedDOFs = 'xyzXYZ' to release them again, replacing 'xyzXYZ' with empty string "" should work cheers Jan 2016-08-23 22:57 GMT+02:00 Weimin Song <question360...@answers.launchpad.net >: > New question #360288 on Yade: > https://answers.launchpad.net/yade/+question/360288 > > Dear all, > > I am trying to do the direct shear simulation of a two-layer > composite. First I compacted the underlying layer, and then I compacted the > upper layer, so that I evaluate the interlock effect on the interface. > Besides, in the process of the compaction of the upper layer, I want to > make the underlying layer particle static, how could I get that? > > -- > 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 > -- 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 #335600]: Parallel simulation
Question #335600 on Yade changed: https://answers.launchpad.net/yade/+question/335600 Status: Open => Answered Jan Stránský proposed the following answer: Hi Mike, currently YADE does not support MPI (distributed memory) parallelization. However, you can still use shared memory parallelization even on server (limited to one computer). cheers Jan 2016-08-29 8:32 GMT+02:00 ntthangvn <question335...@answers.launchpad.net>: > Question #335600 on Yade changed: > https://answers.launchpad.net/yade/+question/335600 > > ntthangvn posted a new comment: > Hi again, > > Following your guidance, I'm ok to run YADE on my computer (with 7 > cores) (Open MP). Is that ok for running YADE in server (I mean run YADE > in MPI)? > > Thank you for your help. > > Mike > > -- > 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 > -- 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 #402658]: add displacement according to their coordinate
Question #402658 on Yade changed: https://answers.launchpad.net/yade/+question/402658 Status: Open => Answered Jan Stránský proposed the following answer: You can prescribe velocity such that during next time step the particle gets to desired position b.state.vel = (desiredPos - b.state.pos) / O.dt cheers Jan 2016-10-01 14:53 GMT+02:00 张文卿 <question402...@answers.launchpad.net>: > New question #402658 on Yade: > https://answers.launchpad.net/yade/+question/402658 > > At the begining of each time step, I want to add displacement to every > particle according to their center coordinate. > just like this:for x direction，Displacement=a*x*timestep. > a is just a coefficient. > But I only find 'addMove'.Is there any way solve this problem.Thanks: ) > > -- > 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 > -- 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 #402592]: Plot Over Line with Paraview
Question #402592 on Yade changed: https://answers.launchpad.net/yade/+question/402592 Jan Stránský proposed the following answer: I *guess* you need a continuous data for this kind of analysis. I did a clip, than filter Deleunay3D and then Plot over line and got some line. This hack worked, although I have no idea if the results are correct or not :-) cheers Jan 2016-09-30 16:32 GMT+02:00 Yor1 <question402...@answers.launchpad.net>: > Question #402592 on Yade changed: > https://answers.launchpad.net/yade/+question/402592 > > Yor1 posted a new comment: > Hi Jan, > > I am sorry. This is the correct link : > https://www.dropbox.com/s/rw8y41chfrwb8fi/20k_sx60sy11_ > excavation-spheres.55302.vtu?dl=0 > > Best regards. > Jabrane. > > -- > 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 > -- 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 #403514]: YADE appears error in the Ubuntu 16.04
Question #403514 on Yade changed:
https://answers.launchpad.net/yade/+question/403514
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello,
this is a bug, probably coming from a different version of Eigen library in
Ubuntu 16.04. We will try to correct it.
Cheers
Jan
2016-10-26 14:07 GMT+02:00 liukeqi <question403...@answers.launchpad.net>:
> New question #403514 on Yade:
> https://answers.launchpad.net/yade/+question/403514
>
> Hello,
> In the past, I used Ubuntu 14.04 to install YADE(1.07.0) and built a
> model which used material "cpmmat". It worked fine. But now, I used Ubuntu
> 16.04 and install YADE(1.20.0) on it by the same way that using the command
> "sudo apt-get install yade", and use it to run the same *.py code(model).
> It does not work any more. If I use YADE on Ubuntu 16.04 to run the
> concrete example in trunk/examples/concrete, it still work fine. I think
> may be the YADE source code changed(from 1.07.0 to 1.20.0) and my model is
> incompatible with the new one. The error information and my code is
> following. Can you give me some suggestions?
>
>
> * the error information
> Welcome to Yade 1.20.0
> TCP python prompt on localhost:9000, auth cookie `cdsusy'
> XMLRPC info provider on http://localhost:21000
> Running script parameter.py
> python: /usr/include/eigen3/Eigen/src/SVD/JacobiSVD.h:648: void
> Eigen::JacobiSVD<MatrixType, QRPreconditioner>::allocate(Eigen::Index,
> Eigen::Index, unsigned int) [with _MatrixType = Eigen::Matrix<double, 3,
> 3>; int QRPreconditioner = 2; Eigen::Index = long int]: Assertion
> `(!(m_computeThinU || m_computeThinV) ||
> (MatrixType::ColsAtCompileTime==Dynamic))
> && "JacobiSVD: thin U and V are only available when your matrix has a
> dynamic number of columns."' failed.
> Aborted (core dumped)
>
> * my model code
> #!/usr/bin/python # This is server.py file
> import string
> from yade import plot,qt
> from yade.pack import *
> from yade import pack, plot
>
> O.materials.append(CpmMat(young=140e9,density=4800,
> frictionAngle=atan(0.8),poisson=.2,sigmaT=180e8,epsCrackOnset=0.033e-4,
> relDuctility=11000))
> initSize=0.02
>
> sp=pack.randomPeriPack(radius=.000625,initSize=Vector3(
> initSize,initSize,initSize))
> O.periodic=True
> sp.toSimulation()
>
> O.dt=1e-7
>
> plot.plots={'ey':('sy',)}
> def plotAddData():
> plot.addData(
> sy=p3d.stress[1],
> ey=p3d.strain[1],
> ex=p3d.strain[0]
> )
>
> EnlargeFactor=1.5
> EnlargeFactor=1.0
> O.engines=[
> ForceResetter(),
> InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=
> EnlargeFactor,label='bo1s')]),
> InteractionLoop(
> [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=
> EnlargeFactor,label='ig2ss')],
> [Ip2_CpmMat_CpmMat_CpmPhys()],
> [Law2_ScGeom_CpmPhys_Cpm()]),
> NewtonIntegrator(),
> Peri3dController(
>
> nSteps=4000,
> # how many time steps the simulation will last
> # after reaching
> nSteps do doneHook action
> doneHook='print
> "Simulation with Peri3dController finished."; O.pause()',
>
> label='p3d'
> ),
> PyRunner(command='plotAddData()',iterPeriod=1),
> ]
>
>
>
> p3d.stressMask=0b11
> p3d.goal=(0,-0.008,0,0,0,0)
> O.step()
> bo1s.aabbEnlargeFactor=ig2ss.interactionDetectionFactor=1.
> O.run(); O.wait()
> plot.plot(subPlots=False)
>
>
>
>
>
>
> --
> 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
>
--
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 #403470]: fatigue failure under cyclic loading using CPM
Question #403470 on Yade changed: https://answers.launchpad.net/yade/+question/403470 Jan Stránský proposed the following answer: In this case, just try some "normal" simulation, get the length of time step, how much real time does it take to compute it and extrapolate how much would 1 cycles take.. Jan 2016-10-25 21:37 GMT+02:00 Weimin Song <question403...@answers.launchpad.net >: > Question #403470 on Yade changed: > https://answers.launchpad.net/yade/+question/403470 > > Weimin Song posted a new comment: > Thank you, Jan. > > The cycles may be above 1 times. > > -- > 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 > -- 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 #403514]: YADE appears error in the Ubuntu 16.04
Question #403514 on Yade changed: https://answers.launchpad.net/yade/+question/403514 Jan Stránský posted a new comment: Hello, the problems is fixed in github version. You can use yadedaily (the change should be available in yadedaily in a few days), compile Yade form source by yourself, or wait until new yade package is added to Ubuntu (have no idea about the process). cheers Jan 2016-10-27 3:48 GMT+02:00 liukeqi <question403...@answers.launchpad.net>: > Question #403514 on Yade changed: > https://answers.launchpad.net/yade/+question/403514 > > Status: Answered => Solved > > liukeqi confirmed that the question is solved: > Thanks Jan Stránský, that solved my question. > > -- > 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 > -- 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 #403901]: Sample of pre-defined shaped particles
Question #403901 on Yade changed: https://answers.launchpad.net/yade/+question/403901 Status: Open => Answered Jan Stránský proposed the following answer: Hi Bruna, I am a bit confused now.. as I understand it, the script from my previous message does exactly what you described here: "a pack formed by particles that are formed by spheres".. Have you tried it? Of course, the script is minimal, so the shape of agglomerate is a cube (could be any shape), but each made of spheres (there are only 3 agglomerates, but it is easy to make arbitrary amount of them).. If I misunderstand something, could you please try it explain somehow differently? :-) maybe using a 2D sketch.. Thanks Jan 2016-11-09 14:03 GMT+01:00 bruna <question403...@answers.launchpad.net>: > Question #403901 on Yade changed: > https://answers.launchpad.net/yade/+question/403901 > > Status: Answered => Open > > bruna is still having a problem: > Hello, Jan. > > Thanks a lot for your answer. > > I have one question about this code. Is it possible that in instead of > the regularHexa I use a pack formed by particles that are formed by > spheres (the particle of my sample is going to be an agglomerate of > spheres). In other words, I want to create a sample that is formed by > particles made from spheres (the spheres that form one particle interact > by means of CohFrictMat interactions) and these spheres are responsable > to give a determined shape to my particles. I am trying to simulate > breakage of particles that has an specific shape. > > Once more thanks for all the help! > > Bruna. > > -- > 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 > -- 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
