yade-users team mailing list archive

[xuanshenyu <question703767@xxxxxxxxxxxxxxxxxxxxx>]

New question #703767 on Yade:
https://answers.launchpad.net/yade/+question/703767

Hi,everyone!

I write a case of merely two particles, in which the CohFritMat and Ip2_CohFrictMat_CohFrictMat_CohFrictPhys were used.
Besides, I added a pair of opposing and periodic forces, which are line increase with O.iter, to the two particles, and the python file as below: 

############
from yade  import pack
from numpy import genfromtxt 
from yade  import utils
from yade  import export
from yade  import ymport
from yade  import wrapper
import csv              
import numpy as np 
import gts, os.path, locale
#locale.setlocale(locale.LC.ALL,'en_US.UTF-8')
rho=920
gravity=0
#mat
mat1 = O.materials.append(CohFrictMat(isCohesive=True, frictionAngle=radians(15.6), density=rho_ice, poisson=0.33, young=5.358e8, alphaKr=0.242, alphaKtw=0.1, normalCohesion=2e5, shearCohesion=2e5))
Mat1=O.materials[mat1]

#sphere
coordinates=numpy.genfromtxt("zuobiao.csv",delimiter=",")
particle=[sphere((x,y,z),radius=0.01,material=Mat1,color=(1,1,0)) for x,y,z in coordinates]
Particle=O.bodies.append(particle)  

#add engines
O.engines=[
      ForceResetter(),
      InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.1), Bo1_Facet_Aabb()]),
      InteractionLoop(
              [Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.1)],
              [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True),Ip2_FrictMat_FrictMat_FrictPhys()],
              [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
      ),
      NewtonIntegrator(gravity=(0, 0, -gravity), damping=0.2),
      PyRunner(command='addForce()',iterPeriod=1),
      VTKRecorder(iterPeriod=100,recorders=['all'],fileName='out-')
]

def addForce():
     a=math.sin(O.iter*0.01)
    O.forces.addF(0,Vector3(0,0,-100*a))
    O.forces.addF(1,Vector3(0,0,100*a))
#timestep
O.dt = .1* PWaveTimeStep() 
O.step()
O.run(5000,True)
O.saveTmp()
#############
At the first timestep of the simulation, the bond of the two particles are created, and the interaction is Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(). At a moment, the bond is broken, and the bond is nerver established again.
However, when the two particle  come into contact again, what kind of contact law yade uses for dealing with the contact again? Ip2_CohFrictMat_CohFrictMat_CohFrictPhys or Ip2_FrictMat_FrictMat_FrictPhys()?




-- 
You received this question notification because your team yade-users is
an answer contact for Yade.