The talks of the third MFront Users Day (21 June 2017)

The talks of the third MFront Users Day are available here:

Cast3m 2017 is out (12 May 2017)

Cast3M 2017 has been released.

A binary version of TFEL compiled for Cast3M 2017 is now available for download on sourceforge:

Non linear constraints with MTest

Arbitrary non linear constraints can be now imposed in MTest using @NonLinearConstraint keyword.


Abritray non linear constraints can be used to:

On the physical meaning of a constraint

A constraint \(c\) is imposed by introducing a Lagrange multiplier \(\lambda\).

Consider a small strain elastic behaviour characterised by its free energy \(\Psi\). In the only loading is the constraint \(c\), the solution satisfies: \[ \underset{\underline{\varepsilon},\lambda}{\min}\Psi-\lambda\,c \]

In this case, the constraint \(c\) is equivalent to the following imposed stress:

\[ -\lambda\,{{\displaystyle \frac{\displaystyle \partial c}{\displaystyle \partial \underline{\varepsilon}}}} \]

If the constraint is \(\sigma_{xx}-\sigma_{0}\), where \(\sigma_{0}\) is a constant value, the previous equation shows that imposing this constraint is not equivalent to imposing an uniaxial stress state \(\left(\sigma_{xx}\,0\,0\,0\,0\,0\right)\).

[MFront Gallery] How to implement an isotropic viscoplastic behaviour with several kinematic variables following the Amstrong-Frederic evolution (27/03/2017)

The implementation of an an isotropic viscoplastic behaviour with several kinematic variables following the Amstrong-Frederic evolution law is described here

The behaviour is described by a standard split of the strain \({{\underline{\varepsilon}}^{\mathrm{to}}}\) in an elastic and a plastic parts, respectively denoted \({{\underline{\varepsilon}}^{\mathrm{el}}}\) and \({{\underline{\varepsilon}}^{\mathrm{vis}}}\):

\[ {{\underline{\varepsilon}}^{\mathrm{to}}}={{\underline{\varepsilon}}^{\mathrm{el}}}+{{\underline{\varepsilon}}^{\mathrm{vis}}}\]

Elastic behaviour

The stress \({\underline{\sigma}}\) is related to the the elastic strain \({{\underline{\varepsilon}}^{\mathrm{el}}}\) by a the standard Hooke behaviour:

\[ {\underline{\sigma}}= \lambda\,{{\mathrm{tr}{{\left({{\underline{\varepsilon}}^{\mathrm{el}}}\right)}}}}\,{\underline{I}}+2\,\mu\,{{\underline{\varepsilon}}^{\mathrm{el}}}\]

Viscoplastic behaviour

The viscoplastic behaviour follows a standard viscoplastic behaviour: \[ {{\underline{\dot{\varepsilon}}}^{\mathrm{vis}}}=\left\langle{{{\displaystyle \frac{\displaystyle F}{\displaystyle K}}}}\right\rangle^{m}\,{\underline{n}}=\dot{p}\,{\underline{n}} \]

where \(F\) is the yield surface defined below, \(<.>\) is Macaulay brackets, \({\underline{n}}\) is the normal to \(F\) with respect to the stress and \(p\) is the equivalent plastic strain.

The yield surface is defined by: \[ F{{\left({\underline{\sigma}},{\underline{X}}_{i},p\right)}}={{\left({\underline{\sigma}}-\sum_{i=1}^{N}{\underline{X}}_{i}\right)}}_{\mathrm{eq}}-R{{\left(p\right)}}=s^{e}_{\mathrm{eq}}-R{{\left(p\right)}} \]


We have introduced an effective deviatoric stress \({\underline{s}}^{e}\) defined by: \[ {\underline{s}}^{e}={\underline{s}}-\sum_{i=1}^{N}{\underline{X}}_{i} \] where \({\underline{s}}\) is the deviatoric part of the stress.

The normal is then given by: \[ {\underline{n}}={{\displaystyle \frac{\displaystyle \partial F}{\displaystyle \partial {\underline{\sigma}}}}}={{{\displaystyle \frac{\displaystyle 3}{\displaystyle 2}}}}\,{{{\displaystyle \frac{\displaystyle {\underline{s}}^{e}}{\displaystyle s^{e}_{\mathrm{eq}}}}}} \]

