Re: [Yade-users] [Question #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Ryota Nakanishi posted a new comment: Thank you Klaus, I understand why that is not suitable for calculating stresses in the tire. -- 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 #636678]: Jupyter notebook support
Question #636678 on Yade changed: https://answers.launchpad.net/yade/+question/636678 ar115 posted a new comment: Hi Jan, I understand from the standpoint of a developer who has much experience maybe jupyter is not important. But I think supporting jupyter would help to promote and advertise Yade more. And I thought that it should be trivial to include support for it as jupyeter is basically python and yade is nothing but python script. -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Status: Open => Answered Klaus Thoeni proposed the following answer: The current implementation of the pFacet is based on gridConnections, i.e., it is a discrete element with "internal contact forces/moments" in its edges. The formulation is not continuum based as the FEM. You could calculate the stresses from the "internal contact forces" but not sure how accurate this is going to be. In terms of computational effort, I think it would be the same. However, the pFacet should be more robust for contact tracking, especially for sliding contacts. HTH Klaus -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Status: Answered => Open Ryota Nakanishi is still having a problem: Hi Klaus, yeah, I know this paper and I'm thinking of using the method as an alternative plan. But why is it unsuitable for watching stresses in the tire by this method? I think representing the tire including the internal parts of it enable to see stresses in it. But of course, this approach needs a lot of DEM elements and then require a lot of computation time. This is one of the reasons why I want to use FEM-DEM coupling. Thanks a lot. Ryota -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Klaus Thoeni proposed the following answer: Hi Ryota, not sure if it helps but yade has deformable elements the tire could be represented with. If you are interested in the soil-tire interaction I would consider this an excellent approach. Nevertheless, if you are interested in stresses in the tire it would not be suitable. Have a look here [1] for examples of using such pFacet elements and here [2] for the theory. HTH Klaus [1] https://github.com/yade/trunk/tree/master/examples/pfacet [2] http://www.sciencedirect.com/science/article/pii/S0266114415001235 -- 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 #635871]: Bonded Particle Model using JCFpmMat
Question #635871 on Yade changed: https://answers.launchpad.net/yade/+question/635871 Status: Open => Answered Luc Scholtès proposed the following answer: Hi Xavier, I am not sure that sending you my script would really help since there is nothing more in it (a line before to launch a simulation I performed previously). I did not exactly get your point about the outFile but you should have a look at python doc about editing text file: https://docs.python.org/2/tutorial/inputoutput.html#reading-and-writing- files Regarding how to run a compresion test with yade, you should have a look at the uniaxial strainer: https://yade- dem.org/doc/yade.wrapper.html?highlight=uniaxial#yade.wrapper.UniaxialStrainer Here is a basic script for simulating such test with JCFPM: # -*- coding: utf-8 -*- # -*- encoding=utf-8 -*- from yade import pack, plot # SIMULATIONS DEFINED HERE packing (previously constructed) OUT='compressionTest_JCFPM' Simulation Control rate=-0.01 #deformation rate iterMax=1 # maximum number of iterations saveVTK=2000 # saving output files for paraview Material microproperties intR=1.1 # allows near neighbour interaction (can be adjusted for every packing) DENS=2500 # could be adapted to match material density: dens_DEM=dens_rock*(V_rock/V_particles)=dens_rock*1/(1-poro_DEM) -> packing porosity as to be computed? YOUNG=20e9 FRICT=7 ALPHA=0.1 TENS=1e6 COH=1e6 material definition def sphereMat(): return JCFpmMat(type=1,density=DENS,young=YOUNG,poisson=ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH) create the specimen pred=pack.inCylinder((0,0,0),(0,1,0),0.25) O.bodies.append(pack.regularHexa(pred,radius=0.025,gap=0.,material=sphereMat)) R=0 Rmax=0 nbSpheres=0. for o in O.bodies: if isinstance(o.shape,Sphere): nbSpheres+=1 R+=o.shape.radius if o.shape.radius>Rmax: Rmax=o.shape.radius Rmean=R/nbSpheres print 'nbSpheres=',nbSpheres,' | Rmean=',Rmean boundary condition (see utils.uniaxialTestFeatures bb=utils.uniaxialTestFeatures() negIds,posIds,longerAxis,crossSectionArea=bb['negIds'],bb['posIds'],bb['axis'],bb['area'] # ENGINES DEFINED HERE O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intR,label='Saabb')]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR,label='SSgeom')], [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1,label='interactionPhys')], [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(recordCracks=True,Key=OUT,label='interactionLaw')] ), UniaxialStrainer(strainRate=rate,axis=longerAxis,asymmetry=0,posIds=posIds,negIds=negIds,crossSectionArea=crossSectionArea,blockDisplacements=1,blockRotations=1,setSpeeds=0,stopStrain=0.1,dead=1,label='strainer'), GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.5, defaultDt=utils.PWaveTimeStep()), NewtonIntegrator(damping=0.4,label='newton'), PyRunner(iterPeriod=int(100),initRun=True,command='recorder()',label='data'), VTKRecorder(iterPeriod=int(saveVTK),initRun=True,fileName=OUT+'-',recorders=['spheres','jcfpm','cracks'],Key=OUT,label='vtk') ] # RECORDER DEFINED HERE def recorder(): yade.plot.addData({'i':O.iter, 'eps':strainer.strain, 'sigma':strainer.avgStress, 'tc':interactionLaw.nbTensCracks, 'sc':interactionLaw.nbShearCracks, 'te':interactionLaw.totalTensCracksE, 'se':interactionLaw.totalShearCracksE, 'unbF':utils.unbalancedForce()}) plot.saveDataTxt(OUT) # if you want to plot during simulation plot.plots={'i':('sigma')} plot.plot() # PREPROCESSING manage interaction detection factor during the first timestep and then set default interaction range ((cf. A DEM model for soft and hard rock, Scholtes & Donze, JMPS 2013)) O.step(); ### initializes the interaction detection factor SSgeom.interactionDetectionFactor=-1. Saabb.aabbEnlargeFactor=-1. coordination number verification numSSlinks=0 numCohesivelinks=0 for i in O.interactions: if isinstance(O.bodies[i.id1].shape,Sphere) and isinstance(O.bodies[i.id2].shape,Sphere): numSSlinks+=1 if i.phys.isCohesive : numCohesivelinks+=1 print "nblinks=", numSSlinks, " | nbCohesivelinks=", numCohesivelinks, " || Kcohesive=", 2.0*numCohesivelinks/nbSpheres # SIMULATION REALLY STARTS HERE strainer.dead=0 O.run(iterMax) Regards 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
Re: [Yade-users] [Question #635935]: creating a cloud of cuboids
Question #635935 on Yade changed: https://answers.launchpad.net/yade/+question/635935 rhaven 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 #635935]: creating a cloud of cuboids
Question #635935 on Yade changed: https://answers.launchpad.net/yade/+question/635935 Status: Answered => Solved rhaven confirmed that the question is solved: Hi Jan, thanks, that works perfect! -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Description changed to: Hi all, I would like to analyze contact problems between soil and tire by applying the coupling between DEM, i.e. YADE, and FEM like Ref[1]. I think what I need is surface coupling and volume coupling. I read the example scripts in Ref[1], but I couldn’t understand how to execute them. Judging from the response in Q site of YADE, I think new documentation is not yet available on the coupling of YADE with FEM, except for the recent example using escrypt [2]. Does anyone know the recent status of development on coupling of DEM and FEM by adding simple python script? Thanks, Ryota Ref: [1] J. Stránský & M.Jirásek. OPEN SOURCE FEM-DEM COUPLING. Proc. 18th International Conference ENGINEERING MECHANICS 2012, pp. 1237-1251. [2] https://yade-dem.org/doc/FEMxDEM.html -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Ryota Nakanishi posted a new comment: Hi Hongyang and Jan, Thank you for your help, Hongyan! I will try your scripts. And Thank you for your reply, Jan! I'll look forward to see your codes on github. Thanks a lot. Ryota -- 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 #635871]: Bonded Particle Model using JCFpmMat
Question #635871 on Yade changed: https://answers.launchpad.net/yade/+question/635871 Status: Answered => Open Xavier Thurman is still having a problem: Hi Luc, Thanks for your help, but I still have a question, what about the outFile ? Is it exported using export.text()? If possible, can you send a complete copy of your code to my email (xavierthur...@163.com)? Besides, in your paper, 'A DEM model for soft and hard rocks: role of grain interlocking on strength', you got a strain-stress curve of compression. I also tried to do this, but I only got displacement-force curve of uniaxial compression. I have no idea to calculate the strain-stress curve of uniaxial compression. Please give me some hints. Thanks a lot, Xavier -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Hongyang Cheng posted a new comment: Hi Ryota, I uploaded the scripts on github some time ago. Not sure it will work for others, but worth mentioning if one would like to have a try. Note that the script actually comes from Jan's contribution years ago... https://github.com/chyalexcheng/multiscale Cheers -- 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 #635935]: creating a cloud of cuboids
Question #635935 on Yade changed: https://answers.launchpad.net/yade/+question/635935 Status: Open => Answered Jan Stránský proposed the following answer: Hi Jesse, yes, periodicity might be the reason, but it is difficult to say without the code. Try for b O.bodies: b.state.pos = O.cell.wrapPt(b.state.pos) before saving Back to Bruno's answer, I see no problem there. If you first generate loose packing of spheres, then you can replace each sphere by any (smaller) shape and will get loose packing of whatever (elongated cuboids, ellipsoids...). 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 #636678]: Jupyter notebook support
Question #636678 on Yade changed: https://answers.launchpad.net/yade/+question/636678 Status: Open => Answered Jan Stránský proposed the following answer: Hello, I am sorry, but this message seems to me like an (actually not very good) advertisement.. Put here by "Ask a question" button, but there is no question, just a (IMO not true) statement. Do you want to create the support yourself? Do you want somebody else to create it? > Will make working with Yade much more user friendly Yade uses plain Python scripts for serious work, IPython is a sugar for testing and playing around, not the main user interface. So improving user friendliness is independent on IPython. IMO, documentation and support from the community is more valuable for the users than IPython or similar > it would be very very useful to support jupyter notebook as well "The Jupyter Notebook is an open-source web application that allows you to create and share documents that contain live code, equations, visualizations and explanatory text" I see no point how it would be "very very useful". But maybe I am just too old school.. anyway, Yade is open source, so if somebody wants/needs to use a different user interface, he/she is free to do so / ask somebody else to do so I apologize if I got something wrong 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 #632453]: Growing radius in randomDensePack
Question #632453 on Yade changed: https://answers.launchpad.net/yade/+question/632453 Jan Stránský posted a new comment: Hi Xavier, you can grow radius until they touch a surrounding particle A specific implementation could be something like this (not tested): # class ParticleGrower: def __init__(self, particles, multiplier=1.05): self.particles = set(p.id for p in particles) self.multiplier = multiplier def go(self): self.check() self.grow() def check(self): for pid in self.particles: p = O.bodies[pid] intrs = p.intrs() if intrs: # nonzero overlap growParticle(p.id, 1./self.multiplier) # decrease particle back to non-overlapping configuration self.particles.discard(pid) # do not deal with this particle any more def grow(self): for pid in self.particles: growParticle(pid, self.multiplier) aloneParticles = [...] # list of particles, that are alone and make problems grower = ParticleGrower(particles) O.engines = [ ... PyRunner(iterPeriod=1,command='grower.go()') ] while grower.particles: O.run(10,True): print len(grower.particles) # 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
[Yade-users] [Question #636712]: import a clump and use it with makeClumpCloud
New question #636712 on Yade: https://answers.launchpad.net/yade/+question/636712 Hi Sorry if this has already been answered, I didnt find the answer anywhere I am importing a clump using c0 = yimport.textClumps('/tmp/clump0.txt') and I want to use it with makeClumpCloud however makeClumpCloud need a SpherePack. I can search through O.bodies to find the clump body, however how do I convert this into a SpherePack object? thanks in advance Jesse -- 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Jan Stránský proposed the following answer: Hi Ryota, I will put my codes to github soon 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 #636692]: soil-tire contact by FEM-DEM coupling
Question #636692 on Yade changed: https://answers.launchpad.net/yade/+question/636692 Status: Open => Answered Hongyang Cheng proposed the following answer: Dear Ryota, I had worked on soil-geosynthetic interaction based on Ning's approach. Please see [1] for the surface coupling scheme. For further details, please contact me via h.ch...@utwente.nl. [1] https://link.springer.com/chapter/10.1007/978-981-10-1926-5_47 Best regards, Hongyang -- 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 #632475]: How to measure the void ratio for a select range?
Question #632475 on Yade changed: https://answers.launchpad.net/yade/+question/632475 Jan Stránský posted a new comment: Hello, please be more specific where the problem is. How to do it theoretically? How to code it? Also a minimum working example would be nice, describing where you don't know how to continue and what you want to achieve. Thanks for more info 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 #632475]: How to measure the void ratio for a select range?
Question #632475 on Yade changed: https://answers.launchpad.net/yade/+question/632475 Status: Open => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Open' state without activity for the last 15 days. -- 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 #632453]: Growing radius in randomDensePack
Question #632453 on Yade changed: https://answers.launchpad.net/yade/+question/632453 Status: Open => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Open' state without activity for the last 15 days. -- 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 #636692]: soil-tire contact by FEM-DEM coupling
New question #636692 on Yade: https://answers.launchpad.net/yade/+question/636692 Hi all, I would like to analyze contact problems between soil and tire by applying the coupling between DEM, i.e. YADE, and FEM like Ref[1]. I think what I need is so-called surface coupling and volume coupling. I read the example scripts in Ref[1], but I couldn’t understand how to execute them. Judging from the response in Q site of YADE, I think new documentation is not yet available on the coupling of YADE with FEM, except for the recent example using escrypt [2]. Does anyone know the recent status of development on coupling of DEM and FEM by adding simple python script? Thanks, Ryota Ref: [1] J. Stránský & M.Jirásek. OPEN SOURCE FEM-DEM COUPLING. Proc. 18th International Conference ENGINEERING MECHANICS 2012, pp. 1237-1251. [2] https://yade-dem.org/doc/FEMxDEM.html -- 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 #635935]: creating a cloud of cuboids
Question #635935 on Yade changed: https://answers.launchpad.net/yade/+question/635935 rhaven gave more information on the question: Hello, Ive made small progress with method 2 and creating the clumps which I will later use with makeClumpCloud, I generate a small aggregate using makeCloud. The result looks like http://pasteall.org/pic/index.php?id=115961 I then generate a clump from these spheres and export. However when I import the geometry again it looks like this http://pasteall.org/pic/index.php?id=115962 Is this occurring due to periodicity? How might I generate a small aggregate? thank you in advance Jesse -- 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