sfepy.terms.terms_fibres module

class sfepy.terms.terms_fibres.FibresActiveTLTerm(*args, **kwargs)[source]

Hyperelastic active fibres term. Effective stress S_{ij} = A f_{\rm max} \exp{\left\{-(\frac{\epsilon - \varepsilon_{\rm opt}}{s})^2\right\}} d_i d_j, where \epsilon = E_{ij} d_i d_j is the Green strain \ull{E} projected to the fibre direction \ul{d}.

Definition:

\int_{\Omega} S_{ij}(\ul{u}) \delta E_{ij}(\ul{u};\ul{v})

Call signature:
dw_tl_fib_a (material_1, material_2, material_3, material_4, material_5, virtual, state)
Arguments:
  • material_1 : f_{\rm max}
  • material_2 : \varepsilon_{\rm opt}
  • material_3 : s
  • material_4 : \ul{d}
  • material_5 : A
  • virtual : \ul{v}
  • state : \ul{u}
arg_shapes = {'virtual': ('D', 'state'), 'state': 'D', 'material_1': '1, 1', 'material_2': '1, 1', 'material_3': '1, 1', 'material_4': 'D, 1', 'material_5': '1, 1'}
arg_types = ('material_1', 'material_2', 'material_3', 'material_4', 'material_5', 'virtual', 'state')
family_data_names = ['green_strain']
get_eval_shape(mat1, mat2, mat3, mat4, mat5, virtual, state, mode=None, term_mode=None, diff_var=None, **kwargs)[source]
get_fargs(mat1, mat2, mat3, mat4, mat5, virtual, state, mode=None, term_mode=None, diff_var=None, **kwargs)[source]
name = 'dw_tl_fib_a'
static stress_function(out, pars, green_strain)[source]
static tan_mod_function(out, pars, green_strain)[source]
sfepy.terms.terms_fibres.fibre_function(out, pars, green_strain, fmode)[source]

Depending on fmode, compute fibre stress (0) or tangent modulus (!= 0).