Hi,
It’s best to provide an example (.inp file) and Mfront library with material (portable)
Hello @mmartin,
Are you able to share the mfront file you are using to build the material definition?
Many thanks,
Jason
I just copied MFRONT.dll’s provided on this forum to my calculix folder.
And I also copied a folder named calculix with lots of *.inp files
@mmartin No worries. I found that in version 4.0 the consistent tangent operator was not calculated and this gave no results after a completed run. I thought you might have been faced with the same problem.
Best wishes,
Jason
1 Like
Thank you very much JWILLIAMS.
I can not wait to download next version.
I hope it will be ready soon.
MANUEL
New version is out Mfront 4.1.0
https://thelfer.github.io/tfel/web/release-notes-4.1.html
and msys2 (mingw, Windows) packages:
https://packages.msys2.org/base/mingw-w64-tfel
2 Likes
Hello everyone,
I’ve also a problem with implementing new behaviour with MFront.
I want to try an orthotropic plastic behaviour, that I found here How to implement an orthotropic plastic behaviour with isotropic linear hardening in MFront (sourceforge.net) .
Whole Implementation: tfel.sourceforge.net/gallery/plasticity/OrthotropicLinearHardeningPlasticity.mfront
Here’s my input data:
**
** GEOMETRY
**
*Node
1, 0., 0., 1.
2, 0., 1., 1.
3, 0., 0., 0.
4, 0., 1., 0.
5, 1., 0., 1.
6, 1., 1., 1.
7, 1., 0., 0.
8, 1., 1., 0.
*Element, type=C3D8, elset=cube
1, 5, 6, 8, 7, 1, 2, 4, 3
*Solid Section, elset=cube, material=@CALCULIXBEHAVIOUR_ORTHOTROPICLINEARHARDENINGPLASTICITY
*Nset, nset=sx2, generate
5, 8, 1
*Nset, nset=sx1, generate
1, 4, 1
*Nset, nset=sy1
1, 3, 5, 7
*Nset, nset=sy2
2, 4, 6, 8,
*Nset, nset=sz1
3, 4, 7, 8
*Nset, nset=sz2
1, 2, 5, 6,
**
** MATERIAL
**
*Material, name=@CALCULIXBEHAVIOUR_ORTHOTROPICLINEARHARDENINGPLASTICITY
*Depvar
1
*User Material, constants=17
7.8e+10,2.64233e+11,3.32e+11,
**YoungModulus1, YoungModulus2, YoungModulus3,
0.13,0.24,0.18,
**PoissonRatio12, PoissonRatio23, PoissonRatio13,
4.8e+10,1.16418e+11,7.8e+10,
**ShearModulus12, ShearModulus23, ShearModulus13,
0.371,0.629,4.052,1.5,1.5,1.5,
**HillTensor,
150e6,
**YieldStress,
150e9
**HardeningSlope
**
** LOADING
**
*Step, nlgeom=YES
*Static
0.02, 1., 1e-05, 0.2
*Boundary
sx1, 1, 1
*Boundary
sy1, 2,2
*Boundary
sz1, 3,3
*Boundary
sx2, 1, 1, 0.2
*EL PRINT, ELSET=cube
E, S
*EL FILE
S
*NODE FILE
U
*End Step
Should I have to modify the MFront data, because there are already values. And is my input data right? Especially with the depvar I’m very unsure. And why is there only one value for yieldstress although it is an orthotropic material?
CalculiX says that I have to define an orientation for orthotropic behaviour, but I think it works with global orientation, if theres nothing defined extra. The Calculation looks like this:
The numbers below are estimated upper bounds
number of:
nodes: 8
elements: 1
one-dimensional elements: 0
two-dimensional elements: 0
integration points per element: 8
degrees of freedom per node: 3
layers per element: 1
distributed facial loads: 0
distributed volumetric loads: 0
concentrated loads: 0
single point constraints: 16
multiple point constraints: 1
terms in all multiple point constraints: 1
tie constraints: 0
dependent nodes tied by cyclic constraints: 0
dependent nodes in pre-tension constraints: 0
sets: 7
terms in all sets: 27
materials: 1
constants per material and temperature: 17
temperature points per material: 2
plastic data points per material: 0
orientations: 0
amplitudes: 4
data points in all amplitudes: 4
print requests: 2
transformations: 0
property cards: 0
*INFO reading *STEP: nonlinear geometric
effects are turned on
STEP 1
Static analysis was selected
Nonlinear material laws are taken into account
Newton-Raphson iterative procedure is active
Nonlinear geometric effects are taken into account
Decascading the MPC’s
Determining the structure of the matrix:
number of equations
8
number of nonzero lower triangular matrix elements
28
increment 1 attempt 1
increment size= 2.000000e-02
sum of previous increments=0.000000e+00
actual step time=2.000000e-02
actual total time=2.000000e-02
iteration 1
Using up to 1 cpu(s) for the stress calculation.
ORTHOTROPICLINEARHARDENINGPLASTICITY:no orientation defined for an orthotropic behaviour
Did you try providing the orientation card?
I tried, but there was no difference.
I did something wrong with the *orientation card. I tried it again and now CalculiX does not complain about the orientation anymore. But it still cant pass the first iteration:
Decascading the MPC’s
Determining the structure of the matrix:
number of equations
8
number of nonzero lower triangular matrix elements
28
increment 1 attempt 1
increment size= 2.000000e-02
sum of previous increments=0.000000e+00
actual step time=2.000000e-02
actual total time=2.000000e-02
iteration 1
Using up to 1 cpu(s) for the stress calculation.
Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.
Factoring the system of equations using the symmetric spooles solver
Using up to 1 cpu(s) for spooles.
Previously, as mentions by Rafal. The source of *.mfront files need to modify to become pointer of parameter values at *.inp files.
Unfortunately, this richest material library is not maintained and supported officially.
What are the richest material libraries, and where are they. I understand the lack of support.
I mean this library (DLL) is not updated, supported and available at CalculiX official website.
Richest is almost complete material i ever found, specifically in damage type. The implementation is growth and many new materials added periodically as seen in their GitHub.
Hi Rafal,
Some time ago (2020) you published a portable version of MFront with a collection of material models which one could call from CalculiX. I guess the materials collection for CalculiX has grown and improve.
Unfortunately I don’t have any idea how to compile this again and update and I would like to kindly ask if you could provide an update of such file with the new portable TFEL/Mfront 4.2.0 or someone could assist on how to compile and update those different dll again.
Hi,
Has someone found this message on the ccx monitor when running MFRONT?
Any possible solution? I’m trying to run that old portable version but only some of them work (ej.DRUCKERPRAGER)
The behaviour ‘IncompressibleMooneyRivlinBehaviour’ has thrown a generic standard exception : tangent operator not yet available
Thanks
I will try to compile a new version by the end of this week, as far as I know, only a few materials have been added: