yade-users team mailing list archive

[Jérôme Duriez <question702267@xxxxxxxxxxxxxxxxxxxxx>]

Question #702267 on Yade changed:
https://answers.launchpad.net/yade/+question/702267

Jérôme Duriez proposed the following answer:
Hi Mehdi,

Does the following script look like what you would like to do ?

##########
O.materials.append(FrictMat(label='frict'))
O.materials.append(CohFrictMat(label='cohfrict'))

for x in [0,1.9]: O.bodies.append(sphere((x,0,0),1,fixed=True,material='frict'))
for x in [0,1.9]: O.bodies.append(sphere((x,3,0),1,fixed=True,material='cohfrict'))

O.engines = [ForceResetter()
             ,InsertionSortCollider([Bo1_Sphere_Aabb()])
             ,InteractionLoop([Ig2_Sphere_Sphere_ScGeom()],
                              [Ip2_FrictMat_FrictMat_CapillaryPhys(),Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
                              [Law2_ScGeom_FrictPhys_CundallStrack()])
             ,Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure = 10)
             ,NewtonIntegrator()
             ]

O.step()
########

I did not test thoroughly but 0-1 bottom row interaction should be
capillary+frictional and 2-3 top row should be cohesive-frictional

Let us know

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