yade-users team mailing list archive

[Ehsan Izadi <question240920@xxxxxxxxxxxxxxxxxxxxx>]

Question #240920 on Yade changed:
https://answers.launchpad.net/yade/+question/240920

    Status: Answered => Open

Ehsan Izadi is still having a problem:
1) This is the materials section : 
-------------------
O.materials.append(FrictMat(young=1e10, poisson=0.3, frictionAngle=radians(30), density=2600,label='particles'))
---------------------

2) Right after I type "O.step()" and enter, all particles will be
vanished. (Normally in the 3D viewport if you zoom into a particular
part of your model and then press shit+c it will be zoomed in such a way
you will be able to see whole scene including all existing bodies)

3) I also changed "defaultDt" and reduced it until 1e-17, still the same
error happen.

4) Reduction of the timeStepSafetyCoefficient did not work; I tried 0.5,
0.2, 0.1, 0.001 values for that.

Note: my code is not so long. actually, I am at the first step of modelling and my code is not very much developed, so, it is not so long. Here is the full script:
---------------------------------------------
#import ymport & pack libraries from YADE
from yade import ymport, pack

#defining material properties

O.materials.append(FrictMat(young=1e10, poisson=0.3, frictionAngle=radians(30), density=2600,label='particles'))
O.materials.append(FrictMat(young=1e10, poisson=0.3, frictionAngle=0, density=0, label='boxWalls'))

#importing facet boundaries to viewport
box=ymport.stl('/home/ehsan/Desktop/try3.stl')
O.bodies.append(box)

'''
O.bodies.append(utils.facetBox((0,0,0.05),(0.2,0.15,0.05),wallMask=31))
'''

#spherepacking


sp=pack.SpherePack()
sp.makeCloud((-0.2,-0.15,0),(0.2,0.15,0.1),rRelFuzz=0.5, num=20000)
sp.toSimulation(material='particles')


triax=TriaxialStressController(
	maxMultiplier=1.+2e4/1e6, # spheres growing factor (fast growth)
	finalMaxMultiplier=1.+2e3/1e6, # spheres growing factor (slow growth)
	thickness = 0,
	stressMask = 7,
	max_vel = 0.005,
	internalCompaction=True # If true the confining pressure is generated by growing particles
)



#engine setup
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	NewtonIntegrator(gravity=(0,0,-9.81),damping=0.2),
	triax,
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8)
	]


triax.goal1=triax.goal3=triax.goal2=1000

O.saveTmp()
--------------------------------------------------------------------

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