TR1: Co-rotated and reduced order time integrators

Technical Report 1 describes a linearly implicit time integration method combined with a family of finite element formulations: (…)

  • Total Lagrangian (TL)
  • body co–rotational full mesh (BC)
  • body co–rotational reduced order (BC–RO)
  • body co–rotational modal (BC–MODAL)

as they are implemented in Solfec. These formulations allow to vary the amount of deformation expressed by a finite element mesh and facilitate saving computational time and storage. They may be helpful in the context of multibody modeling.

The methods described in TR1 can be accessed by creating a BODY object as follows

bdTL = BODY (solfec, 'FINITE_ELEMENT', shape, material, form = 'TL')
bdTL.scheme = 'DEF_LIM'

bdBC = BODY (solfec, 'FINITE_ELEMENT', shape, material, form = 'BC')
bdBC.scheme = 'DEF_LIM'

bdRO = BODY (solfec, 'FINITE_ELEMENT', shape, material,
             form = 'BC-RO', base = pod_base)

dbMD = BODY (solfec, 'FINITE_ELEMENT', shape, material,
             form = 'BC-MODAL', base = modal_base)

In case of the ‘BC–RO’ and ‘BC–MODAL’ the linearly implicit integrator is default. For ‘TL’ and ‘BC’ the default time integration scheme is ‘DEF_EXP’ (Table 7); we change it to ‘DEF_LIM’ after body creation.

The base format is a tuple of two lists

base = ([eval1, eval2, ..., evaln],
        [evec11, evec12, ..., evec1m,
         evec21, evec22, ..., evec2m,
         ..., evecnm])

where we assumed n base vectors of size m. The first list stores eigenvalues and the second list stores eigenvectors. COROTATED_DISPLACEMENTS and RIGID_DISPLACEMENTS commands can be used together to sample displacements and generate a reduced base using Python’s modred package. MODAL_ANALYSIS command can be used to generate a modal base.

Input decks for the rotating bar, pipe impact, and array excitation examples from TR1 can be respectively found in

directories. See also the online version of these examples.


comments powered by Disqus