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` (not `git 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 :)
>>>>>
>>>>
>>>
>>
>>
>