New question #706787 on Yade: https://answers.launchpad.net/yade/+question/706787
In the name of God Hi every one! I encountered an strange problem. When i add some spheres bonded to each other using CohFrictMat in a triaxial box, exactly after addition of bonded particles, the wall material type changes from 'walls' to 'cube'. here is the MWE: from yade import pack ############################################ ### DEFINING VARIABLES AND MATERIALS ### ############################################ # The following 5 lines will be used later for batch execution ConPre=-100000 ConPre1=-5000 key='t_20%_25_100' # put you simulation's name here targetPorosity = 0.38#the porosity we want for the packing compFricDegree = 30 # initial contact friction during the confining phase (will be decreased during the REFD compaction process) finalFricDegree = 30 # contact friction during the deviatoric loading young=1e9 # contact stiffness mn,mx=Vector3(0,0,0),Vector3(.3,.3,.3) # corners of the initial packing O.materials.append(CohFrictMat(young=1e9,poisson=0.2,density=1500,isCohesive=False,frictionAngle=0,label='walls')) ## create walls around the packing walls=aabbWalls([mn,mx],thickness=0) wallIds=O.bodies.append(walls) ## create materials for spheres and plates O.materials.append(CohFrictMat(young=1e6,poisson=0.48,density=1500,label='cube',frictionAngle=radians(compFricDegree), isCohesive=True, normalCohesion=50000000, shearCohesion=50000000, momentRotationLaw=True )) sp=pack.SpherePack() sp.makeCloud((0,0,0),(.3,.3,.3),rMean=.05,distributeMass=False,seed=1) ############################# ### DEFINING ENGINES ### ############################ triax=TriaxialStressController( maxMultiplier=1.+2e4/young, # spheres growing factor (fast growth) finalMaxMultiplier=1.+2e3/young, # spheres growing factor (slow growth) thickness = 0, stressMask = 7, internalCompaction=False, # If true the confining pressure is generated by growing particles ) newton=NewtonIntegrator(damping=.1) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.12),Ig2_Box_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys( setCohesionNow=True, setCohesionOnNewContacts=False, )], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),Law2_ScGeom_FrictPhys_CundallStrack()] ), triax, newton, ] triax.goal1=triax.goal2=ConPre1 triax.goal3=ConPre1 O.dt=0 yade.qt.Controller() for x in [.1,.2]: for y in [.1,.2]: for z in [.1,.2]: s = utils.sphere((x,y,z),.08,material='cube') O.bodies.append(s) ######################################################## #######By adding the following lines to the terminal, the walls disppeared#### ######################################################## #O.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