Accessing Periodic Corner Volume Vectors for Porosity Calcs
Hello,
I would like to use the utils.voxelPorosity function to calculate the porosity during shear with the periTriaxContro
I've looked at how the it's called in the voxelPorosityTr
Any suggestions for the best way to do this?
I've tried to search for in-built methods to get the corner volume vectors. However, I do not see any relevant methods for the periTriaxController engine, nor do I see any in the attributes list:
$ triax.dict()
{'dynCell': True,
'goal': Vector3(
'stressMask': 3,
'maxStrainRate': Vector3(
'maxUnbalanced': 0.0001,
'absStressTol': 1000.0,
'relStressTol': 1e-05,
'growDamping': 0.02,
'globUpdate': 5,
'doneHook': 'startShear()',
'maxBodySpan': Vector3(
'stressTensor': Matrix3(
'stress': Vector3(
'strain': Vector3(
'strainRate': Vector3(
'stiff': Vector3(0,0,0),
'currUnbalanced': 2.9477521076913
'prevGrow': Vector3(
'mass': 15380.193739519347,
'externalWork': 413315.24427129014,
'dead': False,
'ompThreads': -1,
'label': 'triax'}
With `O.cell`, I seem to at least be able to get the dimensions with the `size` attribute, but I feel like this is incomplete without at least one corner. There is also an `hSize` attribute, but as far as I can tell, it seems O.cell.hSize is just O.cell.trsf * O.cell.refHSize.
I think otherwise, I would need to estimate the box based on the strains and initial dimensions, or approximate it based on the min/max limits of the body centers--but wanted to check if there is an out of box solution first.
Cheers,
Kevin
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Kevin Kuei
- Solved:
- Last query:
- Last reply: