← Back to team overview

yade-users team mailing list archive

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.