Evolution of the isotropic hardening

The isotropic hardening is defined by: \[ R{{\left(p\right)}}=R_{\infty} + {{\left(R_{0}-R_{\infty}\right)}}\,\exp{{\left(-b\,p\right)}} \]

Evolution of the kinematic hardenings

\[ {\underline{X}}_{i}={{{\displaystyle \frac{\displaystyle 2}{\displaystyle 3}}}}\,C_{i}\,{\underline{a}}_{i} \]

The evolution of the kinematic variables \({\underline{a}}_{i}\) follows the Armstrong-Frederic rule:

\[ {\underline{\dot{a}}}_{i}={{\underline{\dot{\varepsilon}}}^{\mathrm{vis}}}-g[i]\,{\underline{a}}_{i}\,\dot{p}=\dot{p}\,{{\left({\underline{n}}-g[i]\,{\underline{a}}_{i}\right)}} \]

Third MFront Users Day

CEA and EDF are pleased to announce that the third MFront users meeting will take place on May 30th 2017 and will be organized by CEA DEC/SESC at CEA Saclay in the DIGITEO building (building 565, room 34).

Researchers and engineers willing to present their works are welcome.

They may contact the organizers for information at: .

Registration by sending an email at the previous address before May 12th is required to ensure proper organization.

Announcing the tfel-plot project

The tfel-plot project is meant to create:

This project is based on:

Compared to gnuplot, we wanted to introduced the following new features:

More details can be found on the dedicated github page:

[MFront Gallery] Description of the implementation of a simple orthotropic behaviour (27/01/2017)

A new example has been added in the gallery here.

This example describes a simple orthotropic behaviour.

The behaviour is described by a standard split of the strain \({{\underline{\varepsilon}}^{\mathrm{to}}}\) in an elastic and a plastic parts, respectively denoted \({{\underline{\varepsilon}}^{\mathrm{el}}}\) and \({{\underline{\varepsilon}}^{\mathrm{p}}}\):

\[ {{\underline{\varepsilon}}^{\mathrm{to}}}={{\underline{\varepsilon}}^{\mathrm{el}}}+{{\underline{\varepsilon}}^{\mathrm{p}}}\]

Elastic behaviour

The stress \({\underline{\sigma}}\) is related to the the elastic strain \({{\underline{\varepsilon}}^{\mathrm{el}}}\) by a the orthotropic elastic stiffness \({\underline{\underline{\mathbf{D}}}}\):

\[ {\underline{\sigma}}= {\underline{\underline{\mathbf{D}}}}\,\colon\,{{\underline{\varepsilon}}^{\mathrm{el}}}\]

The plastic part of the behaviour is described by the following yield surface: \[ f{{\left({\sigma_{H}},p\right)}} = {\sigma_{H}}-\sigma_{0}-R\,p \]

where \({\sigma_{H}}\) is the Hill stress defined below, \(p\) is the equivalent plastic strain. \(\sigma_{0}\) is the yield stress and \(R\) is the hardening slope.

The Hill stress \({\sigma_{H}}\) is defined using the fourth order Hill tensor \(H\): \[ {\sigma_{H}}=\sqrt{{\underline{\sigma}}\,\colon\,{\underline{\underline{\mathbf{H}}}}\colon\,{\underline{\sigma}}} \]

The plastic flow is assumed to be associated, so the flow direction \({\underline{n}}\) is given by \({{\displaystyle \frac{\displaystyle \partial f}{\displaystyle \partial {\underline{\sigma}}}}}\):

\[ {\underline{n}} = {{\displaystyle \frac{\displaystyle \partial f}{\displaystyle \partial {\underline{\sigma}}}}} = {{{\displaystyle \frac{\displaystyle 1}{\displaystyle {\sigma_{H}}}}}}\,{\underline{\underline{\mathbf{H}}}}\,\colon\,{\underline{\sigma}}\]

New eigensolvers (23/01/2017)

The default eigen solver for symmetric tensors used in TFEL is based on analitical computations of the eigen values and eigen vectors. Such computations are more efficient but less accurate than the iterative Jacobi algorithm (see (Kopp 2008; Kopp 2017)).

With the courtesy of Joachim Kopp, we have introduced the following algorithms:

The implementation of Joachim Kopp have been updated for C++-11 and make more generic (support of all the floatting point numbers, different types of matrix/vector objects). The algorithms have been put in a separate namespace called fses (Fast Symmetric Eigen Solver) and is independant of the rest of TFEL.

We have also introduced the Jacobi implementation of the Geometric Tools library (see (Eberly 2016; Eberly 2017)).

Those algorithms are available in 3D. For 2D symmetric tensors, we fall back to some default algorithm as described below.

Table 1: List of available eigen solvers.
Name Algorithm in 3D Algorithm in 2D
TFELEIGENSOLVER Analytical (TFEL) Analytical (TFEL)
FSESQLEIGENSOLVER QL with implicit shifts Analytical (FSES)
FSESCUPPENEIGENSOLVER Cuppen's Divide & Conquer Analytical (FSES)

The various eigen solvers available are enumerated in Table  1.

The eigen solver is passed as a template argument of the computeEigenValues or the computeEigenVectors methods as illustrated in the code below:

tmatrix<3u,3u,real> m2;
tvector<3u,real>    vp2;
Table 2: Test on \(10^{6}\) random symmetric tensors in double precision (double).
Algorithm Failure ratio \(\Delta_{\infty}\) Times (ns) Time ratio
TFELEIGENSOLVER 0.000632 7.75e-14 252663338 1
GTESYMMETRICQREIGENSOLVER 0 2.06e-15 525845499 2.08
FSESJACOBIEIGENSOLVER 0 1.05e-15 489507133 1.94
FSESQLEIGENSOLVER 0.000422 3.30e-15 367599140 1.45
FSESCUPPENEIGENSOLVER 0.020174 5.79e-15 374190684 1.48
FSESHYBRIDEIGENSOLVER 0.090065 3.53e-10 154911762 0.61
FSESANALYTICALEIGENSOLVER 0.110399 1.09e-09 157613994 0.62

Some benchmarks

We have compared the available algorithm on \(10^{6}\) random symmetric tensors whose components are in \([-1:1]\).

For a given symmetric tensor, we consider that the computation of the eigenvalues and eigenvectors failed if: \[ \Delta_{\infty}=\max_{i\in[1,2,3]}\left\|{\underline{s}}\,\cdot\,\vec{v}_{i}-\lambda_{i}\,\vec{v}_{i}\right\|>10\,\varepsilon \] where \(\varepsilon\) is the accuracy of the floatting point considered.

The results of those tests are reported on Table  2:

Official port of TFEL for FreeBSD (20/01/2017)

Thanks to the work of Pedro F. Giffuni , an official port of TFEL/MFront is available for FreeBSD:

To install an executable package, you can now simply do:

pkg install tfel-mfront

Alternativel, to build and install TFEL from the ports tree, one can do:

cd /usr/ports/science/tfel/
make install

New hyper(visco)elastic behaviours in the MFront model repository (13/01/2017)

The implementation of the Ogden behaviour is now described in depth in the following page:

This behaviour is interesting as it highlights the features introduced in TFEL-3.0 for computing isotropic functions of symmetric tensors.

This page uses the formal developments detailled in:

Concerning hyperviscoelasticity, a page describing a generic implementation is available here:

If you have particular wishes on behaviours implementation that you would like to see treated, do not hesitate to send a message at .

MFront model repository (7/01/2017)

A page referencing examples of well written mechanical behaviours has been created here:

For each behaviours, we will try to provide tutorial-like pages explaining the implementations details (usage of tensorial objects, special functions of the TFEL library, choice of the algorithms, and so on...)

The first attempt is an hyperelastic behaviour already available in Code_Aster: the Signorini behaviour. You can find the page here:

This is still under review, so corrections and feed-backs would be greatly appreciated. The following behaviours are planned to be addressed:

