yade-users team mailing list archive

[ceguo <question187798@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

ceguo is still having a problem:
Hi Jan,

Actually not. The difference is NOT format difference. Here is an
example (if you remember the script I sent to you before):

>>> packing.consolidate()   # calls Peri3dController and labels as 'cons'
>>> packing.getStress()      # this calls utils.stressTensorOfPeriodicCell()
Matrix3(-48965.9,100.619,0, 100.619,-48972.5,0, 0,0,0)
>>> cons.stress               # labelled as 'cons'
Vector6(-50002.634879037054,-50002.541124107054,0, 0,0,100.14402045344963)
>>> packing.shear()            # calls Peri3dController and labels as 'biax'
>>> packing.getStress()      # this calls utils.stressTensorOfPeriodicCell()
Matrix3(-38680.7,300.731,0, 300.731,-56123,0, 0,0,0)
>>> biax.stress                # labelled as 'biax'
Vector6(-39449.714550942866,-57438.815800571429,0, 0,0,312.58387791336929)

You see, the results are different. May this be due to 2d simulation? To use the current formula in 2d, I set:
O.cell.hSize=utils.Matrix3(3,0,0, 0,3,0, 0,0,1). Here the size in the third direction is set to unity. I think this will work smoothly without modification within the code. Am I right?

Ning

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.