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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
