yade-dev team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

I have done changes in 2362, but it still crashes.

Could you not try script in attachment. It crashes sometimes, not always.

Thanks
______________________________

Anton Gladkyy


2010/7/14 Václav Šmilauer <eudoxos@xxxxxxxx>

> > Vaclav, please, have a look at r2354, it can be very serious, as it
> > calls a  very spontain crash.
> OK, perhaps the issue is that I made States indexable and I missed
> RpmState as it is not exposed with YADE_PLUGIN -- please fix that first,
> then add createIndex() to RpmState ctor and REGISTER_CLASS_INDEX to
> class body (see CpmState how that is done exactly). I don't see why it
> would crash, but it would explain why it is only you who has the
> problem.
>
> v.
>
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
> More help   : https://help.launchpad.net/ListHelp
>

o=Omega()
intRadius = 1.0
frictionAngle = 0.25
young = 37e9
rho=2500 #Grain density, [kg/m^3]
poisson = 0.3 #Poisson coefficient
rad=2.0e-3

RockId=O.materials.append(RpmMat(young=young,frictionAngle=frictionAngle,poisson=poisson,density=rho,exampleNumber=10,initCohesive=True))

SpheresID=[]
id1=O.bodies.append(utils.sphere((0,0,0.0),radius=rad,material=RockId))
id2=O.bodies.append(utils.sphere((0,0,2.0*rad),radius=rad,material=RockId))
id3=O.bodies.append(utils.sphere((2.0*rad,0,0.0),radius=rad,material=RockId))
id4=O.bodies.append(utils.sphere((2.0*rad,0,2.0*rad),radius=rad,material=RockId))

id1=O.bodies.append(utils.sphere((0,0,0.0+5.0*rad),radius=rad,material=RockId))
#id2=O.bodies.append(utils.sphere((0,0,2.0*rad+5.0*rad),radius=rad,material=RockId))
id3=O.bodies.append(utils.sphere((2.0*rad,0,0.0+5.0*rad),radius=rad,material=RockId))
id4=O.bodies.append(utils.sphere((2.0*rad,0,2.0*rad+5.0*rad),radius=rad,material=RockId))

id1=O.bodies.append(utils.sphere((0,0,0.0+10.0*rad),radius=rad,material=RockId))
id2=O.bodies.append(utils.sphere((0,0,2.0*rad+10.0*rad),radius=rad,material=RockId))
id3=O.bodies.append(utils.sphere((2.0*rad,0,0.0+10.0*rad),radius=rad,material=RockId))
id4=O.bodies.append(utils.sphere((2.0*rad,0,2.0*rad+10.0*rad),radius=rad,material=RockId))

id2=O.bodies.append(utils.sphere((0,0,2.0*rad+15.0*rad),radius=rad,material=RockId))
id2=O.bodies.append(utils.sphere((0,0,2.0*rad+20.0*rad),radius=rad,material=RockId))

o.engines=[
	ForceResetter(),
	BoundDispatcher([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius,label='is2aabb'),Bo1_Facet_Aabb()]),
	InsertionSortCollider(),
	InteractionDispatchers(
		[Ig2_Sphere_Sphere_Dem3DofGeom(distFactor=intRadius,label='ss2d3dg'),Ig2_Facet_Sphere_Dem3DofGeom()],
		[Ip2_RpmMat_RpmMat_RpmPhys()],
		[Law2_Dem3DofGeom_RockPMPhys_Rpm()]
	),
	GravityEngine(gravity=[0,0,-9.810]),
	NewtonIntegrator(damping=.3),
	ParticleSizeDistrbutionRPMRecorder(iterPeriod=1,file='psd'),
]

o.dt=utils.PWaveTimeStep()*0.2
O.step()