Hello Robert!
Thank you very very much for your help :).
I will tell you something when I implemented the new active-fibers model :)
Congratulations again for such a nice and useful tool!!
Patricia
El miércoles, 4 de noviembre de 2015, 10:13:44 (UTC+1), Robert Cimrman escribió:
Hello Patricia,
I did some improvements, so now it should be much easier for you to implement the term. Again, it is in my development branch [1]. I force-pushed there, so when getting the code, use
git fetch
(notgit pull
!) and then reset your branch to my version.See also [2] (especially "synchronize your repository with the upstream" part).
r.
[1] https://github.com/rc/sfepy/tree/active-fibres-2 [2] http://sfepy.org/doc-devel/developer_guide.html#contributing-changes
On 11/02/2015 09:10 AM, Patricia Garcia Cañadilla wrote:
Hello Robert,
Thank you very much for your help :). I will look at it, but I can also wait a bit more :), it's not a problem!
El domingo, 1 de noviembre de 2015, 9:08:25 (UTC+1), Robert Cimrman escribió:
Hello Patricia,
I did some prototyping that might help you, see [1]. A minor tweak to HyperElasticBase class was needed so that the time stepper and fibre deformation data can be passed to the stress/tangent modulus evaluation functions, and hence to fibre_function().
So either work with that, or wait a bit more - I will make some proper commits next week (returning to my workplace).
r. [1] https://github.com/rc/sfepy/tree/active-fibres-2
On 10/28/2015 08:12 AM, Robert Cimrman wrote:
On 10/27/2015 12:02 PM, Patricia Garcia Cañadilla wrote:
But for calculating eps inside get_fargs function like in fibre_function:
eps = nm.zeros_like(mat1) omega = nm.empty_like(green_strain) for ii, (ir, ic) in enumerate(iter_sym(fdir.shape[2])): omega[..., ii, 0] = fdir[..., ir, 0] * fdir[..., ic, 0] eps[..., 0, 0] += omega[..., ii, 0] * green_strain[..., ii, 0]
I would need also as a input argument "green_strain"?
I guess not - you can pass both eps and eps0 (or deps) as arguments.
And then return both, fargs and deps not?:
"return deps, fargs"
I see the problem - I will check it soon(ish) :)
Sorry for the amount of silly questions, I am very new using Python
and
finite-element modeling :)
No problem. The above will need some tweaking of other code, to accommodate for user arguments.
Cheers, r.
El viernes, 23 de octubre de 2015, 10:39:07 (UTC+2), Patricia Garcia Cañadilla escribió:
Dear Robert,
I would like to the model the viscoelastic behavior of active
fibers, as
you done in the paper: "Modelling heart tissue using a composite muscle model with blood perfusion" or as here:
http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.488.5795&rep=rep1&type=pdf
and I was wondering if is easy to implemented in sfepy and/or if there is some predefined term in Sfepy, since for example the term "dw_tl_fib_a" doesn't include the "exp{k min(0;de/dt}}" that appears in the paper.
Thank you very much for your help :)