Hi Kid,
you can use
hydraulic_conductivity = ev('ev_integrate_mat.3.Omega(m.k, p)', m=m,
mode='el_avg')...
out['hydraulic_conductivity'] = Struct(name='output_data', mode='cell',
data=hydraulic_conductivity)Essentially, ev_integrate_mat can be used to integrate any material parameter, check [1] to see what the mode argument means.
r.
[1] http://sfepy.org/doc-devel/src/sfepy/terms/terms_basic.html#sfepy.terms.term...
On 12/28/2016 09:01 AM, Kid Guo wrote:
Where Ki is the permeability of the i-th compartment and Kij is the general diffusion permeability. how to ouput the Kij of each cell.
在 2016年12月26日星期一 UTC+8下午8:52:44，Kid Guo写道：
Hi,
I have a problem based on "../sfepy/examples/multi_physics/biot.py". In the attached scripy 'block2.py', I have calculated the stress and strain according to the displacement, at the 111 and 112 lines you can see it.
Now I want to output the K of every cell, I gave a value of k at the 64 line in this script, it is a 3*3 array. I don't know which function or equation to use, so I was stucked at the 113 line. if I change the K value, will where be some difference results?
Thank you advance
K.