Re: [Yade-users] [Question #692420]: Select bodies under post2d.extractor function
Question #692420 on Yade changed: https://answers.launchpad.net/yade/+question/692420 Status: Needs information => Solved Yuxuan Wang confirmed that the question is solved: Thank you so much Jan, I used extractVelocity=lambda b: flattener.normal(b,b.state.vel) if -0.15< b.state.pos[2] <= 0.15 else None and it worked perfectly. Best, Yuxuan -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #692497]: Point source for Thermal Engine?
New question #692497 on Yade: https://answers.launchpad.net/yade/+question/692497 Hello, I'm wondering is there a way to impose point or point-ish heat source for thermal engine, something similar to the point influx in the flow engine? Thanks. Best regards Jiannan -- 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 #692392]: i.geom.penetrationDepth = nan after load yade.gz file
Question #692392 on Yade changed: https://answers.launchpad.net/yade/+question/692392 Leonard posted a new comment: Great, thanks Bruno! Cheers Leonard -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #692392]: i.geom.penetrationDepth = nan after load yade.gz file
Question #692392 on Yade changed: https://answers.launchpad.net/yade/+question/692392 Chareyre posted a new comment: @Leonard The rationale for not saving/loading Un is that it is always unambiguously deduced from positions and sizes. We try to store only independent variables, not derived variables. What you report is not actually a problem, or it is only a problem of the way you output penetrationDepth. Insert a PyRunner to print Un, after NewtonIntegrator. The apparent issue will disappear. B On Thu, 20 Aug 2020 at 17:48, Bruno Chareyre wrote: > > (I did not really get your comment :-) > > I was pointing the fact that the notion of "value of X at iteration N" is > vague in the first place, since X can change at any point within one > time-loop. > Virtually you could insert multiple O.save() or any form of snapshot > within an iteration, interleaved with the usual engines, and you would get > different numbers in each of them depending on their relative positions in > the sequence of engines. > The above save/load issue is nearly the same thing I would say. > B > -- 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 #692392]: i.geom.penetrationDepth = nan after load yade.gz file
Question #692392 on Yade changed: https://answers.launchpad.net/yade/+question/692392 Chareyre posted a new comment: > (I did not really get your comment :-) I was pointing the fact that the notion of "value of X at iteration N" is vague in the first place, since X can change at any point within one time-loop. Virtually you could insert multiple O.save() or any form of snapshot within an iteration, interleaved with the usual engines, and you would get different numbers in each of them depending on their relative positions in the sequence of engines. The above save/load issue is nearly the same thing I would say. B -- 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 #692392]: i.geom.penetrationDepth = nan after load yade.gz file
Question #692392 on Yade changed: https://answers.launchpad.net/yade/+question/692392 Leonard 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
Re: [Yade-users] [Question #692392]: i.geom.penetrationDepth = nan after load yade.gz file
Question #692392 on Yade changed: https://answers.launchpad.net/yade/+question/692392 Status: Answered => Solved Leonard confirmed that the question is solved: Thanks Jan and Bruno for your comments, I tried a few more times in different ways (e.g., set O.dt=0 or not in both scripts) and still can't find the "exact state (i.e., the same penetrationDepth)'' between script1 and script2. What Jan explained make sense. It is good that I can use O.step() to re-gain the penetrationDepth after loading the sample. I will check further and hope this one additional step will not affect the final results compared to the case which doesn't involve this O.step. Thanks Leonard -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #692478]: Interrogation about generation of capillary files byMatlab program
New question #692478 on Yade: https://answers.launchpad.net/yade/+question/692478 Hi, in order to test a different contact angle for my simulation with capillary interaction I am using the Matalb programs writeCapFile.m and its two linked ones. I started first by using the inital case theta=0 and I am surprised by the non continuous evolution of the calculated maximum suction versus the rRatio: see the matlab output: For r = 1, maximum suction found: final uc = 9230 For r = 1.1, maximum suction found: final uc = 16749 For r = 1.25, maximum suction found: final uc = 19036 For r = 1.5, maximum suction found: final uc = 22843 - For r = 1.75, maximum suction found: final uc = 26389 For r = 2, maximum suction found: final uc = 21154 -- For r = 3, maximum suction found: final uc = 31738 For r = 4, maximum suction found: final uc = 41922 For r = 5, maximum suction found: final uc = 51901 For r = 10, maximum suction found: final uc = 88968 The jump down during the evolution between 1.75 and 2 is quite strange it is like two different slopes is it two different ways of calculating the data ? how this fact can influence the output files M(r=...) ? thanks in advance for your explanation Luc -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #692476]: Streamlined body
New question #692476 on Yade: https://answers.launchpad.net/yade/+question/692476 How to define geometry in case of flow around a cylinder? If you set the geometry the way I do it, then the time step becomes nan. ### def __init__(self): O.periodic = True O.cell.setBox(0.4,0.4,0.4) numspheres=1000 young = 5e6 density = 1000 mat1 = FrictMat(young=young,poisson=0.5,frictionAngle=radians(15),density=density,label='spheremat') O.materials.append(mat1) mat2 = FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='wallmat') O.materials.append(mat2) epsilon = 1e-08 minval = 0 + epsilon maxval = 0.4 - epsilon #wall coords, use facets for wall BC: v0 = Vector3(minval, minval, minval) v1 = Vector3(minval,minval,maxval) v2 = Vector3(maxval,minval,minval) v3 = Vector3(maxval,minval,maxval) v4 = Vector3(minval,maxval,minval) v5 = Vector3(minval,maxval,maxval) v6 = Vector3(maxval,maxval,minval) v7 = Vector3(maxval, maxval, maxval) lf0 = facet(vertices=[v0,v1,v2], material='wallmat') O.bodies.append(lf0) lf1 = facet(vertices=[v0,v2,v3], material='wallmat') O.bodies.append(lf1) uf0 = facet(vertices=[v4,v5,v6], material='wallmat') O.bodies.append(uf0) uf1 = facet(vertices=[v4,v6,v7], material='wallmat') O.bodies.append(uf1) ff0 = facet(vertices=[v1,v2,v6], material='wallmat') O.bodies.append(ff0) ff1 = facet(vertices=[v1,v6,v5], material='wallmat') O.bodies.append(ff1) bf0 = facet(vertices=[v0,v3,v7], material='wallmat') O.bodies.append(bf0) bf1 = facet(vertices=[v0,v7,v4], material='wallmat') O.bodies.append(bf1) ##cylinder## cyl = Cylinder(segment=[0,0,1],radius=0.05) b = Body() b.shape = cyl b.mat = O.materials['wallmat'] O.bodies.append(b) #spheres mn, mx= Vector3(minval + epsilon, minval + epsilon, minval + epsilon), Vector3(-0.06, maxval - epsilon, maxval - epsilon) sp = pack.SpherePack(); sp.makeCloud(mn,mx,rMean=0.00075,rRelFuzz=0.10, num=numspheres) O.bodies.append([sphere(center,rad,material='spheremat') for center,rad in sp]) sphereIDs = [b.id for b in O.bodies if type(b.shape)==Sphere] -- 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