Re: [Yade-users] [Question #707069]: How to determine if a cohesive link is broken (Law2_ScGeom6D_CohFrictPhys_CohesionMoment)
Question #707069 on Yade changed: https://answers.launchpad.net/yade/+question/707069 Status: Open => Answered Robert Caulk proposed the following answer: Hello, Have you tried: ``` for i in O.interactions: if hasattr(i.phys, 'cohesionBroken'): if i.phys.cohesionBroken: print(f"Interaction {i} is broken") ``` ? -rob -- 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 #707069]: How to determine if a cohesive link is broken (Law2_ScGeom6D_CohFrictPhys_CohesionMoment)
New question #707069 on Yade: https://answers.launchpad.net/yade/+question/707069 Hi, I use cohesive material (CohFrictMat) and the contact law of Law2_ScGeom6D_CohFrictPhys_CohesionMoment to introduce cohesive bonds in my DEM model. In the description of Law2_ScGeom6D_CohFrictPhys_CohesionMoment [1], it says that if the maximum tensile or maximum shear force is reached and CohFrictPhys::fragile =True (default), the cohesive link is broken. I would like to ask how to know if a cohesive link is broken or not. Thanks Leonard [1]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Law2_ScGeom6D_CohFrictPhys_CohesionMoment -- 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 #706821]: Particle breakage of clump using JCFpm
Question #706821 on Yade changed: https://answers.launchpad.net/yade/+question/706821 etthan posted a new comment: 1) I'd say yes. You could even consider a JCFPM mat for the walls assigning a different type to the sphere to have purely frictional interactions between spheres and walls. For instance, walls' type=0 and spheres' type=1 will define a purely frictional interaction between spheres and walls with a friction angle equal to the minimum value you assigned to both (if walls' friction angle = 0, then you will have frictionless interactions between spheres and walls). The same logic can be applied to the clumps if you wanted to have purely frictional interactions between clumps (e.g., define type=n for each clump). 2) You can assign any loading to the walls but be careful not to apply both displacement and force at the same time to avoid conflict in the control. Also, I think you need to make them nonDynamic to avoid their positions to be adjusted according to Newton's second law of motion. The quasi-staticity is not really related to the way you apply the load but more to the amplitude of the load vs. the system response. Not sure I get your point here... Maybe you refer to oscillations you observe? 3) You'd need to define some custom functions. The JCFPM can track the amount of cracks for each particle (and the degree of damage by relating the amount of broken contacts to the initial number of bonds). If you know which particle is within each clump, I guess you could assess their degree of damage too. For that, you'd need to list your clumps and assign them with the spheres ids. Sorry for not getting into the details of your MWE but it is not exactly Minimal. I hope I could help ANW. -- 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 #707064]: Influence of material
New question #707064 on Yade: https://answers.launchpad.net/yade/+question/707064 I was conducting a triaxial test on the argillite, which contained a matrix and two inclusions, to investigate its small strain stiffness. I tried to use spherical grains and JCF model to simulate this kind of rock. And I seperated the particle into three groups and assigned them different material properties (Young's modulous and Poisson's ratio). But the results were not desired. Because the small strain stiffness of the sample is almost equivalent to the sample containg only one material, which is confusing. Theoretically, the change of materials in the sample will affect the properties of the rock. How can I realize this goal? Could you give me some advice? from yade import ymport, utils , plot O.load('isocompact200kpa2.yade.bz2') O.engines[6].dead=True # ENGINES ARE DEFINED HERE try: os.mkdir('post') except: pass PACKING='121_1k.spheres' OUT=PACKING+'_1MPa_r0.02' iterMax=5 # maximum number of iterations saveVTK=iterMax/5 # Vtk files record interval YOUNG2=100e9 YOUNG3=50e9 ALPHA2=0.06 ALPHA3=0.27 bodies2=[b for b in O.bodies if 5.0 0.001 and e22 < 0.05: #O.save('save'+'{:.4f}'.format(e22)+'yade.bz2') #ini_e22a = e22 #if (e22 - ini_e22b) > 0.5: print ('unbalanced force: %f, mean stress: %f, s11: %f, s22: %f, s33: %f, coordination number: %f, porosity: %f' %(unb,mStress,-triax.stress(triax.wall_right_id)[0],-triax.stress(triax.wall_top_id)[1],-triax.stress(triax.wall_front_id)[2],avgNumInteractions(),Porosity)) #Output() #ini_e22b = e22 plot.saveDataTxt('loadinglog3.txt.bz2') plot.plots={'e22':('s22',None,'ev')} plot.plot() #O.run(iterMax) -- 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