If you have particular wishes on behaviours implementation that you would like to see treated, do not hesitate to send a message at .

TFEL version 3.0 has been released the 16/12/2016.

From a user point of view, TFEL 3.0 brings many game-changing features:

Many improvements for mechanical behaviours have also been made:

A detailed version of the release notes is available here.

A new behaviour implementation, called UnilateralMazars (1 June 2016)

A new behaviour implementation has been submitted by F. Hamon (EDF R&D AMA). This behaviour describes the damaging behaviour of concrete with unilateral effects as described in a dedicated paper:

A new 3D damage model for concrete under monotonic, cyclic and dynamic loadings. Jacky Mazars, François Hamon and Stéphane Grange. Materials and Structures ISSN 1359-5997 DOI 10.1617/s11527-014-0439-8

This implementation can be found in the current development sources of MFront.

Progress in the Abaqus Explicit interface (VUMAT) (28 May 2016)

The Abaqus Explicit interface is becoming quite usable du to the extensive testing efforts of D. Deloison (Airbus). Here is an example a punching test (This test was also modelled using Abaqus Standard).

Documentation of behaviour bricks (27 May 2016)

@Brick "StandardElasticity";

Behaviour bricks will be one of the most important new feature of TFEL 3.0.

A dedicated page has been created here

The EUROPLEXUS interface (27 May 2016)

An interface for the EUROPLEXUS explicit finite element solver has been developed.

EUROPLEXUS (EPX) is a simulation software dedicated to the analysis of fast transient phenomena involving structures and fluids in interaction.

See the dedicated page for more information.

tfel-doc github repository. Talks of the Second MFront Users Day (27 May 2016)

A github repository has been set up to store various documents describing TFEL and MFront usage. The talks of the first and second MFront Users Days are available there:

TFEL github repository (27 May 2016)

The subversion repository used by CEA and EDF are now synchronized with a public githbub repository:

All the branches, commit description and history of TFEL are available. This repository is read-only. Its purpose is to integrate TFEL in continous-integration projects which depends on TFEL.

Announcing the mfront module (May 2016)

A new python module has been introduced to analyse MFront files.

An overview of the module is available here.

Here is a typical usage example:

import mfront

dsl = mfront.getDSL("Chaboche.mfront")

# file description
fd = dsl.getFileDescription()
print("file author:      ", fd.authorName)
print("file date:        ",
print("file descrption:\n", fd.description)

# targets information
tgt = dsl.getTargetsDescription()

# loop over (to be) generated libraries
for l in tgt:

Cast3m 2016 is out (28 April 2016)

Cast3M 2016 has been released.

This version allow even better integration with MFront and can now directly be used to make direct calls to MFront libraries for material properties (mechanical behaviours can be used since Cast3M 2015).

This syntax is now officially supported:

Ty = 'TABLE' ;
Ty.'LIB_LOI'   = '' ;
Ty.'FCT_LOI'   = 'M5_YoungModulus' ;

ma = 'MATERIAU' mo 'YOUN' Ty 'NU' 0.3;

A binary version of TFEL compiled for Cast3M 2016 is now available for download on sourceforge:

Second MFront Users Day

CEA and EDF are pleased to announce that the second MFront users meeting will take place on May 20th 2016 and will be organized by EDF R&D at the EDF Lab Paris Saclay location.

Researchers and engineers willing to present their works are welcome. They may contact the organizers for information at: .

This users day will take place on Friday, May 20th, 2016 at EDF Lab Paris-Saclay (access map).

Registration is required to ensure proper organization. See the dedicated form on the Code_Aster website:

Abaqus interface

The current development version now includes an experimental interface to the Abaqus solver through the umat subroutine.

The following results shows the results obtained on notched beam under a cyclic loading with an isotropic hardening plastic beahviour implemented with mfront:

The mfront results can be compared to the results obtained using Abaqus in-built model on the following figure:

The Abaqus interface is still in its early stage of developments. A full description of its usage and current abilities can be found in the associated documentation.

Feed-back from users would be greatly welcomed.

IMSIA Seminar about MFront on Januar, 27 2016

An IMSIA seminar about MFront will be held on Januar, 27 2016. Here is the official announcement (in french):

Séminaire IMSIA : Implémentation de lois de comportement mécanique à l’aide du générateur de code MFront

Créé dans le cadre de la simulation des éléments combustibles nucléaire au sein d’une plate-forme logicielle nommée PLEIADES, MFront est un générateur de code, distribué en open-source [1], qui vise à permettre aux ingénieurs/chercheurs d’implémenter de manière simple et rapide des lois de comportements mécaniques de manière indépendante du code cible (EF ou FTT) [2,3].

Ce dernier point permet d’échanger les lois MFront entre différents partenaires, universitaires ou industriels. Le lien vers les codes cible se fait via la notion d’interface. A l’heure actuelle, des interfaces existent pour les codes Cast3M, Code-Aster, Abaqus, ZeBuLoN, AMITEX_FFT et d’autres codes métiers, et sont en cours de développement pour d’autres codes tels Europlexus, …

Ce séminaire proposera une description des fonctionnalités de MFront (lois en transformations infinitésimales et en grandes transformations, modèles de zones cohésives) et commentera plusieurs exemples d’applications, discutera des performances numériques obtenues et soulèvera la question de la portabilité des connaissances matériau. Nous montrerons qu’une démarche cohérente allant des expérimentations aux codes de calcul est nécessaire.


[2] Introducing the open-source mfront code generator: Application to mechanical behaviours and material knowledge management within the PLEIADES fuel element modelling platform. Thomas Helfer, Bruno Michel, Jean-Michel Proix, Maxime Salvo, Jérôme Sercombe, Michel Casella, Computers & Mathematics with Applications, Volume 70, Issue 5, September 2015, Pages 994-1023, ISSN 0898-1221,

[3] Implantation de lois de comportement mécanique à l’aide de MFront : simplicité, efficacité, robustesse et portabilité. T. Helfer, J.M. Proix, O. Fandeur. 12ème Colloque National en Calcul des Structures 18-22 Mai 2015, Giens (Var)

MFront and Cyrano3 at the LWR Fuel Performance Meeting 2015 (13 - 17 September 2015, Zurich, Switzerland)

A poster describing the use of MFront in EDF Cyrano3 fuel performance code has been presented at the LWR Fuel Performance Meeting 2015 (13 - 17 September 2015, Zurich, Switzerland):

Release of version 2.0.3 (9 September 2015)

Version 2.0.3 is mostly a bug-fix release:

A full description of the 2.0.3 release can be found here (in french).

Introducting mfront-doc (19 August 2015)

mfront-doc allows the user to extract the documentation out of MFront file. mfront-doc generates files in pandoc markdown format. Those files can be processed using pandoc and be converted to one of the many file format supported by pandoc, including LaTeX, html or various Word processor formats: Microsoft Word docx, OpenOffice/LibreOffice ODT.

mfront-doc is developped in the 3.0.x branche of TFEL. A overview of the mfront-doc functionalities can be found here.

New documentations pages (18 August 2015)

New documentation pages were added to describe the MTest and MFront keywords~:

New interfaces for material properties (July 2015)

The current development version of MFront includes two new interfaces for material properties:

Here is an example of a GNU Octave session used to compute the Young Modulus of uranium-plutonium carbide \(UPuC\) for various porosities over a range of temperatures:

octave:1> T=[300:100:1500]
T =
    300    400    500    600    700    800    900   1000   1100   1200   1300   1400   1500
octave:2> y01=UPuC_YoungModulus(T,0.1)
y01 =
   1.7025e+11   1.6888e+11   1.6752e+11   1.6616e+11   1.6480e+11   1.6344e+11   1.6207e+11   1.6071e+11   1.5935e+11   1.5799e+11   1.5662e+11   1.5526e+11   1.5390e+11
octave:3> y02=UPuC_YoungModulus(T,0.2)
y02 =
   1.1853e+11   1.1758e+11   1.1663e+11   1.1568e+11   1.1474e+11   1.1379e+11   1.1284e+11   1.1189e+11   1.1094e+11   1.0999e+11   1.0905e+11   1.0810e+11   1.0715e+11

A "publications" page has been added to the website (July 2015)

The number of papers in which MFront is used is increasing. A dedicated page has been added here.

If you publish papers which refers to MFront, please consider contributing to this page.

The paper entitled "Introducing the open-source MFront code generator ..." has been accepted for publication in Computers and Mathematics with Applications (24 june 2015)

The first paper dedicated to MFront, written by Thomas Helfer, Bruno Michel, Jean-Michel Proix, Maxime Salvo, Jérôme Sercombe, and Michel Casella, has been accepted Computers and Mathematics with Applications. The paper is available online on the sciencedirect website:

The PLEIADES software environment is devoted to the thermomechanical simulation of nuclear fuel elements behaviour under irradiation. This platform is co-developed in the framework of a research cooperative program between Électricité de France (EDF), AREVA and the French Atomic Energy Commission (CEA). As many thermomechanical solvers are used within the platform, one of the PLEAIADES’s main challenge is to propose a unified software environment for capitalisation of material knowledge coming from research and development programs on various nuclear systems.

This paper introduces a tool called mfront which is basically a code generator based on C++ (Stroustrup and Eberhardt, 2004). Domain specific languages are provided which were designed to simplify the implementations of new material properties, mechanical behaviours and simple material models. mfront was recently released under the GPL open-source licence and is available on its web site:

The authors hope that it will prove useful for researchers and engineers, in particular in the field of solid mechanics. mfront interfaces generate code specific to each solver and language considered.

In this paper, after a general overview of mfront functionalities, a particular focus is made on mechanical behaviours which are by essence more complex and may have significant impact on the numerical performances of mechanical simulations. mfront users can describe all kinds of mechanical phenomena, such as viscoplasticity, plasticity and damage, for various types of mechanical behaviour (small strain or finite strain behaviour, cohesive zone models). Performance benchmarks, performed using the Code-Aster finite element solver, show that the code generated using mfront is in most cases on par or better than the behaviour implementations written in fortran natively available in this solver. The material knowledge management strategy that was set up within the PLEIADES platform is briefly discussed. A material database named sirius proposes a rigorous material verification workflow.

We illustrate the use of mfront through two case of studies: a simple FFC single crystal viscoplastic behaviour and the implementation of a recent behaviour for the fuel material which describes various phenomena: fuel cracking, plasticity and viscoplasticity.

Castem 2015 is out (12 April 2015)

Cast3M 2015 has been released.

This release allow direct call to MFront libraries for mechanical behaviours. The following syntax of the MODELISER operator is now officially supported:

   'LIB_LOI' 'src/'
   'FCT_LOI' 'umatnorton'
   'C_MATERIAU' coel2D
   'C_VARINTER' stav2D
   'PARA_LOI'   para2D
   'CONS' M;

See the dedicated page for more information.

Salome-Meca 2015.1 is out (1O April 2015)

Salome-Meca 2015.1, which combines the Salome platform and Code-Aster, has been released on 1O April 2015. This version is the first to include a pre-packaged version of TFEL (version 2.0.1).

Salome-Meca and Code_Aster Users Day (26 March 2015)

A presentation of MFront was done during the Salome-Meca and Code_Aster Users Day 2015 user meeting at EDF Lab Clamart.

Slides can be found here.

MFront user meeting (6 Februar 2015)

The first MFront user meeting was held in Cadarache on Februar,6 2015. 27 participants from CEA, EDF, Areva and CNRS could discuss and comment about their use of MFront.

Various subjects were discussed:

Those talks are available here

First MFront user meeting (12 December 2014)

We are glad to announce that the first MFront user meeting will be held in Cadarache on Februar,6 2015. Everyone is invited but a registration must be performed before Januar, 16 2015 ().

Various subjects are already planned:

Other talks are welcomed.

Some applications of MFront at the Cast3M user meeting (4 December 2014)

The Cast3M user meeting was held in Paris on November 28, 2014. Jérémy Hure had a talk about the application of MFront in finite strain elasto-plasticity. This talk is available here.

AMITEX_FFTP has its own website (4 December 2014)

AMITEX_FFTP has now its own dedicated webiste.

The main purpose of AMITEX_FFTP is to implement an efficient distributed mechanical solver based on Fast Fourier Transform. AMITEX_FFTP is developped by CEA in the Departement of Nuclear Material.

AMITEX_FFTP is distributed under a free license for research and education purpose.

Release of TFEL version 2.0.1 (2 December 2014)

TFEL version 2.0.1 is now available. This is mainly a bug-fix release after version 2.0.0.

This version is meant to be used in Code-Aster version 12.3 that will be released in January 2015.

Creation of the address (1 December 2014)

A new contact address has been created: .

This address can be used to contact directly the developers of TFEL and MFront for specific issues. However, if your issue may interest a broader audience, you may want to send a post to the TFEL users mailing lists: .

MFront talk at Materiaux 2014 Montpellier (24 November 2014)

An MFront talk was given at Materiaux 2014. Slides (in french) are available here.

Beta release of tfel-2.0.1 for Windows 64bits and Cast3M 2014 (18 November 2014)

TFEL and Cast3M 2014 on Windows 7 64 bits

TFEL and Cast3M 2014 on Windows 7 64 bits

A beta version of tfel-2.0.1 for Windows 64bits and Cast3M 2014 has been released. A binary installer is provided here.

Installing this version requires a functional installation of Cast3M \(2014\) (which shall be patched to call external libraries) and the MSYS shell (It is recommended not to install mingw compilers along with the MSYS shell as Cast3M provides its own version of those compilers).

Installation instructions of those requirements are available here.

Any feedback would be gratefully acknowledge.

Note: The binary provided requires the mingw libraries delivered with Cast3M 2014.

Note: A standalone version of tfel-2.0.1 will be provided shortly.

MFront behaviours can now be used in AMITEX_FFTP (24 October 2014)

AMITEX_FFTP is a massively parallel mechanical solver based on FFT developed at CEA. MFront behaviours can be used in AMITEX_FFTP through the UMAT interface introduced by the Cast3M finite element solver.

Polycrystals computation made with AMITEX_FFTP (1024^{3} voxels)

Polycrystals computation made with AMITEX_FFTP (\(1024^{3}\) voxels)

TFEL/MFront on Cast3M web site (15 October 2014)

A page dedicated to MFront is now available on the Cast3M web site.

The MFront page on the Cast3M web site

The MFront page on the Cast3M web site

Release of TFEL 2.0 (1 October 2014)

Here is the official announcement by Jean-Paul DEFFAIN (in French):


Une version libre de MFront est désormais officiellement disponible sous licence GPL et CECILL-A.

Cette version 2.0 permet entre autres de générer des lois de comportements en transformations infinitésimales et en grandes transformations ainsi que des modèles de zones cohésives. Les lois générées sont utilisables dans les codes aux éléments finis Cast3M, Code-Aster, ZeBuLoN, l’ensemble des applications développées dans la plateforme PLEIADES, notamment Cyrano3, et le solveur FFT de la plate-forme. Des interfaces vers d'autres codes peuvent être rajoutées.

Un projet dédié a été crée sur sourceforge ( et fournit :

Pour les personnes souhaitant contribuer au développement, le dépôt subversion est accessible sur le serveur:

L’accès à ce dépôt est ouvert à tous mais nécessite l’ouverture d’un compte spécifique sur demande au .

Nous remercions chaleureusement tous ceux qui ont contribué à cette version et invitons toutes les personnes intéressées à se joindre au développement de MFront.


Chef du programme SIMU

Commissariat à l'Énergie Atomique


Eberly, David. 2016. “A Robust Eigensolver for 3 × 3 Symmetric Matrices.”

———. 2017. “Geometric Tools.”

Kopp, Joachim. 2008. “Efficient Numerical Diagonalization of Hermitian 3x3 Matrices.” International Journal of Modern Physics C 19 (3): 523–48. doi:10.1142/S0129183108012303.

———. 2017. “Numerical Diagonalization of 3x3 Matrices.”