yade-users team mailing list archive

[bin deng <question593653@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

bin deng gave more information on the question:
Hi Robert,

Thanks for you answer me so quickly.
I also get this erro when set waveAction=false. In fact, as long as I don't add the bndCondIsPressure in the code, it doesn't work.
Here is my code:

*************************************************************************
from yade import pack

num_spheres=5000# number of spheres
young=1.5e7
compFricDegree = 3 # initial contact friction during the confining phase
finalFricDegree = 30 # contact friction during the deviatoric loading
mn,mx=Vector3(0,0,0),Vector3(2,0.8,0.7) # 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)
pred=inAlignedBox((0,0,0),(2,0.8,0.7))

sp=pack.randomDensePack(pred,radius=0.0305,rRelFuzz=0.098,spheresInCell=num_spheres,returnSpherePack=True) 
sp.toSimulation(material='spheres')
newton=NewtonIntegrator(gravity=(0,-9.81,0))

from yade import qt

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(label="flow"),
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
        newton
]

## ______________   FLOW section   _________________
# Activate flow engine and set boundary conditions

flow.bndCondIsPressure=[0,0,0,0,1,0]
flow.waveAction=1
flow.sinMagnitude=500
flow.sineAverage=10
O.dt=0.1e-3
*************************************************************************
I know there should be some problems in my code.  :(
Can you please help me  modify my code or tell me how to modify it? I really want to simulate a sinusoidal-shaped pressure condition at the top boundary of the packing. Thanks for your reply again! : )

Best,
Deng

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.