yade-dev team mailing list archive

[Janek Kozicki <janek_listy@xxxxx>]

noreply@xxxxxxxxxxxxx said:     (by the date of Tue, 06 Jul 2010 08:17:18 -0000)

>  add spheresVolume python wrapper for TriaxialStressController. (I need it to calculate void_ratio=voids/spheresVolume)

hmm, Bruno, actually this spheresVolume wrapper is returning a wrong
value. Forutantely I am able to calculate void ratio without it:

def void_ratio():
	p=triax.porosity  # p=Vo/V
	V=triax.boxVolume
	Vo=p*V  # voids volume
	Vs=V-Vo # spheres volume is total volume without voids
	#Vs=triax.spheresVolume
	#V=triax.boxVolume
	#Vo=V-Vs
	return Vo/Vs

But I thought that using spheresVolume would be nicer.

The are three solutions:

1) We need to either remove that python wrapper for spheresVolume (I do it)

2) Or we need to fix the spheresVolume (TriaxialStressController.cpp
   line 210) and you do it.

3) add a function that actually returns voidRatio.

Which solution do we choose?
-- 
Janek Kozicki                               http://janek.kozicki.pl/  |