yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #14520
Re: [Question #593653]: How to simulate a sinusoidal-shaped pressure condition at the top boundary of the packing
Question #593653 on Yade changed:
https://answers.launchpad.net/yade/+question/593653
bin deng gave more information on the question:
Hi Robert,
I am sorry that there are more lines in the FLOW section.
## ______________ 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
from yade import plot
def history():
plot.addData(p=flow.getPorePressure((1,0.8,0.35)),n=flow.porosity,t=O.time)
O.engines=O.engines+[PyRunner(iterPeriod=200,command='history()',label='recorder')]
O.engines=O.engines+[PyRunner(iterPeriod=200,command='flow.saveVtk()')]
from yade import plot
plot.plots={'t':('p',),'t ':('n',)}
plot.plot()
O.saveTmp()
O.timingEnabled=1
from yade import timing
print "starting oedometer simulation"
O.run(200,1)
timing.stats()
*************************************************************************
--
You received this question notification because your team yade-users is
an answer contact for Yade.