yade-users team mailing list archive

[Diego Abarca <question693568@xxxxxxxxxxxxxxxxxxxxx>]

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

Diego Abarca gave more information on the question:
Thanks Jan, I modified the showValues function so it includes your
comments:

def showValues():    
    # This function is just for the MWE
    print("Current iteration = ", O.iter,
          ", s11 = ", utils.getStress()[0,0],  
          ", s22 = ", utils.getStress()[1,1],
          ", s33 = ", utils.getStress()[2,2],
          ", dvol = ", O.cell.volume-vol0,
          ", coorNum = ", utils.avgNumInteractions(skipFree=False))

And these are the results I get:

Start:  2020-10-22 10:16:24.653917
Cloud created, number of particles =  480
Compacting to target porosity:  2020-10-22 10:16:24.674490
Target porosity passed, porosity =  0.49986744040494263
Starting isocompaction simulation:  2020-10-22 10:16:47.926089
Starting shear simulation:  2020-10-22 10:19:05.654546
Current iteration =  1563876 , s11 =  -99919.50024856672 , s22 =  -100201.48450320616 , s33 =  -100092.81615722248 , dvol =  -6.617444900424222e-24 , coorNum =  2.5125
Current iteration =  1863876 , s11 =  0.0007813920667472948 , s22 =  -0.018439345804731364 , s33 =  -0.0008659837088851886 , dvol =  -2.0696058926076752e-19 , coorNum =  0.004166666666666667
Current iteration =  2163876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -3.751303782597383e-18 , coorNum =  0.0
Current iteration =  2463876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -8.511132637799019e-18 , coorNum =  0.0
Current iteration =  2763876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -1.3271298982690577e-17 , coorNum =  0.0
Current iteration =  3063876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -1.7382916022671158e-17 , coorNum =  0.0
Current iteration =  3363876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -2.150842307949773e-17 , coorNum =  0.0
Current iteration =  3663876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -2.686728951686537e-17 , coorNum =  0.0
Current iteration =  3963876 , s11 =  0.0 , s22 =  0.0 , s33 =  0.0 , dvol =  -3.423459680362606e-17 , coorNum =  0.0
Finished simulation phase:  2020-10-22 10:22:12.170592

So indeed, the number of interactions drops to zero at the same time that stresses drop to zero, while volume stays "constant".  
dvol goes from around e-24 to e-19, and I would think that is realistically like going from zero to zero (volume is of the order of 1e-8), so I still don't understand why the particles in my simulation have the space to loosen up so dramatically. Is there another way to measure volume of the cell that would show that the cell is actually growing?

Many thanks.

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