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 :)
>>>
>>
>