yade-users team mailing list archive

[Christian Jakob <question175394@xxxxxxxxxxxxxxxxxxxxx>]

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

Christian Jakob posted a new comment:
I am using hertz model with following properties:

poisson_ratio			= 0.3
young_modulus			= 2*shear_modulus*(1+poisson_ratio)
friction_coeff			= 2
angle					= math.atan(friction_coeff)
rho_p					= 2650

My materials and engines:

id_FacetMat=O.materials.append(FrictMat(young=young_modulus,poisson=poisson_ratio,frictionAngle=angle))
id_SphereMat=O.materials.append(FrictMat(young=young_modulus,poisson=poisson_ratio,density=rho_p,frictionAngle=angle))

FacetMat=O.materials[id_FacetMat]
SphereMat=O.materials[id_SphereMat]

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_CapillaryPhys(label='ContactModel')],
		[Law2_ScGeom_FrictPhys_CundallStrack()],
		),
	Law2_ScGeom_CapillaryPhys_Capillarity(CapillaryPressure=10000),
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=5,timestepSafetyCoefficient=0.8, defaultDt=4*utils.PWaveTimeStep()),
	GravityEngine(gravity=(0,0,-9.81)),
	NewtonIntegrator(damping=local_damping),
]
ContactModel.betan=viscous_normal # normal direction
ContactModel.betas=viscous_shear # shear direction
ContactModel.useDamping=True

I think this is ok, or did I do a mistake?

In PFC I use the same properties shear mod.= 3e7, poisson = 0.3 and rho = 2650.
The only difference is the wall, where I use linear model with stiffnesses kn=1e10 and ks=0 in PFC.

Can I do the same in Yade? ... Using hertz model for particles and
linear model for the facets?

christian

p.s. with timeStepUpdateInterval=5 calculation is stable ... PFC uses 1
as default update interval, I did not try to set it to 100 ...

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.