yade-users team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

Script is in attachment.
______________________________
[ENG] Best Regards
[GER] Mit freundlichen Grüßen
[RUS] С наилучшими пожеланиями
[UKR] З найкращими побажаннями

Anton Gladkyy


2009/6/25 Václav Šmilauer <eudoxos@xxxxxxxx>

>
>  If I wait a little bit more I get  messages included to attach file.
>>
> Backtrace clearly shows crash in the collider (the <signal handler called>
> line). Please post your whole script (engines etc).
>
> V
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users>
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users>
> More help   : https://help.launchpad.net/ListHelp
>

from yade import pack
from yade import utils
from euclid import *

#_______________________________________________________________________
PI = 3.14159265;
  
a = 62 #Size a, [mm] X
b = 40 #Size b, [mm] Y
h = 33 #Size h, [mm] Z

grainSize = 1 #Grain size, [mm]
distBetweenGrains = 1 #Distance between grains, [mm]

rho=1000 #Grain density, [kg/m^3]
poisson = 0.3 #Poisson coefficient
young = 30e9 #Young module

boxThick = 3 #Box thickness, [mm]
boxCoef = 5 #In how many times box larger, than an example
#_______________________________________________________________________
a = 0.001*a
b = 0.001*b
h = 0.001*h
grainSize = 0.001*grainSize
distBetweenGrains = 0.001*distBetweenGrains
boxThick = 0.001*boxThick

#_______________________________________________________________________



kw={'density':rho,'young':young,'poisson':poisson}

O.bodies.append(utils.alignedFacetBox((0,0,0),(a/2*boxCoef,b/2*boxCoef,boxThick*boxCoef),31,color=[0.98,0.984,0.576]))

O.bodies.append(pack.regularHexa (pack.inAlignedBox((-a/2,-b/2,boxThick/2),(+a/2,+b/2,h+boxThick/2)),radius=grainSize/2,gap=distBetweenGrains,color=(0.282,0.875,1),**kw))

#_______________________________________________________________________

try:
	from yade import qt
	qt.Controller()
except ImportError: pass

O.engines=[
	BexResetter(),
	BoundingVolumeMetaEngine([InteractingSphere2AABB(),InteractingFacet2AABB(),MetaInteractingGeometry2AABB()]),
	InsertionSortCollider(),
	InteractionDispatchers(
		[ef2_Sphere_Sphere_Dem3DofGeom(),ef2_Facet_Sphere_Dem3DofGeom()],
		[SimpleElasticRelationships()],
		[ef2_Dem3Dof_Elastic_ElasticLaw()],
	),
  GravityEngine(gravity=[0,0,-9.81]),
	NewtonsDampedLaw(damping=.2)
]

O.dt=2*utils.PWaveTimeStep()
O.saveTmp('init')
#for b in O.bodies: