yade-dev team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

There is an attached script.

Bruno, please, do not do direct uncommit. I have committed already something.
Thanks

Anton




On Wed, Nov 23, 2011 at 5:00 PM, Bruno Chareyre
<bruno.chareyre@xxxxxxxxxxx> wrote:
> How did you reproduce the bug Anton please?
>

#!/usr/bin/python
# -*- coding: utf-8 -*-

"""
To run this script you need to have all 10 text files from https://yade-dem.org/wiki/CapillaryTriaxialTest
in the same folder as this script!
"""

model_type		= 0	#1=Hertz model, else capillary model

shear_modulus 	= 1e5
poisson_ratio 	= 0.3
young_modulus	= 2*shear_modulus*(1+poisson_ratio)
friction		= 0.5
angle			= atan(friction)
local_damping 	= 0.01
lowercorner		= Vector3(0,0,0)
uppercorner		= Vector3(0.002,0.002,0.004)

id_SphereMat=O.materials.append(FrictMat(young=young_modulus,poisson=poisson_ratio,density=2500,frictionAngle=angle))
SphereMat=O.materials[id_SphereMat]

sp=pack.SpherePack()
sp.makeCloud(lowercorner,uppercorner,.0002,rRelFuzz=.3)
O.bodies.append([utils.sphere(c,r,material=SphereMat) for c,r in sp])

O.bodies.append(geom.facetBox(uppercorner/2,uppercorner/2,wire=True,fixed=True,material=SphereMat))

if model_type == 1:
	O.engines=[
		ForceResetter(),
		InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
		InteractionLoop(
			[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
			[Ip2_FrictMat_FrictMat_FrictPhys()],
			[Law2_ScGeom_FrictPhys_CundallStrack()],
		),
		GravityEngine(gravity=(0,0,-9.81)),
		NewtonIntegrator(damping=local_damping),
	]
else:
	O.engines=[
		ForceResetter(),
		InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
		InteractionLoop(
			[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
			[Ip2_FrictMat_FrictMat_CapillaryPhys()],
			[Law2_ScGeom_FrictPhys_CundallStrack()],
		),
		GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=5,timestepSafetyCoefficient=0.8,defaultDt=4*utils.PWaveTimeStep()),
		Law2_ScGeom_CapillaryPhys_Capillarity(CapillaryPressure=10000),	#pressure difference between air and fluid in [Pa]
		GravityEngine(gravity=(0,0,-9.81)),
		NewtonIntegrator(damping=local_damping),
	]

#O.dt=0.5*utils.PWaveTimeStep()
O.run(10000,True)

O.bodies.erase(0)

O.run(100)