yade-users team mailing list archive

[liucheng83 <question247705@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,
I want to try simulate some particles fall down in fluid under gravity force. But it get the wrong message the same as the question #241071  posted by Ehsan Izadi , can anybody can help me?  https://answers.launchpad.net/yade/+question/241071
The script is here:

# basic simulation showing sphere falling ball gravity,
# bouncing against another sphere representing the support
from yade import pack
# DATA COMPONENTS
num_spheres=1000# number of spheres
young=1e6
compFricDegree = 3 # initial contact friction during the confining phase
finalFricDegree = 30 # contact friction during the deviatoric loading
mn,mx=Vector3(0,0,0),Vector3(1,1,1) # corners of the initial packing

O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)

sp=pack.SpherePack()
sp.makeCloud(mn,mx,-1,0.3333,num_spheres,False, 0.95,seed=1) #"seed" make the "random" generation always the same
sp.toSimulation(material='spheres')

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

# simulation loop -- see presentation for the explanation
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()],label="iloop"
	),
	FlowEngine(dead=1,label="flow"),#introduced as a dead engine for the moment, see 2nd section
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
	triax,
	NewtonIntegrator(damping=0.1,gravity=(0,0,-1*9.8)),
]


triax.wall_bottom_activated=False
triax.wall_top_activated=False
triax.wall_left_activated=False
triax.wall_right_activated=False
triax.wall_back_activated=False
triax.wall_front_activated=False



# set timestep to a fraction of the critical timestep
# the fraction is very small, so that the simulation is not too fast
# and the motion can be observed
O.dt=.5e-4*utils.PWaveTimeStep()

O.run(1,1)
print "###------  ###"

# save the simulation, so that it can be reloaded later, for experimentation
O.saveTmp()

###---------------------------------------------------------------------------------------------------add by lc

####A. Check bulk modulus of the dry material from load/unload cycles
###triax.stressMask=2
###triax.goal1=triax.goal3=0

###triax.internalCompaction=False
###triax.wall_bottom_activated=False
####load
###triax.goal2=11000; O.run(2000,1)
####unload
###triax.goal2=10000; O.run(2000,1)
####load
###triax.goal2=11000; O.run(2000,1)
###e22=triax.strain[1]
####unload
###triax.goal2=10000; O.run(2000,1)

###e22=e22-triax.strain[1]
###modulus = 1000./abs(e22)

#B. Activate flow engine and set boundary conditions in order to get permeability
flow.dead=0
flow.defTolerance=0.3
flow.meshUpdateInterval=200
flow.useSolver=3
flow.permeabilityFactor=1
flow.viscosity=1
flow.bndCondIsPressure=[0,0,0,0,0,0]
flow.bndCondValue=[0,0,0,0,0,0]
flow.boundaryUseMaxMin=[0,0,0,0,0,0]
O.dt=0.1e-3
O.dynDt=False

O.run(1,1)


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