yade-users team mailing list archive

[Christian Jakob <question165458@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I created a model with 1000 spheres and one facet. There I cannot fix two problems:

1.) The facet is in xy-plane, but I want to have it normal to (.02,.02,1), not (0,0,1)!
How can I do this?
I tried different things with facet, facetBox, wall and facetPolygon, but I could not fix it.

2.) Spheres are not moving and I dont know why. timing.stats() gives:

Name                                                    Count                 Time            Rel. time
-------------------------------------------------------------------------------------------------------
ForceResetter                                         0                  0us                           
InsertionSortCollider                                 0                  0us                           
InteractionLoop                                       0                  0us                           
GravityEngine                                         0                  0us                           
NewtonIntegrator                                      0                  0us                           
TOTAL                                                                    0us    



Now the script:

# script START ------------------------------------------------------------------------------------

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

# speed test file for YADE


num_balls1D = 10

num_balls = num_balls1D*num_balls1D*num_balls1D

logfilename = 'speed-test-output%i-balls.log' % (int(num_balls))

origin_wall = num_balls1D/2

x_pos = 0

y_pos = 0

z_pos = 0
friction=0.5
angle=atan(friction)

id_FacetMat=O.materials.append(ViscElMat(kn=1e8,ks=1e8,cn=0.0,cs=0.0,frictionAngle=angle))
id_SphereMat=O.materials.append(ViscElMat(kn=1e6,ks=1e6,cn=0.0,cs=0.0,density=1000,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_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
		[Ip2_ViscElMat_ViscElMat_ViscElPhys()],
		[Law2_ScGeom_ViscElPhys_Basic()],
	),
	GravityEngine(gravity=(0,0,-9.81)),
	NewtonIntegrator(damping=0.7),
]

id_facet=O.bodies.append(geom.facetBox((origin_wall,origin_wall,0),(100,100,0),material=FacetMat))

for ii in range(1,num_balls1D):

	x_pos = x_pos + 1

	y_pos = 0

	z_pos = 0

	for jj in range(1,num_balls1D):

		y_pos = y_pos + 1

		z_pos = 0

		for kk in range(1,num_balls1D):

			z_pos = z_pos + 1
			O.bodies.append(utils.sphere([x_pos,y_pos,z_pos], material=SphereMat, radius=0.5, wire=False, highlight=False))#, color=(50,50,0)


from yade import timing

O.dt=1e-3

start_time=O.realtime
print start_time

O.run(100,True)

stop_time=O.realtime
time_needed=stop_time-start_time
print time_needed

O.timingEnabled==True
timing.stats()

# script END   ------------------------------------------------------------------------------------

Christian.

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