Re: [Yade-users] [Question #293342]: Problem with simulation
Question #293342 on Yade changed: https://answers.launchpad.net/yade/+question/293342 Status: Answered => Open Alex is still having a problem: Hello Jérôme and Klauss, Thank you for the answers. They have helped a lot. I'm sorry for the delay in replying and for the title... The next time I will be more specific and more original. ;-) I rewrote my code and I continue with the same problem, that is, when I run the simulation the cylinder goes inside the spheres and that's not what I want. Once again, I ask if someone could help in this problem. Below is my new code. import sys from yade import ymport from yade import qt from numpy import linspace material definition def sphereMat(): return CohFrictMat(young=8e5,poisson=0.3,density=4e3, frictionAngle=radians(30),normalCohesion=4e4,shearCohesion=1e5,momentRotationLaw=False,label='sphereMat') def gridNodeMat(): return CohFrictMat(young=8e5,poisson=0.3,density=4e3,frictionAngle=radians(30), normalCohesion=1e5,shearCohesion=1e5,momentRotationLaw=True,label='gridNodeMat') def gridCoMat(): return CohFrictMat(young=8e5,poisson=0.3,density=4e3, frictionAngle=radians(30), normalCohesion=1e5,shearCohesion=1e5,momentRotationLaw=True,label='gridCoMat') end material definition ### O.engines=[ ForceResetter(), InsertionSortCollider([ Bo1_GridConnection_Aabb(), Bo1_Sphere_Aabb() ]), InteractionLoop( Geometric interactions [ Ig2_GridNode_GridNode_GridNodeGeom6D(), Ig2_Sphere_GridConnection_ScGridCoGeom() ], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True), ], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment(label='law'), Law2_ScGridCoGeom_CohFrictPhys_CundallStrack() ] ), Motion equation NewtonIntegrator(damping=0.4,gravity=[0,0,0]), PyRunner(command='main()',iterPeriod=10) ] O.dt=1e-7 rCyl=0.025 nL= 2 L=0.3 Create all nodes first : nodesIds=[] for i in linspace(0,L,nL): nodesIds.append( O.bodies.append( gridNode([0.3,i + 0.4,2.02],rCyl,wire=False,fixed=False,material=gridNodeMat) ) ) Now create connection between the nodes for i,j in zip( nodesIds[:-1], nodesIds[1:]): O.bodies.append( gridConnection(i,j,rCyl, material=gridCoMat ) ) read a sphere's package PACKING = '10K_P1' O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0),material=sphereMat)) qt.View(); def main(): global Fn,Ft O.bodies[0].dynamic=False O.bodies[0].state.vel[2]=-500.0 O.bodies[1].state.vel[2]=-500.0 O.bodies[2].state.vel[2]=-500.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
[Yade-users] [Question #293342]: Problem with simulation
New question #293342 on Yade: https://answers.launchpad.net/yade/+question/293342 Hello friends, I created the simulation below and I had a problem. My objective is to make the cylinder collide with spheres. But, I can't do this. When I run the simulation the cylinder goes inside the spheres. I ask if someone could help in this problem. import sys from yade import ymport, utils , plot import math Material microproperties -> SHALE! intR=1.22955 DENS=2659.6453 YOUNG=40e9 FRICT=25 ALPHA=0.25 TENS=10e5 COH=10e6 aRes=1e-3 damp=0.2 jointFRICT=25 jointDIL=15 material definition def sphereMat(): return JCFpmMat(type=1,density=DENS,young=YOUNG,poisson = ALPHA,frictionAngle=radians(FRICT), tensileStrength=TENS,cohesion=COH,jointNormalStiffness=YOUNG/10.,jointShearStiffness=0.2*YOUNG/10., jointTensileStrength=0.,jointCohesion=0.,jointFrictionAngle=radians(jointFRICT),jointDilationAngle=radians(jointDIL)) O.engines=[ ForceResetter(), InsertionSortCollider([ Bo1_ChainedCylinder_Aabb(), Bo1_Sphere_Aabb() ]), InteractionLoop( # Geomeetric interactions [Ig2_ChainedCylinder_ChainedCylinder_ScGeom6D(), Ig2_Sphere_ChainedCylinder_CylScGeom6D(), Ig2_Sphere_Sphere_ScGeom6D(), Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom() ], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True), Ip2_FrictMat_FrictMat_FrictPhys() ], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment(label='law'), Law2_ScGeom_FrictPhys_CundallStrack(), Law2_L3Geom_FrictPhys_ElPerfPl(), ] ), ## Motion equation NewtonIntegrator(damping=0.4,gravity=[0,0,0]), PyRunner(command='main()',iterPeriod=10) ] O.dt=1e-7 # create a cylinder cil = chainedCylinder(begin=Vector3(0.3,0.4,2.02), end=Vector3(0.3,0.6,2.02),radius=0.025,color=Vector3(0.6,0.5,0.5)) # read a sphere's package PACKING = '10K_P1' O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0),material=sphereMat)) yade.qt.View(); def main(): global Fn,Ft O.bodies[0].dynamic=False O.bodies[0].state.vel[2]=-500.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
Re: [Yade-users] [Question #275973]: Mannually Update Location of Components of Clumps
Question #275973 on Yade changed: https://answers.launchpad.net/yade/+question/275973 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: HI, It looks you are doing FEM/DEM modeling pioneered by Munjiza‘s (FEs within DEM). I am guessing the so-called PFacet element in Yade [1] would work in this case, but It seems you have to wait until it is ready in trunk. [1] https://www.researchgate.net/publication/281203746_Geotextiles_and_Geomembranes Cheers, Alex On Tue, Dec 1, 2015 at 4:07 AM, MikeZhao < question275...@answers.launchpad.net> wrote: > Question #275973 on Yade changed: > https://answers.launchpad.net/yade/+question/275973 > > MikeZhao gave more information on the question: > Actually I just tried to use the simplest example to describe the > problem I've encountered: > > I am working on FEM-DEM coupling now, and I want to use DEM to build a > contact detector for different FE meshes. So I first transferred > elements of FEM into utils.tetraPoly in DEM. To make calculation more > efficient, I made tetraPoly that belongs to the same FE mesh a clump. > And I need to update the positions of those tetraPoly in DEM using the > positions of FEM's elements, so I need to change the positions of those > tetraPoly mannually.. > > -- > 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 #264174]: how is shear force related to poisson ratio?
Question #264174 on Yade changed: https://answers.launchpad.net/yade/+question/264174 Hongyang Alex Cheng proposed the following answer: Hi, Shenyang Like Bruno said. you should take a look at the functors you use. The contact springs are defined by the law2 and the Ip2 functors. The meaning of Young's modulus and poisson ratio varies from different set of law2 and lp2 functors in use (see linear contact and mindlin contact for instance). If you use the linear contact law changing poisson ratio will only affect the shear stiffness. Note that adjusting material constants doesn't have effect on existing contacts. Alex On Fri, Mar 27, 2015 at 10:41 AM, Shenyang Cai < question264...@answers.launchpad.net> wrote: > Question #264174 on Yade changed: > https://answers.launchpad.net/yade/+question/264174 > > Shenyang Cai posted a new comment: > Hi Bruno, > > I'm still confused. When two spheres are in contact, I suppose there > must be a normal stifness and a shear stiffness(though it may equal > zero) between them, regardless of the IPhys and Law functors. > > My idea is to change some parameters in the material class to enlarge > the shear stiffness between the spheres. And I thought "poisson" is a > factor. If changing "poisson" does not lead to the change of shear > stiffness, then what paremeters should I focus on to achieve my goal? > > > Thanks! > Shenyang > > -- > You received this question notification because you are a member of > yade-users, which 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 you are a member of yade-users, which 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 #263780]: Changing the material while running
Question #263780 on Yade changed: https://answers.launchpad.net/yade/+question/263780 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: Hi, I don't know how to do it while running the simulation if you are referring to change material types. But you could export the geometry of your "upper cube" then start over with newly defined material types. For how to change material constants, you might want to take a look at https://answers.launchpad.net/yade/+question/238797. BTW to fix the error message, try def change(): for b in O.bodies: if b.state.pos[2] <= 1: b.mat.label = 'something' # this will not change the material type alex -- You received this question notification because you are a member of yade-users, which 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 #261685]: Function ExportSpheres
Question #261685 on Yade changed: https://answers.launchpad.net/yade/+question/261685 Status: Answered => Solved Alex confirmed that the question is solved: Hello, Thank you for your example, it allowed me to find a solution. I changed method using: export.textExt () export.textExt('../../tmp/data.txt',format='x_y_z_r_matId',attrs=('b.mat.label')) resultat.txt #format x_y_z_r_matId 0.040452-0.624095 -0.583232 0.1 1 0.333673-0.849889 -0.113725 0.1 1 -0.382878 0.13828 -0.819718 0.1 1 -0.0962031 -0.787703 -0.382730.1 1 0.26245 -0.638839 -0.619982 0.1 1 0.173884-0.710586 -0.429939 0.1 1 0.0083 -294.141137.46 0.1 2 -342.899339.512 690.185 0.1 2 -0.10.0876388 -31.35280.1 2 0.0083 0.181458-31.17590.1 2 342.915 339.512 690.185 0.1 2 0.1166670.0876388 -31.35280.1 2 -0.256126 -0.666963 -0.347073 0.1 3 0.36668 -0.814118 0.1215170.1 3 -108.708-51.2151-728.1960.1 3 0.0887508 -0.416072 -0.757411 0.1 3 0.399858-0.794111 -0.386862 0.1 3 Thank you for this help Alex, -- You received this question notification because you are a member of yade-users, which 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 #261685]: Function ExportSpheres
Question #261685 on Yade changed: https://answers.launchpad.net/yade/+question/261685 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: Hi, For that purpose, I always assign duplicates of material type to groups of spheres, sphMat1 = O.materials.append(FrictMat(young=Esph,poisson=vsph,frictionAngle=radians(Phisph),density=Dsph)) sphMat2 = O.materials.append(FrictMat(young=Esph,poisson=vsph,frictionAngle=radians(Phisph),density=Dsph)) ... then export b.mat.id with VtkExporter. b.mask may also help [1]. In fact with VTKExporter you can define any attribute you like and visualize it in Paraview. [1]https://yade- dem.org/doc/yade.wrapper.html#yade.wrapper.Body.groupMask On Thu, Feb 5, 2015 at 7:07 PM, Alex wrote: > Question #261685 on Yade changed: > https://answers.launchpad.net/yade/+question/261685 > > Status: Answered => Open > > Alex is still having a problem: > Hello, > > Thank you for the quick response. > > With your response I get color in my data export, thank you... by cons > the position of the spheres and color are separate and my goal is to > identify the cluster of spheres by their colors. > > In my simulation, I test the homogeneity of a mixture of croquette (5 > samples) > the color is a way to identify a cluster of spheres (symbolizing the > croquettes). > > is there another way to identify a group of spheres? > > -- > You received this question notification because you are a member of > yade-users, which 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 you are a member of yade-users, which 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 #261685]: Function ExportSpheres
Question #261685 on Yade changed: https://answers.launchpad.net/yade/+question/261685 Status: Answered => Open Alex is still having a problem: Hello, Thank you for the quick response. With your response I get color in my data export, thank you... by cons the position of the spheres and color are separate and my goal is to identify the cluster of spheres by their colors. In my simulation, I test the homogeneity of a mixture of croquette (5 samples) the color is a way to identify a cluster of spheres (symbolizing the croquettes). is there another way to identify a group of spheres? -- You received this question notification because you are a member of yade-users, which 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 #261685]: Function ExportSpheres
Question #261685 on Yade changed: https://answers.launchpad.net/yade/+question/261685 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: Hi Alex You may try VtkExporter.exportSpheres (what = [('pos', 'b.state.pos'), ('color', 'b.shape.color')]) to export positions and colors. BTW, why do you want colors? Alex -- You received this question notification because you are a member of yade-users, which 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 #261685]: Function ExportSpheres
New question #261685 on Yade: https://answers.launchpad.net/yade/+question/261685 Hello, via a function in a simulation script I would like to export the data of spheres with their colors (or login if feasible) I use exportSpheres () function that works perfectly, for cons I can not allow specifying the color argument? here is my function: def ExportData (): print "Export data ..." vtkExporter export.VTKExporter = ('../../ tmp / vtkExporterTesting') # VtkExporter.exportSpheres (what = [('dist', 'b.state.pos.norm ()')]) # <--- this line works vtkExporter.exportSpheres (what = [('all')]) # <--- this line return error = (ValueError: too many values to unpack) My goal is to have a function that allows me to export the position of the spheres with their colors or identifier (IDS). thank you for your return. Alexandre, -- You received this question notification because you are a member of yade-users, which 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 #261082]: Save/compute average strain 0.5*(∇u+∇u^t)
Question #261082 on Yade changed: https://answers.launchpad.net/yade/+question/261082 Status: Answered => Solved Hongyang Alex Cheng confirmed that the question is solved: Thanks, Bruno. I'll try TW.getVolPoroDef(). I am simulating triaxial test on a soilbag in the periodic boundary condition. The macro strain is given as Vector3. I am just curious whether it has non zero off-diagonal parts because the boundary is not smooth. -- You received this question notification because you are a member of yade-users, which 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 #261082]: Save/compute average strain 0.5*(∇u+∇u^t)
Question #261082 on Yade changed: https://answers.launchpad.net/yade/+question/261082 Hongyang Alex Cheng confirmed that the question is solved: Thanks Bruno Chareyre, that solved my question. -- You received this question notification because you are a member of yade-users, which 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 #261082]: Save/compute average strain 0.5*(∇u+∇u^t)
Question #261082 on Yade changed: https://answers.launchpad.net/yade/+question/261082 Status: Answered => Open Hongyang Alex Cheng is still having a problem: Hi, Bruno >This is exactly what we are doing (with regular triangulation), else the concept of "closest neighbours" would be undefined. Thanks. Now I know how the TW, MMA and KLA work. >Since you can read it (as in the line above), you can also write it into a file. Is that what you mean by "save"? I would like to compute the macro deviatoric strain. The function KLA.computeParticlesDeformation() [1] returns 'ParticleDeformation' but only ouputs 'Total volume ', 'grad_u', and 'sym_grad_u' on the screen. If it is saved in memory can you tell me how to get it? [1] https://github.com/yade/trunk/blob/283afb47ed4fdb86047618f0cb3ac892c040ef88/lib/triangulation/KinematicLocalisationAnalyser.cpp#L837 Alex -- You received this question notification because you are a member of yade-users, which 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 #261082]: Save/compute average strain 0.5*(∇u+∇u^t)
New question #261082 on Yade: https://answers.launchpad.net/yade/+question/261082 Dear all, The TesselationWrapper is able to compute deformation from given chronological states/positions of a packing. Though it writes particle-based strain tensors to in vtk format and returns the average strain (sym_grad_u), I didn't find a way to save the strain in Matrix3 format so that I can compute the deviatoric part. An alternative way is to triangulate the packing (Regular/Delaunay) following Bagi's method such that the products would be similar to what TesselationWrapper produces. The Delaunay triangulation can be easily imported from scipy [1], but the regular triangulation [2] looks more preferable because it is also used in Yade. Unfortunately I don't know how to import it to Yade. Is it possible to save the average strain returned from the TesselationWrapper? If not, could anyone tell me how to import the regular triangulation package to Yade without touching the source code? [1] http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.spatial.Delaunay.html [2] http://doc.cgal.org/latest/Triangulation_3/ Cheers Alex -- You received this question notification because you are a member of yade-users, which 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 #253692]: Assign frictAngle for Ip functor MindlinPhys (MatchMaker instance)
Question #253692 on Yade changed: https://answers.launchpad.net/yade/+question/253692 Status: Answered => Solved Hongyang Alex Cheng confirmed that the question is solved: Thanks Klaus Thoeni, that solved my question. -- You received this question notification because you are a member of yade-users, which 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 #253692]: Assign frictAngle for Ip functor MindlinPhys (MatchMaker instance)
New question #253692 on Yade: https://answers.launchpad.net/yade/+question/253692 Hi all, By default, Ip2_FrictMat_FrictMat_MindlinPhys takes the minimal of the Phis of two bodies in contact as the frictional angle. I found that MatchMaker [1] might be able to change that default because frictAngle is an instance of MatchMaker according to [2]. I tried to simply assign 'avg' to the frictAngle of this Ip functor, but it gives the following error message. Not convertible to MatchMaker TypeError: No registered converter was able to produce a C++ rvalue of type boost::shared_ptr from this Python object of type str Could anyone teach me how to use MatchMaker class and how to assign values for the frictionAngle of Ip2_FrictMat_FrictMat_MindlinPhys? Is there any example that I can refer to? Thank you, Alex [1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.MatchMaker [2] https://yade-dem.org/doc/yade.wrapper.html?highlight=mindlinphys#yade.wrapper.Ip2_FrictMat_FrictMat_MindlinPhys.frictAngle -- You received this question notification because you are a member of yade-users, which 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 #252337]: Segmentation fault on reload
Question #252337 on Yade changed: https://answers.launchpad.net/yade/+question/252337 Hongyang Alex Cheng proposed the following answer: Hi Anna, I tried on my computer (also ubuntu 14.04 LTS). Same problem appeared first, but after I deleted qt.View() in line 89 reload worked normally. By the way, I am using git-2d27f33 version. Cheers Alex -- You received this question notification because you are a member of yade-users, which 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 #252191]: Saving facets to individual .vtu files
Question #252191 on Yade changed: https://answers.launchpad.net/yade/+question/252191 Hongyang Alex Cheng posted a new comment: Using groupMask is easier [1], I vote for Klaus's method. https://yade- dem.org/doc/yade.wrapper.html?highlight=vtkrecorder#yade.wrapper.VTKRecorder.mask Alex -- You received this question notification because you are a member of yade-users, which 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 #252191]: Saving facets to individual .vtu files
Question #252191 on Yade changed: https://answers.launchpad.net/yade/+question/252191 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: Hi, Colin You might want to take a look at the VTKExporter instead [1]. I haven't used facet before, but I did manage to export into Paraview some spheres with specified Ids. You can specify your target bodies in 'ids' in and the data you wanna extract in 'what'. 'ids' and 'what' are params of [1]. [1] https://yade- dem.org/doc/yade.export.html?highlight=vtkexporter#yade.export.VTKExporter.exportFacets Regards, Alex -- You received this question notification because you are a member of yade-users, which 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 #251712]: Two problems for post-processing
Question #251712 on Yade changed: https://answers.launchpad.net/yade/+question/251712 Hongyang Alex Cheng posted a new comment: Thanks, Bruno. This is really good to know and the eigen value&vector fuction as well. :-) chyalexch...@gmail.com From: Bruno Chareyre Date: 2014-07-16 22:36 To: chyalexcheng Subject: Re: [Question #251712]: Two problems for post-processing Your question #251712 on Yade changed: https://answers.launchpad.net/yade/+question/251712 Bruno Chareyre posted a new comment: >I found the products of bodyStressTensors(), per-particle stress matrices, not symmetric. Is this supposed to be normal? It tends to symmetry only as particles tend to static equilibrium. BC -- You received this question notification because you asked the question. You received this question notification because you are a member of yade-users, which 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 #251712]: Two problems for post-processing
Question #251712 on Yade changed: https://answers.launchpad.net/yade/+question/251712 Hongyang Alex Cheng posted a new comment: @ Bruno Chareyre Thanks for explaining more on getting p and q. I use 0.5*Sij*Sji to calculate deviatoric stress invariant J2 and then q. [1] I mentioned about Eigen because I was expecting some nonlinear equations solver in the library Eigen3. Sorry for making it confusing. BTW, I found the products of bodyStressTensors(), per-particle stress matrices, not symmetric. Is this supposed to be normal? [1] http://www.rockmechs.com/stress-strain/stress/deviatoric-stress-and- invariants/ -- You received this question notification because you are a member of yade-users, which 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 #251712]: Two problems for post-processing
Question #251712 on Yade changed: https://answers.launchpad.net/yade/+question/251712 Status: Answered => Solved Hongyang Alex Cheng confirmed that the question is solved: Hi Jerome Yes, I want mean and deviatoric stresses to render particles (although tessellation seems the best choice, but I haven't fully understand how it works). With principal stresses, mean and deviatoric stresses can always be easily derived. Glad to know we are targeting the same thing. Please commit. Thank you. Alex -- You received this question notification because you are a member of yade-users, which 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 #251712]: Two problems for post-processing
Question #251712 on Yade changed: https://answers.launchpad.net/yade/+question/251712 Hongyang Alex Cheng confirmed that the question is solved: Thanks Jérôme Duriez, that solved my question. -- You received this question notification because you are a member of yade-users, which 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 #251712]: Two problems for post-processing
Question #251712 on Yade changed: https://answers.launchpad.net/yade/+question/251712 Hongyang Alex Cheng posted a new comment: Thanks, Anton Due to my limited knowledge I haven't seen any formulation for contact stress like in [1], but I'll check it. If I am not mistaken, [1] returns normal stress and tangential stress at particle contact level. I've been thinking whether it's better to export mean stress p and deviatoric stress q for post-processing instead of normStress and shearStress. That's why I raised my second question. Here is the nonlinear equations I wanna solve where si (i=1,2,3) are principal stresses, I1, I2 and I3 are three stress invariants. s1+s2+s3 = I1 s1*s2+s2*s3+s3*s1=I2 s1*s2*s3=I3 [1] https://github.com/yade/trunk/blob/master/pkg/dem/Shop_02.cpp#L176 -- You received this question notification because you are a member of yade-users, which 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 #251712]: Two problems for post-processing
New question #251712 on Yade: https://answers.launchpad.net/yade/+question/251712 Hello, 1. With VTKRecorder it is possible to export normal stress and shear stress of a particle but what are the physical meanings of normStress and shearStress and how are these stresses computed? I checked the source code and found the following sentences might give some hints, but I failed to fully understand them. In VTKRecorder.cpp: Line 454 and 455 const Vector3r& stress = bodyStates[b->getId()].normStress; const Vector3r& shear = bodyStates[b->getId()].shearStress; 2. The function bodyStressTensors() calculates per-particle micro-stress which is a 3 by 3 matrix. I think it would be better to calculate mean principal stress p and deviatoric stress q from it, which will need a nonlinear equations solver. I tried to import scipy.optimize [1] but it led to core dumped. Does eigen library has some modules to solve nonlinear equations? Are there any software that is importable to Yade will work for this? I am using git version 2014-06-29.git-de4c01a. [1] http://docs.scipy.org/doc/scipy/reference/optimize.html#module-scipy.optimize Thank you, Alex -- You received this question notification because you are a member of yade-users, which 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 #251125]: Question about units
Question #251125 on Yade changed: https://answers.launchpad.net/yade/+question/251125 Hongyang Alex Cheng proposed the following answer: @Bruno Chareyre: Oops. sorry for getting it wrong. In that sense, we can use quite a variety of unit systems. Thank you for correcting me. On Fri, Jul 4, 2014 at 12:46 AM, Bruno Chareyre < question251...@answers.launchpad.net> wrote: > Question #251125 on Yade changed: > https://answers.launchpad.net/yade/+question/251125 > > Bruno Chareyre proposed the following answer: > @Alex-Cheng: no, no, no. Yade is not using any unit. > > @Emon, google is your friend if you don't like to do the math [1]: > force: F dyne | dyn | g·cm/s2 | 10−5 N > > [1] > > http://en.wikipedia.org/wiki/Centimetre%E2%80%93gram%E2%80%93second_system_of_units > > -- > You received this question notification because you are a member of > yade-users, which 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 you are a member of yade-users, which 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 #251125]: Question about units
Question #251125 on Yade changed: https://answers.launchpad.net/yade/+question/251125 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: Hi, Yade is using meter-kilogram-second-newton system. Check the attributes in O.Materials and O. interactions might help. regards, alex -- You received this question notification because you are a member of yade-users, which 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 #250922]: Bug report in function VTKExporter.exportInteractions()
Question #250922 on Yade changed: https://answers.launchpad.net/yade/+question/250922 Status: Answered => Solved Hongyang Alex Cheng confirmed that the question is solved: Hi Jan, Thanks. I am installing yade from the newest source directory now. export.degrees helped. Guess I gotta learn more about python. Alex -- You received this question notification because you are a member of yade-users, which 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 #250922]: Bug report in function VTKExporter.exportInteractions()
New question #250922 on Yade: https://answers.launchpad.net/yade/+question/250922 Dear all, I've been playing with VTKExporter lately and found that there might be a sentence "i = O.interaction[ii,jj]" missing in the func exportInteractions [1] for the looping which outputs Vector3 data. [1] https://yade-dem.org/doc/_modules/yade/export.html#VTKExporter.exportInteractions It appears that all the Vector3 data (i.phys.normalForce) exported with this function are the same (same as in the first interaction). The following is part of the script where I think the problem might be: --- # Line 34 for i in O.interactions: if i.isReal: break --- This loop returns i as the first interaction. It probably is the reason that I got same value for all i.phys.normalForce. --- # Line 45 elif isinstance(test,Vector3): outFile.write("\nVECTORS %s double\n"%(name)) for ii,jj in intrs: (should we add "i = O.interaction[ii,jj])" here?) v = eval(command) outFile.write("%g %g %g\n"%(v[0],v[1],v[2])) In addition to that, does someone know where to find this python script in the source code files? If I am right about this bug, I think I have to modify the code file and install Yade again. I have a question about VTKExporter. It is irrelevant to above. When we call VTKExporter class functions, exportInteractions(params) for instance, is it possible to derive some data from O.interactions's attributes like defining the param what as "what=[('contactPhi','degrees(acos(i.geom.normal[2]))')]"? Every time I do that periodically in PyRunner. It reports "name "degrees" ( or acos) is not defined bla bla bla (math module didn't work?). Could someone help me out? :-) Regards, alex -- You received this question notification because you are a member of yade-users, which 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 #250783]: some questions about 'O.dt' and 'PyRunner'
Question #250783 on Yade changed: https://answers.launchpad.net/yade/+question/250783 Status: Open => Answered Hongyang Alex Cheng proposed the following answer: Hi Zuoqiang, Thanks for explaining. If I am not mistaken, what you want is a series of simulation states at which the y_strain is specified beforehand so that you can compare different triax simulation results (e.g. different psd) at the same strain level. I've never thought about before, but if a y_strain-specified simulation state is what you are after it's always possible to use PyRunner Engine to check y_strain periodically and save your variables. regards, alex chyalexch...@gmail.com From: Fu zuoguang Date: 2014-06-28 15:46 To: yade-users Subject: Re: [Yade-users] [Question #250783]: some questions about 'O.dt' and 'PyRunner' Question #250783 on Yade changed: https://answers.launchpad.net/yade/+question/250783 Status: Answered => Open Fu zuoguang is still having a problem: Dear Hongyang Alex Cheng: I am sorry for answering your question so late.Yes, you are right,‘GlobalStiffnessTimeStepper’ has done very well in ensuring the stability of simulations. My reason of choosing ‘dynDt’ has nothing to do with controlling the stability but has close relationship with postprocess.Take the tri-axial simulation in Yade as a simple example, after the isotropic compression, the process of deviatoric shearing can be carried out just by these two conditions: (1) Laternal stress which is equal to isotropic pressure.(always choose x and z axis) (2) Axial strain which is only determined by a certain strain rate.(always choose y axis) And the final state of this simulation is always checked by a special axial strain. It is obvious that the relationship between the unit simulation step delta-t and the final axial strain can be written as follows: (Delta-t)*(total-steps)*(y-strainrate)=final axial strain (1) The y-strainrate and the final axial strain are always treated as the prediction of the simulation and are known. For this reason (Delta-t )*(total-steps) should be a constant to establish the formula(1). In yade, a changed ‘timestepper’ can hardly provide us a certain constant (total-steps) for controlling the total steps(cost) in simulation with a given smallest particle size and a total number. So I hope that there exists a constant delta-t to offer much help to not only the stable of simulation but also the fully controlled number of steps. In this condition, as the absolutely necessary portion of total result, each step or each interval of steps in data can represent a specific state into the whole process and then be as one of the specific representations of all the simulation objectives. In the other hand, the attributes of particles may not the same in different simulations. If they can use the same delta-t as the unit step, each element of the two result data output by ‘iterPeriod’ may have the same meaning. It can bring convenience to comparison of the data derived from different simulations. I think the chosen delta-t, which is smaller than ‘timestepper’ can meet what I want, even though it may bring a little more cost of simulation. I am not very right and there may be a better way for my purpose. Please tell me a little more. Fu -- You received this question notification because you are a member of yade-users, which 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 you are a member of yade-users, which 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 #250783]: some questions about 'O.dt' and 'PyRunner'
Question #250783 on Yade changed: https://answers.launchpad.net/yade/+question/250783 Hongyang Alex Cheng requested more information: I am curious that isn't the GlobalStiffnessTimeStepper doing well for making simulation stable? Why do you want a dynamic dt? Could you explain some more? Cheng -- You received this question notification because you are a member of yade-users, which 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 #250783]: some questions about 'O.dt' and 'PyRunner'
Question #250783 on Yade changed: https://answers.launchpad.net/yade/+question/250783 Hongyang Alex Cheng proposed the following answer: I am curious that isn't the GlobalStiffnessTimeStepper doing well for making simulation stable? Why do you want a dynamic dt? Could you explain some more? Cheng -- You received this question notification because you are a member of yade-users, which 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 #250783]: some questions about 'O.dt' and 'PyRunner'
Question #250783 on Yade changed: https://answers.launchpad.net/yade/+question/250783 Hongyang Alex Cheng proposed the following answer: For details on dynDt, please refer to [1]: https://yade- dem.org/doc/yade.wrapper.html?highlight=dt#yade.wrapper.Omega.dynDt -- You received this question notification because you are a member of yade-users, which 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 #246502]: Build remote interaction for WireMat
Question #246502 on Yade changed: https://answers.launchpad.net/yade/+question/246502 Hongyang Alex Cheng confirmed that the question is solved: Thanks Jérôme Duriez, that solved my question. -- You received this question notification because you are a member of yade-users, which 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 #246502]: Build remote interaction for WireMat
Question #246502 on Yade changed: https://answers.launchpad.net/yade/+question/246502 Status: Answered => Solved Hongyang Alex Cheng confirmed that the question is solved: Hi, jduriez Thanks for the remarks. I am simulating uniaxial compression on a woven textile sand bag. I choose WireMat for modeling the bag and FrictMat for the soil inside the bag. Here is my script for the engine parts (sorry for not pasting the whole script, it's too long). My script worked this time. Problem seems to be that I built remote links before defining the interaction loop. Anyway, thank you very much. ## ## Defining engines## ## O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(label='Sphaabb'),Bo1_Box_Aabb(label='Wallaabb')]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(label='Ig2ssGeom'),Ig2_Box_Sphere_ScGeom(label='wall')], [Ip2_WireMat_WireMat_WirePhys(label='wovenbag'), Ip2_FrictMat_FrictMat_MindlinPhys( krot=kr, ktwist=kr, label='soilsphere')], [Law2_ScGeom_WirePhys_WirePM(label='Law_wvbg'), Law2_ScGeom_MindlinPhys_Mindlin( includeMoment=True, label='Law_ssph')] ), NewtonIntegrator(damping=0.2,gravity=[0,0,-9.81],label='newton'), GlobalStiffnessTimeStepper(timeStepUpdateInterval=100) ] -- You received this question notification because you are a member of yade-users, which 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 #246502]: Build remote interaction for WireMat
Question #246502 on Yade changed: https://answers.launchpad.net/yade/+question/246502 Status: Solved => Open Hongyang Alex Cheng is still having a problem: Hi Klaus, I still have some questions about WireMat-FrictMat interaction. I understand WireMat inherits FrictMat so that CundallStrack law can work for the interactions between WireMat and FrictMat. What if I want to use Hertz-Mindlin law instead? I tried to do that by simply replacing the lp2 and law2 functors by Ip2_FrictMat_FrictMat_MindlinPhys and Law2_ScGeom_MindlinPhys_Mindlin, but it didn't work and gave me the following error report: FATAL FATAL /home/alexcheng/myYade/trunk/pkg/common/InteractionLoop.cpp:142 action: None of given Law2 functors can handle interaction #19133+19132, types geom:ScGeom=1 and phys:FrictPhys=3 (LawDispatcher::getFunctor2D returned empty functor) Could you please give some advices on this? btw, what do geom:ScGeom=1 and phys:FrictPhys=3 mean? Regards, Alex -- You received this question notification because you are a member of yade-users, which 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 #247095]: filling STL surface with spheres
Question #247095 on Yade changed: https://answers.launchpad.net/yade/+question/247095 Hongyang Alex Cheng proposed the following answer: Hi, Mayank Here is my approach. Hope it could give you some hint. I used autocad to built the predicate and output it as an stl file. Using 'stl2gts' my stl file was converted to the gts file which can be imported into Yade. As for preparing spheres using predicate as boundary, you can always use the function 'pack.filterSpherePack()'. Bruno Chareyre's suggestion in [1] will definitely help. [1] https://answers.launchpad.net/yade/+question/246971 Regards, Alex On Thu, Apr 17, 2014 at 8:01 AM, Jan Stránský < question247...@answers.launchpad.net> wrote: > Question #247095 on Yade changed: > https://answers.launchpad.net/yade/+question/247095 > > Status: Open => Answered > > Jan Stránský proposed the following answer: > Hi Mayank, > > your approach is correct. STL could not be used directly in Yade as a > predicate (to be filled by spheres), but can be converted to gts and with > gts it is possible. Are you sure you use stl2gts correctly? Also have a > look at discussion [1]. > cheers > Jan > > [1] https://answers.launchpad.net/yade/+question/246833 > > > > 2014-04-16 20:26 GMT+02:00 Mayank Sabharwal < > question247...@answers.launchpad.net>: > > > New question #247095 on Yade: > > https://answers.launchpad.net/yade/+question/247095 > > > > I am trying to fill spheres into a STL surface. > > My only objective is to use this STL as a predicate to be able to use the > > different sphere packings. > > Is it possible with the current yade? > > > > I tried converting the stl to gts but for some reason the stl2gts tool is > > not working (it doesnot give errors but doesnot give any output either) > > > > It would be great help if someone could help me out here > > > > thanks, > > Mayank > > > > -- > > You received this question notification because you are a member of > > yade-users, which 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 you are a member of > yade-users, which 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 you are a member of yade-users, which 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 #246971]: Attain a specified porosity in a predicate
Question #246971 on Yade changed: https://answers.launchpad.net/yade/+question/246971 Status: Answered => Solved Hongyang Alex Cheng confirmed that the question is solved: Thanks Bruno Chareyre, that solved my question. -- You received this question notification because you are a member of yade-users, which 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 #246971]: Attain a specified porosity in a predicate
Question #246971 on Yade changed: https://answers.launchpad.net/yade/+question/246971 Hongyang Alex Cheng posted a new comment: Thanks, Bruno. I take your advice and it actually worked excepted that the porosity is slightly different from I want, probably because of the rebound of the particles. I have been thinking that the function 'pack.filterSpherePack' does nothing but clip from the packing a sample using the predicate as boundary. Therefore the stress state inside the predicate must be uniform. Do you think in a geotextile bag the stress dis-uniformity in soil is negligible? If taking that into account, how to adjust the packing process? It is getting beyond the scope of this topic. Sorry about that. Alex Cheng -- You received this question notification because you are a member of yade-users, which 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 #246971]: Attain a specified porosity in a predicate
New question #246971 on Yade: https://answers.launchpad.net/yade/+question/246971 Hi, guys How can I add sphere packing in a predicate with a specified porosity? I've been trying to model a woven geotextile bag filled with Toyoura sand inside (standard sand in Japan). I built a predicate for the geotextile bag and imported it into yade-git, constructing a closed triagulated surface. To create fill the predicate with spheres in a given particle size distribution. I use the following commands. # The packing made by makeCloud is bigger in size than the predicate, in order to use the module pack.filterSpherePack(). sp = pack.SpherePack() sp.makeCloud(minCorner, maxCorner, rMean=-1, rRelFuzz=0.3, num=1000, perodic=False, porosity=0.95, psdSizes=[.1, .15, .2, .25, .4], psdCumm=[0, .1, .8, .96, 1], seed=1) pack.filterSpherePack(pred, sp, True) By now I have spheres in the predicate. I guess next my move is to expand the spheres radii and reducing frictional angle as the REFD method does in the triax-tutorial. But I think I can't use the triax engine in my case. Because there shouldn't be any confining pressure applied on the spheres, the internal compaction must be disabled. Could anyone give me some comments and advices? -- You received this question notification because you are a member of yade-users, which 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 #246502]: Build remote interaction for WireMat
Question #246502 on Yade changed: https://answers.launchpad.net/yade/+question/246502 Status: Answered => Solved Alex confirmed that the question is solved: Thank you Klaus. Looking forward to seeing your recent progress on Yade-DEM development in the upcoming Yade seminar. Regards, Alex -- You received this question notification because you are a member of yade-users, which 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 #246502]: Build remote interaction for WireMat
Question #246502 on Yade changed: https://answers.launchpad.net/yade/+question/246502 Status: Answered => Open Alex is still having a problem: Hi Klaus Thank you for your detailed explanation. I have some other questions. In your example wirecontacttest.py, I notice that there are two Ip2 functor and two Law2 functors defined for the WireMat and the FrictMat (which obeys Cundall’s constitutive law) respectively. Is that correct that the Ip2 and Law2 functors of FrictMat also work for the interaction between the dropping sphere block and the wire net? -- You received this question notification because you are a member of yade-users, which 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 #246507]: yade.wrapper.hideBody
Question #246507 on Yade changed: https://answers.launchpad.net/yade/+question/246507 Status: Open => Answered Alex proposed the following answer: Hi Hicham The following commands from yade import qt qt.Renderer().hideBody(id) might help. regards, alex chyalexch...@gmail.com From: Hicham BENNIOU Date: 2014-04-04 18:31 To: yade-users Subject: [Yade-users] [Question #246507]: yade.wrapper.hideBody New question #246507 on Yade: https://answers.launchpad.net/yade/+question/246507 Hello everyone, I was looking for hideBody [1] in OpenGLRenderer class, but it appears that it doesn't exist : Yade [1]: yade.wrapper.OpenGLRenderer.hideBody(0) --- AttributeErrorTraceback (most recent call last) /usr/bin/yadedaily in () > 1 yade.wrapper.OpenGLRenderer.hideBody(0) AttributeError: type object 'OpenGLRenderer' has no attribute 'hideBody' Help ? Thanks ! Hicham [1] : https://yade- dem.org/doc/yade.wrapper.html?highlight=hidebody#openglrenderer -- You received this question notification because you are a member of yade-users, which 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 you are a member of yade-users, which 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 #246502]: Build remote interaction for WireMat
New question #246502 on Yade: https://answers.launchpad.net/yade/+question/246502 Hi, In the example wiretensiltest.py, if I am not mistaken, the purpose of setting enlarge interaction radius is to build remote contact between wire spheres. I want to know which functor in the interaction loop is doing that. If it is the Ig2 functor, will it also work for the spheres of the non-WireMat material types while creating links? what is the attribute linkThresholdInteraction doing in Ip2 and Law2 functors? What will happen if its value is not one as default? -- You received this question notification because you are a member of yade-users, which 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