Scene radius must be positive - Ignoring value
Hello:
I'm trying to generate a cylinder sample with target porosity, the MWE comes from [1], after I run the MWE, I click "show 3D" then click "Reference" or "Center" to show the model in different direction, there will be:
In [1]: Scene radius must be positive - Ignoring value
With this notification, the scene moves strange when I adjust it through mouse.
My yade version is 2018.02b with Ubuntu 18.04, here is the MWE:
################
from yade import pack
num_spheres=1000
key='_triax_base_'
targetPorosity = 0.45
compFricDegree = 30
finalFricDegree = 30
rate=-0.02
damp=0.2
stabilityThresh
young=5e6
mn,mx=Vector3(
O.materials.
O.materials.
walls=aabbWalls
wallIds=
sp=pack.
sp.makeCloud(
O.bodies.
Gl1_Sphere.
triax=TriaxialS
maxMultiplier=
finalMaxMultip
thickness = 0,
stressMask = 7,
internalCompac
)
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
TriaxialStateR
newton
]
triax.goal1=
while 1:
O.run(1000, True)
unb=unbalance
print 'unbalanced force:',unb,' mean stress: ',triax.meanStress
if unb<stabilityTh
break
print "### Isotropic state saved ###"
import sys
while triax.porosity>
compFricDegree = 0.95*compFricDegree
setContactFric
print "\r compFrictionDegree: ",compFricDegree," porosity:
sys.stdout.flush()
O.run(500,1)
print "### Compacted state saved ###"
#
aabb=utils.
predX=(
predY=(
predR=min(
pred=pack.
# ### erase walls and spheres which is not in cylinder
def cylinderShape():
for i in O.bodies:
if not isinstance(
O.bodies.
if isinstance(
if not pred(i.state.pos, i.shape.radius):
O.
cylinderShape()
###############
Could you please help me with this problem?
Thanks in advance!
Leonard
[1]https:/
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: