Re: [Yade-users] [Question #657232]: Need suggestions on the computer server for YADE computation.
Question #657232 on Yade changed: https://answers.launchpad.net/yade/+question/657232 Status: Open => Answered Christian Jakob proposed the following answer: Hi, I would definitely prefer more single CPU power and less number of CPUs because DEM is bad to parallelize (in general, not just YADE). Regards, Christian -- 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 #644873]: getroundness
Question #644873 on Yade changed: https://answers.launchpad.net/yade/+question/644873 Christian Jakob posted a new comment: Hi, I developed the code for getRoundness function. It should work for any mixed clump-sphere-packing. You can have a look at the source code [1]. I dont know why it is not working in your case. Can you provide a minimal working example [2]? [1] https://github.com/yade/trunk/blob/master/py/wrapper/yadeWrapper.cpp#L355 [2] https://en.wikipedia.org/wiki/Minimal_Working_Example -- 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 #599179]: Capillary contact law
Question #599179 on Yade changed: https://answers.launchpad.net/yade/+question/599179 Christian Jakob posted a new comment: Hi, I think the best way to understand capillary effect on particle forces is a force-distance-plot. Fortunately I have one in my PhD. You can have a look at page 35 (pic. 3.10). http://nbn-resolving.de/urn:nbn:de:bsz:14-qucosa-218045 btw, the plot was done with output from Law2_ScGeom_CapillaryPhys_Capillarity ;) Regards, Christian -- 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 #536811]: clump generation for sphere larger than one given size
Question #536811 on Yade changed: https://answers.launchpad.net/yade/+question/536811 Status: Open => Answered Christian Jakob proposed the following answer: hi, using O.bodies.replaceByClumps() within a loop is not a good idea. to solve your problem i see two ways: 1.) you generate particles with radius < 0.01 first, then you run replaceByClumps() and generate the rest of the particles afterwards 2.) you can modify replaceByClumps() function in c++ [1], so that you can exclude particles with an excludeList like in updateClumpProperties() [2]. [1] https://github.com/yade/trunk/blob/master/py/wrapper/yadeWrapper.cpp#L210 [2] https://yade-dem.org/doc/yade.wrapper.html?highlight=updateclump#yade.wrapper.BodyContainer.updateClumpProperties -- 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 #428045]: introduction
Question #428045 on Yade changed: https://answers.launchpad.net/yade/+question/428045 Status: Open => Answered Christian Jakob proposed the following answer: Hi You just have to type "yade tutorial" into google... Regards, Christian -- 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 #402952]: Using a gts surface as a shell and filling it using radius expansion method
Question #402952 on Yade changed: https://answers.launchpad.net/yade/+question/402952 Christian Jakob proposed the following answer: I am not familiar with polyhedrons, but examples can be found here: https://github.com/yade/trunk/tree/master/examples/polyhedra -- 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 #402952]: Using a gts surface as a shell and filling it using radius expansion method
Question #402952 on Yade changed: https://answers.launchpad.net/yade/+question/402952 Status: Open => Needs information Christian Jakob requested more information: Hi, Do you speak about clumps (rigid/non-deformable bodies consisting of several spherical particles) or about clusters (deformable bodies consisting of several spherical particles bonded together via cohesive contacts)? Yade can handle both, but generation method may be different. Please also note, that Yade can also simulate polyhedral particles like shown here: http://dx.doi.org/10.1016/j.powtec.2014.05.052 Regards, Christian -- 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 #394789]: Using radius expansion method to bulid a DEM model
Question #394789 on Yade changed: https://answers.launchpad.net/yade/+question/394789 Status: Open => Answered Christian Jakob proposed the following answer: Hi, Please have a look at the triax-tutorial. The TriaxialStressController provides the option "internalCompaction", which will do what you want. https://github.com/yade/trunk/blob/master/examples/triax-tutorial /script-session1.py#L94 Hope it helps, Christian -- 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 #384899]: VoxelPorosity
Question #384899 on Yade changed: https://answers.launchpad.net/yade/+question/384899 Status: Open => Answered Christian Jakob proposed the following answer: I investigated your script a little bit and tried to use another range in voxelPorosity() function: VP=utils.voxelPorosity(200,(-g,-g,-g),(g,g,g)) So now I get VP 0.757457375 which means the function is working. I am not sure what g and h means in your script, but it seems you just did a wrong input. Also I dont know if voxelPorosity() is a good choice in this case since you have a spherical shaped packing ... Christian -- 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 #384899]: VoxelPorosity
Question #384899 on Yade changed: https://answers.launchpad.net/yade/+question/384899 Status: Open => Answered Christian Jakob proposed the following answer: Hi, For voxelPorosity you need to specify a start point/vector and an end point/vector, not a volume. It will divide the cuboid range into res*res*res cells and counts cells inside and outside particles, where res is the resolution (default: 200). The function returns the ratio, which is porosity. https://yade- dem.org/doc/yade.utils.html?highlight=voxelporo#yade._utils.voxelPorosity Hope it helps, Christian -- 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 #371126]: particles fly out of boundary
Question #371126 on Yade changed: https://answers.launchpad.net/yade/+question/371126 Status: Open => Answered Christian Jakob proposed the following answer: I am not familiar with randomDensePack, but as the name indicates it will generate a dense packing. So maybe particles are overlapping and high forces will act on particles. If time step is too high particles can penetrate through walls (facets). If you want to calm your model down to adequate velocities you may also consider calm() function [1] in combination with a PyRunner[2]. An example can be found here: https://github.com/yade/trunk/blob/master/examples/clumps /replaceByClumps-example.py#L85 Regards, Christian [1] https://yade-dem.org/doc/yade.utils.html?highlight=calm#yade._utils.calm [2] https://yade-dem.org/doc/yade.wrapper.html?highlight=pyrunner#yade.wrapper.PyRunner -- 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 #371126]: particles fly out of boundary
Question #371126 on Yade changed: https://answers.launchpad.net/yade/+question/371126 Status: Open => Answered Christian Jakob proposed the following answer: Hi, Did you try to use a lower time step? -- 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 #204526]: delete particles with active capillary law
Question #204526 on Yade changed: https://answers.launchpad.net/yade/+question/204526 Christian Jakob posted a new comment: O.bodies.erase(id) is working fine in this case. -- 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 #295676]: Pack.inCylinder with specific clump
Question #295676 on Yade changed: https://answers.launchpad.net/yade/+question/295676 Status: Open => Answered Christian Jakob proposed the following answer: Hi, You can try replaceByClumps() method (see example in clumps folder). The procedure is to create spheres first and replace them with specific clumps (via clump templates) in a second step. The disadvantage of this method is, that there might be some overlaps between the clumps. So you maybe have to reduce overlaps via calm() function. Hope it helps, Christian -- 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 #294232]: change in volume of meniscus as a function of capillary pressure
Question #294232 on Yade changed: https://answers.launchpad.net/yade/+question/294232 Status: Open => Answered Christian Jakob proposed the following answer: I think you do not need i.phys.isBroken at all. What you need is at least 1 step between parameter change and measurement. Good luck -- 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 #294232]: change in volume of meniscus as a function of capillary pressure
Question #294232 on Yade changed: https://answers.launchpad.net/yade/+question/294232 Status: Open => Answered Christian Jakob proposed the following answer: Hi, O.engines = O.engines[:3]+[Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=cur_capillary,label='capillary')]+O.engines[3:] Maybe this line is the problem in your script (not tested, not sure what it does). If the engine label is 'capillary' you can try this line instead: capillary.capillaryPressure = cur_capillary -- 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 #294215]: how computer configuration to run discrete element method in YADE
Question #294215 on Yade changed: https://answers.launchpad.net/yade/+question/294215 Status: Open => Answered Christian Jakob proposed the following answer: Please have a look to the installation instructions: https://yade-dem.org/doc/installation.html#prerequisities -- 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 #291202]: Clumping - angl of repose
Question #291202 on Yade changed: https://answers.launchpad.net/yade/+question/291202 Status: Open => Answered Christian Jakob proposed the following answer: Hi, I did not investigate your script very deeply, but I guess the explosion arises from replaceByClumps() function. A workaround for avoiding explosions is given in the original example line 81 and following ... https://github.com/yade/trunk/blob/master/examples/clumps /replaceByClumps-example.py#L81 Cheers, Christian -- 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 #292663]: capillary pressure
Question #292663 on Yade changed: https://answers.launchpad.net/yade/+question/292663 Status: Open => Answered Christian Jakob proposed the following answer: Hi, I get the same message, when there is a contact between particles with radius quotient R1/R2 > 10. The solution for this contacts are simply out of data range in the capillary files. In my case this had nearly no influence on the results of my calculations. Regards, Christian -- 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 #279359]: MakeCloud() single-threaded?
Question #279359 on Yade changed: https://answers.launchpad.net/yade/+question/279359 Status: Open => Answered Christian Jakob proposed the following answer: Hi, makeCloud() is used for particle generation and has no influence on DEM calculation cycle. It is not parallelized atm [1]. Christian [1] https://github.com/yade/trunk/blob/2fd300ded158f3fe938ab6b7688aaff9b9ad35a6/pkg/dem/SpherePack.cpp#L86 -- 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 #279270]: Problem with efficiency
Question #279270 on Yade changed: https://answers.launchpad.net/yade/+question/279270 Status: Open => Needs information Christian Jakob requested more information: Hi, Does the generation of clumps takes longer in the 1.20 version, or does the calculation takes longer? Can you please check and report? Thank in advance, Christian -- 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 Christian Jakob proposed the following answer: If you really need to modify your clumps, I think O.bodies.releaseFromClump(0,3) O.bodies.addToClump([0],3) is the only way I know. -- 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 Christian Jakob posted a new comment: Hi, I do not understand the reason this: > ## here is the example: > a = utils.sphere([0,0,0],1) > b = utils.sphere([0,0,1],1) > c = utils.sphere([0,0,2],1) > O.bodies.appendClumped([a,b,c]); > > ##Now, I want to mannually change the position of 'a': > a.state.pos = Vector3(0,0,100) Why dont you do the following? a = utils.sphere([0,0,100],1) b = utils.sphere([0,0,1],1) c = utils.sphere([0,0,2],1) O.bodies.appendClumped([a,b,c]) -- 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 #270627]: translate problem
Question #270627 on Yade changed: https://answers.launchpad.net/yade/+question/270627 Christian Jakob posted a new comment: Hi, We (Yade community) really try to help people with answers on questions or solutions for problems. Unfortunately sometimes the questions can not be answered and/or the problems can not be solved. In your case the problem for us is to understand what is the question or where is the problem. With comments like there is something I don’t understand in the model behavior we can not help you, sorry. Additionally I think you did not read [1] carefully. There it is said: 3. Together with problem statement, send also a code which will let others reproduce your problem. Try to send it as a MWE (Minimal Working Example™). It should be: * minimal = short (it is difficult to find mistakes in long pieces of code and it is usually much more annoying). * working = before sending, test it to be sure it is working (it has no syntax errors etc.). Regards, Christian [1] https://yade-dem.org/wiki/Howtoask -- 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 #270747]: Clumps
Question #270747 on Yade changed: https://answers.launchpad.net/yade/+question/270747 Status: Open = Answered Christian Jakob proposed the following answer: Hi, I do not fully understand your question, but I guess you want to avoid the explosion of your model after the replacement step?! ... In this case, please have a look at the example line 81 and following [1]. replaceByClumps() can produce overlaps between newly created clumps. If you have a high stiffness, the model will explode. To avoid this you can use calm() function. calm() is used to set particle velocities and rotations to zero (see example). Along with PyRunner it will slowly reduce the overlaps. If your model still explodes, decrease iterPeriod of the PyRunner. If you want to replace just a specific part of your model you can try to generate this specific part first, then use replaceByClumps() and generate the rest of the model afterwards. Hope it helps, Christian [1] https://github.com/yade/trunk/blob/master/examples/clumps /replaceByClumps-example.py#L81 -- 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 #270488]: Clumping
Question #270488 on Yade changed: https://answers.launchpad.net/yade/+question/270488 Status: Open = Answered Christian Jakob proposed the following answer: Hi, Just a short answer regarding 1- and 2-: 1- Please have a look at [1]. There is also an example, see [2]. 2- I did not test your script. Maybe an indention problem? I hope someone else can help ... Regards, Christian [1] https://yade-dem.org/doc/yade.wrapper.html?highlight=replaceby#yade.wrapper.BodyContainer.replaceByClumps [2] https://github.com/yade/trunk/blob/master/examples/clumps/replaceByClumps-example.py -- 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 #269942]: How calculate the repose angle
Question #269942 on Yade changed: https://answers.launchpad.net/yade/+question/269942 Status: Open = Answered Christian Jakob proposed the following answer: I think you need to click on the mouse wheel ... (as mentioned in the script) -- 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 #269942]: How calculate the repose angle
Question #269942 on Yade changed: https://answers.launchpad.net/yade/+question/269942 Status: Open = Answered Christian Jakob proposed the following answer: Hi, I had the same problem some years ago. I used pythons scatter function for this purpose. You need to give some input data at the beginning of the script and it should work (not tested). A scatter-plot window should open and you can get the angle by clicking at two points in the window. Here is the script: #!/usr/bin/python # -*- coding: utf-8 -*- center_model = Vector3(0,0,0) #put the center of your model here O.load(yourSaveFile) #put the name of your save file here outputFilename = 'yourOutputName.out' #put your output filename here height = [] rad_dist = [] for ii in range(0,len(O.bodies)): b = O.bodies[ii] if b and isinstance(b.shape,Sphere): pos = b.state.pos dist_vec = pos - center_model rad_dist.append(sqrt( (dist_vec[0])**2 + (dist_vec[1])**2 ))#get radial distance from center axis height.append(pos[2]) #get height and write into list from pylab import * scatter(rad_dist,height,marker='+') data_for_angle = ginput(0,0)#get the coordinates by clicking two points, then click on the mouse wheel zdiff = abs(data_for_angle[0][1] - data_for_angle[1][1])#y entries of picked points (respectively spheres height) rdiff = abs(data_for_angle[0][0] - data_for_angle[1][0])#x entries of picked points (respectively spheres radial distance from center) alpha = atan(zdiff/rdiff) alpha = 180*alpha/math.pi #in [°] ''' z |- d | - i |- f |( -- tan(alpha) = zdiff / rdiff f -- rdiff ''' # print result and write output: print 'angle of repose is ',alpha,' degree' f = open(outputFilename,'a') f.write('result for %s: %.2f degree\n' % (part,alpha)) f.close() exit() -- 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 #269912]: Problem with ClumpTemplate()
Question #269912 on Yade changed: https://answers.launchpad.net/yade/+question/269912 Status: Open = Answered Christian Jakob proposed the following answer: Hi, Did you had a look at the example [1] after line 83? HIH, Christian [1] https://github.com/yade/trunk/blob/master/examples/clumps /replaceByClumps-example.py#L83 -- 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 #253562]: Status of bug #1362090 changed to 'Fix Released' in Yade
Bug #1362090 status changed in Yade: Confirmed = Fix Released https://bugs.launchpad.net/yade/+bug/1362090 boundary thickness is 0.001 by default in FlowEngine This bug is linked to #253562. PFV compressibility in undrained triaxial test https://answers.launchpad.net/yade/+question/253562 -- 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 #261529]: a cylinder as a single body
Question #261529 on Yade changed: https://answers.launchpad.net/yade/+question/261529 Christian Jakob posted a new comment: Hi, updateProperties function assumes, that clump members are spheres. it will not work for facets atm ... regards, christian -- 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 #258987]: updating O.bodies
Question #258987 on Yade changed: https://answers.launchpad.net/yade/+question/258987 Christian Jakob proposed the following answer: Hi, There is already a calm function in yade, that you can use: https://yade-dem.org/doc/yade.utils.html?highlight=calm#yade._utils.calm Regards, Christian -- 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 #258859]: Clump mass issue with replaceByClumps()
Question #258859 on Yade changed: https://answers.launchpad.net/yade/+question/258859 Christian Jakob proposed the following answer: Hi, Thank you for your feedback on yade's clump logic. Your results look good even for a low discretization (in your case 6) ... updateClumpProperties() can be used when the mass/volume/inertia of a clump was changed manually by the user (and it should be updated to the real value afterwards). I think you do not need it atm. But of course you can check this function ... Cheers, Christian -- 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 #258859]: Clump mass issue with replaceByClumps()
Question #258859 on Yade changed: https://answers.launchpad.net/yade/+question/258859 Status: Open = Answered Christian Jakob proposed the following answer: Hi, For inertia/mass/volume calculation, use O.bodies.replaceByClumps(templates,[1.0], discretization=5) where discretization is explained here: https://yade- dem.org/doc/yade.wrapper.html#yade.wrapper.BodyContainer.clump So for more accurate inertia/mass/volume calculation increase discretization. But however increasing it will lead to higher initialization time for clumps. Fortunately it is parallelized, so you can use -j option to reduce calc. time on creation of clumps ;) Please give us the results for mass for discretization=5, 10 and 15. Thanks. Regards, Christian -- 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 #257388]: Creating a pack with Nc more than 6
Question #257388 on Yade changed: https://answers.launchpad.net/yade/+question/257388 Christian Jakob posted a new comment: Hi, You can use clumps instead of spheres to increase number of contacts. But the packing will not be denser by doing so... Christian -- 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 #257388]: Creating a pack with Nc more than 6
Question #257388 on Yade changed: https://answers.launchpad.net/yade/+question/257388 Christian Jakob posted a new comment: Of course compacting spheres before clumping leads to arbitrary large Nc depending on the number of spheres per clump, but I'm not really sure that clumps will give Nc5 if they are compacted as clumps (I mean starting from a clump cloud). Do you guys have results about that? Yes, I see, that the difference between 100% spheres and 100% clumps (two spherical members) is not much (around +0.1 ... +0.2). On one hand there are more contacts, but on the other hand the porosity is higher and hence less contacts. A denser packing can be achieved by using lower friction while consolidation. Also using more members per clump will lead to more contacts, but I did notdo this, yet. -- 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 #255486]: Confusion about tutorial: creating and running simulation
Question #255486 on Yade changed: https://answers.launchpad.net/yade/+question/255486 Status: Open = Answered Christian Jakob proposed the following answer: Hi Robert, In 'Creating simulation', why does the example include an error message? I don't know. The command is working at my machine, so the problem is on buildbot side ... You can ignore the error. And in the next section, shouldn't 'As explained above' be 'As explained below'? fixed ;) If you find any more mistakes in the documentation, please give us a note and it will be fixed. Regards, Christian -- 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 #255199]: CohFrictMat engine
Question #255199 on Yade changed: https://answers.launchpad.net/yade/+question/255199 Christian Jakob posted a new comment: hi, afaik you have set both factors to 1.2 (aabbEnlargeFactor and interactionDetectionFactor) for correct implementation. christian -- 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 #253562]: PFV compressibility in undrained triaxial test
Question #253562 on Yade changed: https://answers.launchpad.net/yade/+question/253562 Linked to bug: #1362090 https://bugs.launchpad.net/bugs/1362090 problem with pressure calculation in PFV cells -- 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 #253562]: PFV compressibility in undrained triaxial test
New question #253562 on Yade: https://answers.launchpad.net/yade/+question/253562 Hi, I want to get the pore pressure of a saturated soil sample in an undrained triaxial test when the sample volume is decreasing. ... TriaxialStressController(thickness=0,stressMask=7,internalCompaction=False,label='triax') ... triax.goal1 = 35000 triax.goal2 = 35000 triax.goal3 = 7 triax.wall_back_activated = False ... flow.viscosity = 1.3e-3 flow.fluidBulkModulus = 2e9 flow.bndCondValue=[0,0,0,0,0,0] flow.bndCondIsPressure=[0,0,0,0,0,0] ... For one of my simulations a decrease of ca. 1,8 % in porosity (=-dV/V) is obtained. When I calculate corresponding pressure increase I get dp = K*(-dV/V) = 2e9*0.018 = 3600 Pa But the output from PFV shows maximum average pressures around of 16 Pa (not MPa!). As long as I use SI units the output of flow.averagePressure() should be Pa, right? So is there a problem in the code or did I miss something in the compressibility option of PFV? regards, christian -- 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 #252374]: Porosity of a mixture of clumps and spheres
Question #252374 on Yade changed: https://answers.launchpad.net/yade/+question/252374 Status: Open = Answered Christian Jakob proposed the following answer: Hi, If you use makeCloud() it is recommended to call O.bodies.updateClumpProperties() It will correct volume/masses/inertia of clumps [1]. If you use TriaxialStressController() you get the porosity [2] very cheap via triax.porosity I updated the code, so that clumps volume is considered correctly. Another tool is voxelPorosity() It can also be used for your purpose [3], but every call costs a little bit computation time... Hope it helps, Christian [1] https://yade-dem.org/doc/yade.wrapper.html?highlight=updateclump#yade.wrapper.BodyContainer.updateClumpProperties [2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.TriaxialStressController.porosity [3] https://yade-dem.org/doc/yade.utils.html?highlight=voxelporosity#yade._utils.voxelPorosity -- 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 #252063]: PFV with clumps?
Question #252063 on Yade changed: https://answers.launchpad.net/yade/+question/252063 Status: Answered = Solved Christian Jakob 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 #252053]: Yade in PBS
Question #252053 on Yade changed: https://answers.launchpad.net/yade/+question/252053 Christian Jakob posted a new comment: Hi Carl, Nice to hear that you get yade running on your cluster. We are using yade on HPC in Freiberg too. It is hard to say what the problem is. Your pbs script looks fine. But you should check where the problem is: * error in the script? * problem with your yade version? * problem on the cluster (missing/wrong libraries, compilation problem, ...)? * problem with pbs? Try to run the script on another computer (same yade version) Try to run the script on the cluster without pbs Try to run the script with another yade version Good luck! Christian -- 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 #252063]: PFV with clumps?
New question #252063 on Yade: https://answers.launchpad.net/yade/+question/252063 Hi, I am thinking about a simulation with clumps and PFV method. I know clumps are not compatible with PFV, yet. I simulate a sand at microscale and want to have a more realistic grain shape using clumps (maybe 2 or more overlapping members). I see there is the first problem in triangulation algorithm. But it should be possible to use clump center and equivalent radius instead of sphere center and sphere radius for triangulation. I think this is ok as long as clumps are nearly spherical, right? The next problem is the calculation of fluid volume inside the cells. But isn't it possible to use clump centers and equivalent radius at this stage as well? I know the volumes will not be exactly, but I could live with a small error as long as calculation stays stable... What do you think about this? Did I miss something, that makes it completely impossible? Regards, Christian -- 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 #250157]: minimum computer configuration
Question #250157 on Yade changed: https://answers.launchpad.net/yade/+question/250157 Christian Jakob proposed the following answer: Hi Anna, With my experience I can say, that for DEM calculations following specifications are important: - high CPU frequency - disable multithreading (just use physical cores) RAM depends on what you need (more particles = more RAM). 16 GB should be enough ... I prefer Intels i7 technology. I am working with an i7-2600K (4 physical cores with 3,4 GHz) and very happy with this. Regards, Christian -- 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 #250196]: undrained triax test of fully saturated material (using TriaxialStressController and FlowEngine)
New question #250196 on Yade: https://answers.launchpad.net/yade/+question/250196 Hi, For my simulations it is no problem to do a drained triax test with this PFV options: flow.bndCondIsPressure=[0,0,1,0,0,0] flow.bndCondValue=[0,0,0,0,0,0] I also want to do an undrained test, like this: flow.bndCondIsPressure=[0,0,0,0,0,0] flow.bndCondValue=[0,0,0,0,0,0] Then I get this warning: CHOLMOD warning: matrix not positive definite something went wrong in Cholesky factorization, use LDLt as fallback this time To solve this, I have to use at least one fixed pressure boundary condition or impose pressure at one (or more) cells, e.g.: flow.imposePressure(Vector3(triax.width/2,triax.height/2,triax.depth/2),0.001) So the pressure of this cell will be constant during the calculation and in the end it is still a drained test, right? If I do flow.clearImposedPressure() after some steps were done, the same warning (see above) appears again. Is there a trick to simulate an undrained test with FlowEngine? Christian -- 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 #250196]: undrained triax test of fully saturated material (using TriaxialStressController and FlowEngine)
Question #250196 on Yade changed: https://answers.launchpad.net/yade/+question/250196 Christian Jakob posted a new comment: Hi Jerome, Ok, I see the problem is that FlowEngine supports just incompressible flow. Hm, then the question is how to make FlowEngine compressible? Christian -- 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 #249457]: Results
Question #249457 on Yade changed: https://answers.launchpad.net/yade/+question/249457 Status: Open = Answered Christian Jakob proposed the following answer: Hi Diego, For storing values, you can create a txt file and write data as ascii in it: f = open('myFile.txt','w') #'w' = write, 'a' = append, 'r' = read f.write('%i %e %f' % (intValue,floatValue1,floatValue2)) f.close() There are several other possibilities like tables, please have a look at python manual: https://docs.python.org/2/tutorial/inputoutput.html Hope it helps, Christian -- 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 #247705]: PFV problems with sedimentation in fluid
Question #247705 on Yade changed: https://answers.launchpad.net/yade/+question/247705 Christian Jakob posted a new comment: hi, just a short answer to question (2): stressMask is not good explained in documentation, have a look at examples/triax-tutorial for a better documentation https://github.com/yade/trunk/blob/master/examples/triax-tutorial /script-session1.py#L87 hope it helps, christian -- 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 #247158]: Yade rookie
Question #247158 on Yade changed: https://answers.launchpad.net/yade/+question/247158 Status: Open = Answered Christian Jakob proposed the following answer: Hi rookie ;) Assuming you use facets, FrictMat and ScGeom: It seems you forgot to add Ig2_Facet_SphereScGeom() to InteractionLoop engine. Just have a look at this code from examples/funnel.py: O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], ... ), ... ] If are using boxes, just replace all Facet by Box in the code above. Good luck, Christian -- 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 #245755]: Setting manually the inertia and mass center of clumps
Question #245755 on Yade changed: https://answers.launchpad.net/yade/+question/245755 Christian Jakob posted a new comment: Hi Pascal, I can not really answer your question, but here some notes: - If you use complex clump geometries it is recommended to use the actual Yade version since there some bug fixes were done in clump logic a few weeks ago. - It would be very nice to get a verification if inertia calculation procedure works fine in Yade ... Regards, Christian -- 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 #245755]: Setting manually the inertia and mass center of clumps
Question #245755 on Yade changed: https://answers.launchpad.net/yade/+question/245755 Christian Jakob posted a new comment: Do you use actual for current also in german? ;) right, I am speeking germish, not english ^^ -- 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 #245072]: exciting box
Question #245072 on Yade changed: https://answers.launchpad.net/yade/+question/245072 Christian Jakob posted a new comment: Hi, I am not sure, but I think you forgot to set the materials to your bodies (spheres and facets). Changing these line may solve it (not tested): Mat = O.materials.append(FrictMat(young=30e9,density=1000,poisson=.2,frictionAngle=.5)) O.bodies.append(utils.geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=63,material=Mat)) O.bodies.append([sphere([0.5,0.5,0.5],.05,material=Mat)]) Regards, Christian -- 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 #245072]: exciting box
Question #245072 on Yade changed: https://answers.launchpad.net/yade/+question/245072 Christian Jakob posted a new comment: sorry, wrong code. it has to be like... Mat = FrictMat(...) O.materials.append(Mat) -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Christian Jakob posted a new comment: What does it do the new feature? Well for your script replacing 50 spheres by clumps I got a generation time about 35 seconds. The new feature allows you to run it at all your cpu cores, which leads to less generation time. For your script with 50 particles I achieved approx. 8 seconds using 8 cores (command: yade -j8 scriptname.py). Imagine you want to create a model consisting of 1000 particles or more ... cheers, christian -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Christian Jakob posted a new comment: Hi Behzad, I hope you enjoy this new feature too ;) https://github.com/yade/trunk/commit/828893babba373d7c376c084a288e6af2f910512 Cheers, Christian -- 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 #244634]: only 1 core working out of 32
Question #244634 on Yade changed: https://answers.launchpad.net/yade/+question/244634 Christian Jakob proposed the following answer: Definitely better! But still only 8 out of 16 CPU are used, between 30 and 50%. Instead 8 cores were never used for the entire time that yade was running. You have 8 physical cores, which can be hyperthreaded to achieve 16 virtual cores. Openblas is still installed, but now HyperThreading is disabled. So now you can only use your 8 physical cores (most efficient for DEM, perfect). Definitely better! But still only 8 out of 16 CPU are used, between 30 and 50%. Instead 8 cores were never used for the entire time that yade was running. Is there anything else I could do? Nope, you did it, everything is fine. cheers, christian -- 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 #244634]: only 1 core working out of 32
Question #244634 on Yade changed: https://answers.launchpad.net/yade/+question/244634 Christian Jakob proposed the following answer: forgot to say that you should not use a higher value than 8 in the -j option: yade-Unknown -j8 checkPerf.py -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Status: Open = Needs information Christian Jakob requested more information: setting integrateInertia=False doesns't fix the problem. confirmed One problem I detected with my own script is, that I don't think it's the same problem. Because the number of balls in the clump I have is not that much. I have around 30 balls for clumps in my work. But, Francois is talking about 256 balls. ok Christian and Klaus, what is the exact version of yade and ubuntu you guys are using and do not have any issue which the script I posted? I am using latest trunk version (development version) [1]. Latest change in clump logic has been done in october 2013 [2]. So there was no change for versions = 1.0.0 ... [1] https://github.com/yade/trunk [2] https://github.com/yade/trunk/commit/8202468e7b483c936dbd099b2b03b79ee7150ab4 how do you get your script working? can you post a combination of relRadList and relPosList, that is working? christian -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Christian Jakob posted a new comment: sorry forgot to end my sentence: One problem I detected with my own script is, that when using replaceByClumps in a periodic space sometimes balls are created outside periodic cell. To fix this one has to decrease generation area for sphere before using replacebyClumps, but it is not the problem here... -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Christian Jakob posted a new comment: when using replaceByClumps in a periodic space sometimes balls are created outside periodic cell It should not be a problem. Is it? In yade's PBC particles can be everywhere. I do not remember the error message, but it is breaking the code somewhere. @behzad: While investigating the code I found a bug, but not the one I was looking for ^^ I propose to continue this conversation at the bug report, I am quite sure it is related somehow: https://bugs.launchpad.net/yade/+bug/1273172 -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Status: Needs information = Answered Christian Jakob proposed the following answer: Hi again, I have a good and a bad message for you: The good one is: I found the problem and fixed it [1]. The bad one is: For multiple overlapping clump members an additional space grid integration is needed for getting volume of the clump templates. It is done once for each template. I checked the volumes of spheres and corresponding clumps (using discretization=5) and it is ok now: Volume of shpere before replacement: 0.110894 Volume of clump: 0.110879 Volume of shpere before replacement: 0.106574 Volume of clump: 0.106577 Volume of shpere before replacement: 0.102368 Volume of clump: 0.102379 ... While deactivating integration scheme (integrateInertia=False) properties (volume/mass/inertia) are not valid, but initialization is much faster: Volume of shpere before replacement: 0.110894 Volume of clump: 0.330293 Volume of shpere before replacement: 0.106574 Volume of clump: 0.317426 Volume of shpere before replacement: 0.102368 Volume of clump: 0.304897 ... Have fun with this new feature! Now I am thinking about a parallelizing integration scheme ... Cheers, Christian Oh I forgot: It does not fix the other bug [2]. [1] https://github.com/yade/trunk/commit/f60e0c0315595842c853a0f154121a5d4e1449a2 [2] https://bugs.launchpad.net/yade/+bug/1273172 -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Christian Jakob proposed the following answer: Would you please let me know how can I update my code to get change you made? I appreciate it. There are different possibilities: https://yade-dem.org/doc/installation.html 1) The fastest option is to download trunk.zip at github https://github.com/yade/trunk unzip, compile, use 2) If you use yade-daily, the update should be available in the coming hours/days/weeks? I do not know... 3) Wait for the next official release ... Cheers, Christian -- 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 #244312]: Clump templates
Question #244312 on Yade changed: https://answers.launchpad.net/yade/+question/244312 Christian Jakob posted a new comment: Hi Behzad, Some time ago Francois Kneib detected this bug in clump generation: https://bugs.launchpad.net/yade/+bug/1273172 Please have a look at this. Is it the same in your case? Christian -- 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 #243019]: loop over cells in PFV using id
Question #243019 on Yade changed: https://answers.launchpad.net/yade/+question/243019 Christian Jakob posted a new comment: There is no python wrapping of the cell objects , hence we need one getter per variable. wouldnt it be better create wrapper for cell objects, so that users can do something like: for c in O.flowCells: print c.pressure c -- 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 #243019]: loop over cells in PFV using id
Question #243019 on Yade changed: https://answers.launchpad.net/yade/+question/243019 Status: Answered = Solved Christian Jakob 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
[Yade-users] [Question #243019]: loop over cells in PFV using id
New question #243019 on Yade: https://answers.launchpad.net/yade/+question/243019 Hi, How can I loop over all cells in PFV method? I noticed there is an id of each cell, but how can I find out all ids? -- 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 #226679]: Clump with wall
Question #226679 on Yade changed: https://answers.launchpad.net/yade/+question/226679 Christian Jakob proposed the following answer: But the Programme stop and warnt NameError: name 'makeClumpTemplate' is not defined. The command is clumpTemplate(), not makeClumpTemplate() ... hih c -- 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 #240620]: tetra +facet
Question #240620 on Yade changed: https://answers.launchpad.net/yade/+question/240620 Christian Jakob posted a new comment: Polyhedra code was implemented with 3715, so the examples should work at 3730 http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/revision/3715 -- 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 #240557]: error during make install
Question #240557 on Yade changed: https://answers.launchpad.net/yade/+question/240557 Christian Jakob posted a new comment: I try to install Yade on the high performance cluster of our university. I like to hear that. Do you really need GUI at the cluster? Does this cluster provide GUI applications? I found this mail in the archive, hope it helps ... http://www.mail-archive.com/yade-users@lists.launchpad.net/msg05936.html christian -- 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 #240620]: tetra +facet
Question #240620 on Yade changed: https://answers.launchpad.net/yade/+question/240620 Christian Jakob posted a new comment: hi martin, there are several examples in examples/tetra and examples/polyhedra folders. good luck, christian -- 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 #240628]: get fluid velocity in PFV cells
New question #240628 on Yade: https://answers.launchpad.net/yade/+question/240628 Hi, I am using PFV and want to know how I can get fluid velocity in the cells. While visualisation in paraview one can see fluid velocities, but how can i get this information in python code? I tried following codes: 1) vel = flow.avFlVelOnSph(b.id) NameError: global name 'avFlVelOnSph' is not defined Maybe this happens, because I am using an old yade version from 16.10.2013 ...? 2) vel = flow.getCellFlux(flow.getCell(Vector3(xDim/2,yDim/2,zDim/2))) ERROR /home/christian/YADE/trunk/pkg/dem/FlowEngine.cpp:189 getCellFlux: Getting flux with cond higher than imposedP size. 0.0 Is this just for cells, with imposed pressure? 3) vel = flow.getBoundaryFlux(0) This one is working, but I get much smaller values, than in Paraview (factor 1e6!) christian -- 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 #238797]: Update interaction physics after changing material of spheres
Question #238797 on Yade changed: https://answers.launchpad.net/yade/+question/238797 Christian Jakob proposed the following answer: You can avoid explosions by using calm() function, which sets velocities of all bodies to zero. In addition with a pyrunner, it can help you ... You can try something like this: O.engines=O.engines+[PyRunner(iterPeriod=1,command='calm()',label='calmRunner')] O.run(1000,True) calmRunner.iterPeriod = 9 O.run(5000,True) calmRunner.iterPeriod = O.run(1,True) calmRunner.iterPeriod = 9 O.run(2,True) calmRunner.dead=True#deactivate calm runner -- 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 #233232]: F
Question #233232 on Yade changed: https://answers.launchpad.net/yade/+question/233232 Christian Jakob posted a new comment: Finally I could fix the bouyancy example. The problem was in saturatedList: When a clump was added to this (python) list, its permanent forces added by addF() disappeared. When a sphere was added to this (python) list, everything works fine ... I dont know why, but it is fixed now by removing saturatedList from the script ^^ Christian -- 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 #233320]: compiling with quad precision
Question #233320 on Yade changed: https://answers.launchpad.net/yade/+question/233320 Christian Jakob posted a new comment: maxF = 0.0 minF = 1e9 for b in O.bodies: F = O.forces.f(b.id) maxF = max(maxF,F.norm()) minF = min(minF,F.norm()) -- 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 #233320]: compiling with quad precision
Question #233320 on Yade changed: https://answers.launchpad.net/yade/+question/233320 Christian Jakob posted a new comment: std 0,0313071957 0,0460952889 hi matthias, thanks for this interesting results. that means quad precision shows nearly no effect on standard deviation in multithreading? amazing! Cheers, Christian -- 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 #233320]: compiling with quad precision
Question #233320 on Yade changed: https://answers.launchpad.net/yade/+question/233320 Christian Jakob posted a new comment: this standard deviation is the standard deviation of the posted data columns. yes, i followed this discussion from the beginning. so of course, also long double multithreaded simulation is affected by floating point indeterminism, but fewer than double or single floats right. i meant, that the difference 0,046 - 0,031 = 0,015 is not as high as expected (at least by me). could you (as far as your time allows) please check indeterminism differences (double vs. long double) with maximum/minimum body forces at particles? e.g. at end of your simulation ... -- 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 #233232]: F
Question #233232 on Yade changed: https://answers.launchpad.net/yade/+question/233232 Christian Jakob posted a new comment: Thanks Bruno for this new feature. Unfortunately it is not possible to use O.forces.addF() for clumps ... -- 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 #233232]: F
Question #233232 on Yade changed: https://answers.launchpad.net/yade/+question/233232 Christian Jakob posted a new comment: https://github.com/yade/trunk/blob/master/examples/clumps/apply- buoyancy-clumps.py replace lines 99 - 104 by this line: O.forces.addF(b.id,F_buo,permanent=True) and you will see, that buoyancy is not acting on clumps ... -- 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 #233232]: F
Question #233232 on Yade changed: https://answers.launchpad.net/yade/+question/233232 Christian Jakob posted a new comment: Just to be sure, we are talking about the same script. This part is working fine: F_buo = -1*(pi/3)*dh*dh*(3*rad - dh)*rho_f*integrator.gravity # = -V*rho*g #apply buoyancy force (will overwrite old forces applied with addF command) if b.isStandalone and isinstance(b.shape,Sphere): O.forces.addF(b.id,F_buo,permanent=True) if b.isClump: keys = O.bodies[b.id].shape.members.keys() for ii in range(0,len(keys)): O.forces.addF(keys[ii],(O.bodies[keys[ii]].state.mass/b.state.mass)*F_buo,permanent=True) But this affects only spheres (no clumps): F_buo = -1*(pi/3)*dh*dh*(3*rad - dh)*rho_f*integrator.gravity # = -V*rho*g #apply buoyancy force (will overwrite old forces applied with addF command) O.forces.addF(b.id,F_buo,permanent=True) I can imagine, that you used this (which explains differences between your and my results): F_buo = -1*(pi/3)*dh*dh*(3*rad - dh)*rho_f*integrator.gravity # = -V*rho*g #apply buoyancy force (will overwrite old forces applied with addF command) if b.isStandalone and isinstance(b.shape,Sphere): O.forces.addF(b.id,F_buo,permanent=True) #indents are lost, sorry! -- 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 #233232]: F
Question #233232 on Yade changed: https://answers.launchpad.net/yade/+question/233232 Christian Jakob posted a new comment: ... to be more precise: I also tested this 2-liner, which works equally well: O.bodies.appendClumped([sphere((0,0,0),1),sphere((1,0,0),1)]) O.forces.addF(2,(1,0,0),True) this works well, but in buoyancy example spheres and clumps gets separated. sphere are floating, clumps are falling ... -- 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 #236206]: Capillary force as a function of interparticle distance
Question #236206 on Yade changed: https://answers.launchpad.net/yade/+question/236206 Status: Open = Answered Christian Jakob proposed the following answer: Hi, Fixed via: - setting a label to capillary law: Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=40,createDistantMeniscii=True,label='Cap') - adding these lines at the end of the script (before O.run() command): O.run(1,True) Cap.createDistantMeniscii=False -- 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 #236323]: Yade-dem.org problem?
Question #236323 on Yade changed: https://answers.launchpad.net/yade/+question/236323 Christian Jakob posted a new comment: Confirm! -- 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 #236721]: Export periodic cell to vtk
Question #236721 on Yade changed: https://answers.launchpad.net/yade/+question/236721 Christian Jakob proposed the following answer: afaik, there is no such tool. i would definitely use this new feature! thanks in advance -- 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 #236206]: Capillary force as a function of interparticle distance
Question #236206 on Yade changed: https://answers.launchpad.net/yade/+question/236206 Christian Jakob posted a new comment: Hi, Sorry for this late reply. Are you still facing this problem? Christian. -- 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 #235985]: Highest clump position in clump cloud
Question #235985 on Yade changed: https://answers.launchpad.net/yade/+question/235985 Status: Open = Answered Christian Jakob proposed the following answer: hi sergio, it is a bit tricky, but possible. one idea: - give your clump clouds different color - loop over clumps - find maximum example: for b in O.bodies: if b.color= ???: maxHeight1 = max(maxHeight1,b.state.pos[2]) if b.color= ???: maxHeight2 = max(maxHeight2,b.state.pos[2]) ... regards, christian p.s. i know, that i should use launchpad, but i can not access it right now. blocked in china?! -- 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 #235552]: How to pass a python class to a C++ class
Question #235552 on Yade changed: https://answers.launchpad.net/yade/+question/235552 Christian Jakob posted a new comment: Hi, I tried to pass it as a python::object Yes this is tricky. An example how you can pass python lists and objects is here: https://github.com/yade/trunk/blob/master/py/wrapper/yadeWrapper.cpp#L201 Can you explain what your ClumpsFactory should do? Creating/modifying clumps is already implemented in several methods. Please also have a look at the manual section about clumps https://yade-dem.org/doc/user.html#clumping-particles-together regards, christian -- 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 #235259]: How to count number of interactions
Question #235259 on Yade changed: https://answers.launchpad.net/yade/+question/235259 Status: Open = Answered Christian Jakob proposed the following answer: Maybe the problem is, that the material id can not be used for specify a material (not sure, if this is the problem) instead of frictMatId = O.materials.append(NormalInelasticMat(density=300.0,poisson=0.5,young=1e9)) try this myMat = NormalInelasticMat(density=300.0,poisson=0.5,young=1e9) frictMatId = O.materials.append(myMat) Then you can set it to your sphere: ballId=O.bodies.append([sphere(center=(0.5,0.5,0.5),radius=.05,material=myMat)]) Getting number of interactions can be done like this: counter = 0 for i in O.interactions: if i.isReal: counter+=1 (did not check) -- 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 #235259]: How to count number of interactions
Question #235259 on Yade changed: https://answers.launchpad.net/yade/+question/235259 Christian Jakob proposed the following answer: My actual problem is, that only 1 interaction does appear in the container O.interactions. How many interactions would you expect from a simple model with one sphere and facets? But none of the collisions between the ball and the floor. Why is that? There are no interactions between facets, so your interaction must be between your sphere and one of the facets. You can check which bodie ids are involved in interactions by: for i in O.interactions: if i.isReal: print 'found interaction between ',i.id1,' and ',i.id2 -- 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 #235138]: how to define friction at the boundaries in triaxial test
Question #235138 on Yade changed: https://answers.launchpad.net/yade/+question/235138 Christian Jakob proposed the following answer: Hi, how ever, i have some additional development during the simulation. an error named FATAL /build/buildd/yade-0.70.0/core/ThreadRunner.cpp:31 run: Exception occured: Body #37 has velocity==NaN! occured and this automatically terminates the simulation. do you know anything that could be the cause of this problem and how to fix it?. 1.) follow statements/remarks of jduriez and do something like: idSand=O.materials.append(..,label='sandLabel')) Then you can do: sandLabel.young=1e6 or O.materials[idSand].young=1e6 which is doing the same ... But changing the property of a material will not change material properties of existing bodies. Do you want to change material properties during the run? 2.) If your simulation is dynamical, time step could be too high. Decrease it either directly with a fixed O.dt= ...? or decrease timestepSafetyCoefficient in GlobalStiffnessTimeStepper engine. Additionally I want to say, that GravityEngine is deprecated. You can set it directly in the NewtonIntegrator NewtonIntegrator(gravity=(0,0,-9.81),damping=0.7,label='integrator') and remove this line: GravityEngine(gravity=(0,0,-9.81)), Next time please send a short example script, that shows the problem. It is quite time consuming investigating a 1000-line script including 950 irrelevant lines. Regards, Christian -- 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 #235138]: how to define friction at the boundaries in triaxial test
Question #235138 on Yade changed: https://answers.launchpad.net/yade/+question/235138 Status: Open = Answered Christian Jakob proposed the following answer: Hi stephen, Welcome to yade community. In yade you have to define materials (e.g. one material for particles and one for boundaries): ### properties: density = 2650 shear_modulus = 2e5 poisson_ratio = 0.15 young_modulus = 2*shear_modulus*(1+poisson_ratio) ### define materials: SphereMat = FrictMat(young=young_modulus,poisson=poisson_ratio,density=density,frictionAngle=0.5) BoundaryMat = FrictMat(young=young_modulus,poisson=poisson_ratio,frictionAngle=1) O.materials.append(SphereMat) O.materials.append(BoundaryMat) Then you can use the materials in the generation process: O.bodies.append(sphere([0,0,0],radius=R, material = SphereMat, fixed=True)) O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, material = BoundaryMat, dynamic=1)) Hope it helps, Christian -- 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 #235138]: how to define friction at the boundaries in triaxial test
Question #235138 on Yade changed: https://answers.launchpad.net/yade/+question/235138 Status: Open = Answered Christian Jakob proposed the following answer: in your script i see some inconsistencies, that can lead to unexpected results: # 1. problem: idSand=O.materials.append(FrictMat(young=young,poisson=pr,frictionAngle=radians(fa),label='idSand')) #frictionAngle=fa !!! # 2. problem: O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, material = BoundaryMat, dynamic=1)) #this is ok O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, dynamic=1)) #missing material = BoundaryMat !!! O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=16,dynamic=0)) #missing material = BoundaryMat !!! O.bodies.append(utils.geom.facetBox((0.035,0.035,0.105),(0.035,0.035,0.035),wallMask=15)) #missing material = BoundaryMat # 3.problem: sp=pack.SpherePack() sp.makeCloud(mincorner,maxcorner,psdSizes=psdSizes,num=number,psdCumm=psdCumm,distributeMass=1) sp.toSimulation() you did not set material for spheres. but i do not know how to do this with sp.toSimulation() i prefer to use: O.bodies.append([sphere(c,r,material=idSand) for c,r in sp]) #instead of sp.toSimulation() regards, christian -- 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 #233232]: F
Question #233232 on Yade changed: https://answers.launchpad.net/yade/+question/233232 Christian Jakob posted a new comment: Actually this force will be erased at the next step, right? right, i tested that already. It gives me an idea on a possible variant that would store user defined forces and apply them forever. It would be a usefull replacement of the painful ForceEngine and it doesn't sound difficult to implement. this sounds good for me too. if this new feature is available clumps/apply-buoyancy-clumps.py example could be rewritten ... -- 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 #233241]: how to get the body numbers who is contact with a sphere
Question #233241 on Yade changed: https://answers.launchpad.net/yade/+question/233241 Christian Jakob posted a new comment: did you try this? count = 0 for i in O.interactions.withBody(i): if i.isReal: count += 1 print count -- 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 #233241]: how to get the body numbers who is contact with a sphere
Question #233241 on Yade changed: https://answers.launchpad.net/yade/+question/233241 Status: Open = Answered Christian Jakob proposed the following answer: hm, personally I prefer using proofed and implemented methods: #if you have clumps in the model: coord_num = avgNumInteractions(considerClumps=True) #else: coord_num = avgNumInteractions(considerClumps=False) -- 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 #233241]: how to get the body numbers who is contact with a sphere
Question #233241 on Yade changed: https://answers.launchpad.net/yade/+question/233241 Christian Jakob proposed the following answer: Yes, a simple way to get the coordination numbers in whole, but is there something wrong with my script to get coordination numbers in a small regin. So write your tool, and check out if it gives same result as avgNumInteractions(). If yes, it works fine ;) Does avgNumInteractions(considerClumps=True) consider the contacts from the wall? Yes, once for sphere-wall contact, twice for sphere-sphere contacts. There was a similar discussion about coordination number, see here: https://answers.launchpad.net/yade/+question/207350 cheers, christian -- 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 #232941]: Capillary pressure between two particles
Question #232941 on Yade changed: https://answers.launchpad.net/yade/+question/232941 Christian Jakob proposed the following answer: Sorry for my totally dump suggestion in the last post. I checked out your script today and recognized, that you tried to plot a vector b0displ over iteration number. Small workaround to fix that: def addPlotData(): d0 = O.bodies[0].state.displ() d1 = O.bodies[1].state.displ() plot.addData(i=O.iter,b0displ=d0[2],b1displ=d1[2]) But you will get zero all the time, because meniscii is not created. So this you can do by changing these lines: r = 1e-4 #particle radius h = 1e-6 #praticle distance #create two sphere paticles# O.bodies.append([ utils.sphere(center=(0,0,0),radius=r,fixed=False), utils.sphere((0,0,2*r-h),r) ]) hih, c -- 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 #232941]: Capillary pressure between two particles
Question #232941 on Yade changed: https://answers.launchpad.net/yade/+question/232941 Status: Open = Answered Christian Jakob proposed the following answer: None of them don't affect the movement of particles. Did you play with suction parameter in the script I send you? It definitely changes the forces and rupture distance, just check it out. In addition, I have also changed the particle distant and it did not work. In your script aabbs are exactly touching, if particles are touching. So you need to enlarge aabb in this case: In order to allow capillary forces between distant spheres, it is necessary to enlarge the bounding boxes using Bo1_Sphere_Aabb::aabbEnlargeFactor and make the Ig2 define define distant interactions via interactionDetectionFactor. It is also necessary to disable interactions removal by the constitutive law (Law2). - https://yade-dem.org/doc/yade.wrapper.html?highlight=capillary#yade.wrapper.Law2_ScGeom_CapillaryPhys_Capillarity Cheers, christian -- 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 #232941]: Capillary pressure between two particles
Question #232941 on Yade changed: https://answers.launchpad.net/yade/+question/232941 Christian Jakob posted a new comment: ... , where aabb = axis aligned bounding box: https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Aabb -- 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 #232941]: Capillary pressure between two particles
Question #232941 on Yade changed: https://answers.launchpad.net/yade/+question/232941 Christian Jakob posted a new comment: Sorry if my suggestions are wrong, I can not check out your script right now. But I think it is just a problem in order of commands. If you plot before calculation, you will never get data in your plot. So this part... plot.plots={'i':('b0displ')} plot.plot() plot.live=True plot.autozoom=True O.dt=0.5*PWaveTimeStep() qt.View() ...should be at the end of your script (after O.run command) -- 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