yade-users team mailing list archive

[Jan Stránský <question677865@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,

2)
###
bottom = [b for b in O.bodies if b.state.pos[2] < 0.1] # for example
for b in bottom: b.state.blockedDOFs = 'xyzXYZ'
###

1) it is possible to prescribe directly displacement by changing position:
###
for b in top: b.state.pos = b.state.pos + desiredDisplacement
###
but it is not recommended in YADE. Better prescribe velocity (together with some DOFs blocked):
###
for b in top:
   b.state.blockedDOFs = 'z' # e.g. only z displacement blocked
   b.state.vel = desiredDisplacement / O.dt
O.step()
for b in top:
   b.state.vel = Vector3.Zero
###

or instead of sudden jump, you can prescribe velocity for linear
increase of displacement.

cheers
Jan

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