← Back to team overview

yade-users team mailing list archive

[Question #289705]: Minimal cohesive simulation

 

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

Hello. Can I modify the following script (modified from https://answers.dogfood.paddev.net/yade/+question/254316), in order to the upper sphere drags the lower one, by means of a cohesive contact? Thanks a lot in advance

from yade import *
import sys, time, random, os, gts, math
from yade import ymport
from yade import pack
from yade import qt
from yade import bodiesHandling

O.materials.append(CohFrictMat(
	young=3e8,
	poisson=0.3,
	frictionAngle=radians(30),
	density=2600,
	isCohesive=True,
	alphaKr=0.031,
	alphaKtw=0.031,
	momentRotationLaw=False,
	etaRoll=5.0,
	label='granular_material',
    normalCohesion = 10.**9.,
    fragile = False))

s = utils.sphere(center=[0,0,0],radius=1)
p = utils.sphere(center=[0,0,2],radius=1)

vel = 0.0001

O.bodies.append(s)
O.bodies.append(p)

O.bodies[1].state.vel=(0,0,vel)

O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb()]),
    InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom6D()],
        [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
        [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]),
    NewtonIntegrator(damping=0.2)
]

O.dt = utils.PWaveTimeStep()
O.step()

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