yade-users team mailing list archive

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

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

    Status: Open => Answered

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

> if so, maybe i just need more time to see the difference? i didnt
print anything so the only way i can be sure that my code works is to
check 3D view by my eyes…i thought it should be visible as my specimen
is actual incompact.

depneding on "to see" is a bad approach. The only way how to be sure is
NOT checking 3D view, the only way how to be sure are the numbers in the
model.

Thanks for the code, could you try to make it minimal and also not-so-
long-running? I have started it a few minutes ago, still waiting for the
initial configuration to finish..

O.forces.setPermF(f.id,preStress*a*n) is the right approach

> f.state.blockedDOFs = 'XYZz'

check if normal is not parallel to z, which is blocked

Otherwise, do
###
... your code

f = facets[0] # or any other valid index
print(f.state.pos)
O.step()
print(f.state.pos) # should be different from the first one
###

cheers
Jan

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