FEniCS complete code for magneto-visco-mechanical homogenization: 3D full-field magneto-mechanical 
homogenization that accounts for: finite strains, viscoelasticity, incompressibility and mixed 
stress/magnetic flux control.
Copyright (C) 2021  Sergio Lucarini, Miguel Angel Moreno-Mateos, Kostas Danas, Daniel Garcia-Gonzalez

If using this code for research or industrial purposes, please cite:
S. Lucarini, M.A. Moreno-Mateos, K. Danas, D. Garcia-Gonzalez
Insights into the viscohyperelastic response of soft magnetorheological
elastomers: competition of macrostructural versus microstructural players
International Journal of Solids and Structures, 2022.

- Fenics_magneto_visco_homogenization.py
This python script contains the solver for the magneto-viscoelastic homogenization.
Input variables are the tests conditions, such as mechanical and magnetic loads and material properties.
Outputs of the solves are some macroscopic results and the .xdmf and .h5 with the full field results.

- aux_homog_nl.py
This is a python code that contains some functions used by the solver: periodic boundary conditions, 
input and output preparation.

- s30.xdmf and s30.h5
These files are an input geometry example that contains a sphere in the middle of un unitary cube.

- Usage instructions
Installation of Fenics module is needed.
The python modules numpy and matplotlib packages are required.
