yade-users team mailing list archive

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

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Ning,

try tu use utils.stressTensorOfPeriodicCell(True), this is the only
difference in both cases, as in both approaches the same c++ function is
called (but with different parameter by default).

Jan


ceguo píše v St 15. 02. 2012 v 09:11 +0000:
> 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.