# Using MFront Mazars Material Model **URL:** \ **Category:** Uncategorized\ **Created:** [February 7, 2022, 9:25pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987 "2022-02-07T21:25:15Z")\ **Posts on this page:** 20\ **Page:** 1
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [February 7, 2022, 9:25pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/1 "2022-02-07T21:25:15Z")
Hello Everyone, Has anyone had any luck using the MFront Mazars model in Calculix? I am not sure how to use the @ModellingHypothesis ‘Tridimensional’ keywords in my input file. I also had a go using the Unilateral Mazars model but no luck in getting a result. To clarify, the model runs but does not give any results regarding the forces (the input file is provided below). Any advice will be greatly appreciated. Many thanks, JW ```auto *HEADING MODEL: UNILATERAL MAZARS TEST MODEL 01, DATE: 07 FEBRUARY 2022, UNIT MM_TON ** *Node 1, 0., 0., 1. 2, 0., 0., 0. 3, 0., 1., 1. 4, 0., 1., 0. 5, 1., 0., 1. 6, 1., 0., 0. 7, 1., 1., 1. 8, 1., 1., 0. *ELEMENT, TYPE=C3D8, ELSET=CUBE 1, 1, 2, 4, 3, 5, 6, 8, 7 *NSET, NSET=LOAD 3, 1, 5, 7 *NSET, NSET=SUPPORT 2, 4, 8, 6 *ELSET, ELSET=CUBE 1 ************************************** ** MATERIALS ************************************** **USER MATERIAL (MFRONT) *MATERIAL, NAME=@CALCULIXBEHAVIOUR_UNILATERALMAZARS *DEPVAR 6 *USER MATERIAL, CONSTANTS=9 30000, 0.21, 1., 1.e4, 1.25, 517, 1.1e-4, 3e-4, 0.2 ** YoungModulus>, , , , , , , , ** ************************************** *SOLID SECTION, ELSET=CUBE, MATERIAL=@CALCULIXBEHAVIOUR_UNILATERALMAZARS ************************************** *TIME POINTS,NAME=T1,GENERATE 1E-6,1.,0.025 ************************************** ** STEPS ************************************** *STEP, NLGEOM=YES, INC=1000 *STATIC, SOLVER=PARDISO 1.0, 1.0, 1e-12, 0.1 ************************************** ** BOUNDARY CONDITIONS ************************************** *BOUNDARY SUPPORT, 1, 3, 0 *BOUNDARY LOAD, 3, 3,-.0141114515308986 *NODE PRINT, NSET=SUPPORT, TIME POINTS=T1, TOTALS=ONLY RF *NODE FILE, TIME POINTS=T1 RF, U *EL FILE, TIME POINTS=T1 S, E *END STEP ``` ---
**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\ **Post date:** [February 7, 2022, 10:10pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/2 "2022-02-07T22:10:20Z")
hi, i have succeed to test and ran Mazars material models with CalculiX Mfront integrations, but for a classic ones (scalar, quasi-static) not newest two (unilateral,cyclic-dynamic). Seems the library could not properly linked, throw an error as not loaded functions. `*ERROR: invalid material name '@CALCULIXBEHAVIOUR_UNILATERALMAZARS' (unable to load function)` this problem also has been reported to @rafal.brzegowy and unfortunately did not solved yet. best, ![CalculiX_Mfront_Test1](https://global.discourse-cdn.com/free1/uploads/calculix/original/1X/90172128bb516b83fb32839ffcea47dce831bc46.png) ---
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [February 7, 2022, 10:41pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/3 "2022-02-07T22:41:34Z")
Many thanks @xyont, adding 0.0 to the second line of the user material constants got the Mazars material model case up and running. I am able to load the Unilateral Mazars but for the same general set-up as the Mazars test case, the Unilateral version runs but does not provide any results. Best wishes, JW ---
**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\ **Post date:** [February 7, 2022, 11:13pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/4 "2022-02-07T23:13:18Z")
thanks for the hints at nine of constants, for now is loaded and running but your example problems are divergence using PaStiX solver. i do testing also on my example, it’s finished to ran but shown no results also as you mentions. ![CalculiX_Mfront_Test2](https://global.discourse-cdn.com/free1/uploads/calculix/original/1X/bf9c14ea6b3e9714bd2284cae4ea47432086cea9.png) ---
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [February 9, 2022, 5:13am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/5 "2022-02-09T05:13:16Z")
Hi @xyont, I got the Unilateral Mazars material model up and running today but using the UniliateralMazars.mfront found at the link below. I was able to run the example given by “A new 3D damage model for concrete under monotonic, cyclic and dynamic loadings” (Mazars et al 2015) on a unit cube. I have attached a picture of the stress-strain plot. The run was done using a Dynamic, Explicit step. [TFEL/MFRONT Discussion](https://sourceforge.net/p/tfel/discussion/general/thread/e7cffa8942/) Best wishes, JW ![Plot01](https://global.discourse-cdn.com/free1/uploads/calculix/original/1X/c24704748622ce2bd3c28622e8c069375ff154f8.jpeg) ---
**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\ **Post date:** [February 9, 2022, 6:17am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/6 "2022-02-09T06:17:04Z")
hi @JWilliams great news it’s working. how speed is compared to Mazars scalar at monotonic loading? i’m curious to try and know even for simple models. did you deploy or recompile an update of UnilateralMazars from the discussion (Helfer Thomas) and generate dll by yourself ? submit a patch to Rafal and distributed officially can be a great contribution to completeness and rich of CalculiX material libraries. best, ---
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [February 9, 2022, 11:24am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/7 "2022-02-09T11:24:38Z")
Hi @xyont, Sorry, I have not made any comparisons to Mazars scalar. With help from @rafal.brzegowy, I was able to deploy and recompile UnilateralMazars on my machine. Based on a query raised in the TFEL/Mfront discussion ([TFEL Discussion 02](https://sourceforge.net/p/tfel/discussion/mechanicalbeahviours/thread/deb3321424/)) the version found in tfel.4.0 does not compute the consistent tangent operator and fails when called using CalculiX. I believe the line of code below is required ```auto } else if ((smt == SECANTOPERATOR) || (smt == CONSISTENTTANGENTOPERATOR)) { Dt = (1 - d) * Hooke; ``` I am not certain how to submit a patch though. Best wishes, JW ---
**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\ **Post date:** [February 10, 2022, 2:36am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/8 "2022-02-10T02:36:07Z")
hi @JWilliams thanks again for the hints of caused problems, reading the discussion you linked it’s a patch from the main developer (Helfer Thomas) glad he is helping. unfortunately, i’m not have any experienced in deploy or compiled MFront libraries. may i start to learning in deploying or waiting in hopefully @rafal.brzegowy to have a time and could give an updates the library of working ones. best regards, sn ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [July 2, 2022, 6:24pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/9 "2022-07-02T18:24:29Z")
Good afternoon, We have managed to set up an analysis calling mechanical behaviors defined in shared libraries (@). The inp contains material models Mazars, UnilatteralMazzars and an elastic isotropic cube just for comparison. Results don’t seem to have any sense for different reasons. 1. there is no recover of displacement even if the pressure is much smaller than the tensile strength. 2. no stress is computed; 3. no difference in displacements is obtained if the Young Modulus and the tensile / compression strength are specified in MPa instead of Pa. In fact, changing material parameters don’t change the result except for density which need to be insanely high to get some comparable result. We hare using the ccx compilation provided with the Mecway distribution as the one from Dhont give the error: \*ERROR: invalid material name ‘@CALCULIXBEHAVIOUR\_UNILATERALMAZARS@UNILATERALMAZARS’ (unable to load function) ¿Is there any advance about the portable version of MFRONT for windows? Thanks ![screenshot.22](https://global.discourse-cdn.com/free1/uploads/calculix/original/1X/e8886206a8d598091dbb03007dbead481eb1ac73.png) \*\* Generated by Mecway 15.0 \*NODE 1,0.2,-0.4,0 2,0.8,0.7,0.6 3,1,-0.6,0.2 4,0.8,-0.8,0.4 5,0.4,-2.1,0 6,0,-0.4,0 7,0.4,-0.8,0.4 8,0,-1.7,0.6 9,0.4,-0.8,0.2 10,0,-2.1,0.4 11,0,-0.2,0.6 12,0,-0.6,0.6 13,0,-0.6,0.4 14,0,-0.6,0.2 15,0.8,1.1,0.6 16,0,-0.6,0 17,0,-0.2,0.4 18,0,-0.4,0.6 19,0.6,1.1,0.4 20,1,-0.6,0.6 21,1,-2.3,0.2 22,1,-0.6,0 23,0,-1.7,0.2 24,0.4,-0.2,0 25,0.6,-0.2,0 26,0.8,-0.2,0.2 27,0,-1.7,0.4 28,0.4,1.3,0.6 29,1,-0.6,0.4 30,0.6,0.9,0.4 31,1,-2.3,0.4 32,0,-0.4,0.4 33,0.4,1.3,0.4 34,0,-0.8,0.6 35,0.6,0.7,0.6 36,0,-2.3,0.6 37,1,-2.1,0.2 38,0.4,1.3,0.2 39,0.8,-0.8,0.2 40,0,1.3,0.6 41,0,-1.9,0.6 42,0,1.3,0.4 43,0,-0.8,0.4 44,0,1.3,0.2 45,0.6,-0.4,0.2 46,0,1.3,0 47,1,1.3,0.6 48,1,-2.1,0.4 49,0.6,0.7,0 50,0.4,1.3,0 51,0.8,-0.4,0.2 52,0.8,0.9,0.2 53,0.8,1.1,0.2 54,1,-1.9,0.2 55,0.6,-1.9,0 56,0.6,-1.9,0.4 57,0,1.1,0.6 58,0,-2.1,0.6 59,0,1.1,0.4 60,0,-0.8,0.2 61,0,1.1,0.2 62,0.8,0.7,0.4 63,0,1.1,0 64,0.8,-0.8,0 65,1,-1.9,0.4 66,0.4,-0.8,0 67,0,-0.4,0.2 68,0.4,-2.1,0.4 69,0.8,-0.4,0 70,0,-0.2,0.2 71,1,-1.7,0.2 72,0.4,1.1,0.6 73,0,-0.8,0 74,0,0.9,0.6 75,0.6,-0.6,0.6 76,0,0.9,0.4 77,0.4,-2.3,0.2 78,0,0.9,0.2 79,0.4,1.1,0.4 80,0,0.9,0 81,0.8,-0.6,0.2 82,1,-1.7,0.4 83,1,-0.2,0.2 84,0.4,1.1,0.2 85,1,0.9,0.6 86,1,0.9,0 87,0.6,-0.8,0.2 88,0.8,1.3,0.2 89,0.4,1.1,0 90,0.2,-1.7,0.2 91,0,0.7,0.6 92,0.2,-0.6,0.6 93,0,0.7,0.4 94,0.6,1.1,0.2 95,0,0.7,0.2 96,0,-0.2,0 97,0,0.7,0 98,0.8,-0.2,0.4 99,0.2,-2.1,0.6 100,0.8,-0.2,0 101,0.6,-1.7,0.6 102,1,1.3,0.2 103,1,-0.2,0.6 104,0.2,-0.6,0.4 105,0.4,-0.4,0.6 106,1,-0.8,0.6 107,0.4,-0.4,0.4 108,0.6,-0.8,0 109,0.2,-1.7,0.6 110,0.6,1.3,0.6 111,0.4,-0.4,0.2 112,0.6,0.7,0.4 113,1,-2.3,0 114,0.6,-1.7,0 115,1,-2.3,0.6 116,0.8,-0.6,0.4 117,0.2,-2.1,0 118,0.6,0.9,0.2 119,1,-2.1,0 120,0.8,-0.6,0 121,1,-2.1,0.6 122,0.6,-0.2,0.6 123,0.4,-0.4,0 124,0.2,-1.9,0 125,1,-1.9,0 126,0.6,-0.8,0.4 127,1,-1.9,0.6 128,0.2,-2.3,0.2 129,0.6,-2.3,0.4 130,0,-2.3,0.2 131,1,-1.7,0 132,0.2,-2.3,0.4 133,1,-1.7,0.6 134,1,0.7,0.2 135,0.2,-1.9,0.2 136,0.4,-2.3,0 137,0.4,0.9,0.6 138,0.8,-0.4,0.6 139,0.6,1.3,0.4 140,0.4,-1.7,0.4 141,0.2,-1.7,0.4 142,1,1.1,0.4 143,1,1.3,0.4 144,0.6,-0.4,0.4 145,0.6,1.3,0.2 146,0.6,0.9,0 147,0.2,-0.6,0.2 148,0.6,-0.6,0.2 149,0.4,0.9,0.4 150,1,1.1,0 151,1,-0.8,0 152,0,-1.9,0 153,0.6,-0.4,0.6 154,0.4,0.9,0.2 155,1,-0.8,0.2 156,0.2,-2.1,0.2 157,0.6,-0.2,0.4 158,0.4,0.9,0 159,1,-0.4,0.4 160,1,-0.2,0.4 161,0.2,-2.3,0 162,1,1.1,0.2 163,0.6,-0.2,0.2 164,0.4,-1.9,0.4 165,0.2,-1.9,0.4 166,0.2,-0.6,0 167,0.6,-0.6,0 168,0.8,0.9,0.6 169,0.6,-2.1,0.4 170,0.2,1.3,0.6 171,0,-2.1,0 172,0.2,1.3,0.4 173,0,-1.9,0.4 174,0.2,1.3,0.2 175,0.2,-1.7,0 176,0.2,1.3,0 177,0.8,1.1,0.4 178,1,-0.4,0 179,0.6,1.1,0 180,0.4,-0.6,0.6 181,0.6,-2.3,0.2 182,0.6,-1.9,0.2 183,0.8,0.9,0 184,0.8,-2.3,0.2 185,0.4,0.7,0.6 186,0,-2.3,0.4 187,0.2,1.1,0.6 188,0.4,0.7,0.4 189,0.2,1.1,0.4 190,0.4,-0.2,0.6 191,0.2,1.1,0.2 192,0.4,-0.6,0.4 193,0.2,1.1,0 194,1,0.9,0.2 195,0.8,-2.3,0.4 196,0.4,-0.6,0.2 197,0.6,-2.1,0 198,0,-2.3,0 199,0.4,-0.6,0 200,0.4,0.7,0.2 201,0.8,-2.1,0.2 202,0.4,-2.3,0.4 203,0.6,0.7,0.2 204,0.2,0.9,0.6 205,0.4,0.7,0 206,0.2,0.9,0.4 207,0.8,1.1,0 208,0.2,0.9,0.2 209,0.4,-0.2,0.4 210,0.2,0.9,0 211,0.6,-2.3,0.6 212,0.8,-2.1,0.4 213,1,-0.4,0.2 214,0.8,0.7,0.2 215,0.6,-1.7,0.4 216,0.6,-2.1,0.2 217,0.6,-1.9,0.6 218,0.8,-1.9,0.4 219,0.6,-2.1,0.6 220,0.2,-0.8,0.6 221,0.2,0.7,0.6 222,0.6,-2.3,0 223,0.2,0.7,0.4 224,0.8,-1.7,0.2 225,0.2,0.7,0.2 226,0.8,-0.8,0.6 227,0.2,0.7,0 228,0.4,-1.9,0.6 229,0,-1.7,0 230,0.4,-0.2,0.2 231,1,1.3,0 232,0.4,-2.1,0.2 233,0.2,-0.8,0.4 234,0.8,-0.6,0.6 235,0.8,-1.9,0.2 236,0.8,-0.2,0.6 237,0,-2.1,0.2 238,0.2,-0.8,0.2 239,0.4,-1.7,0.2 240,0.2,-2.3,0.6 241,0.8,-1.7,0 242,0.2,-0.8,0 243,0.4,-1.9,0.2 244,0.2,-0.2,0.6 245,0.6,1.1,0.6 246,0.2,-2.1,0.4 247,0.8,0.7,0 248,0.2,-0.2,0.4 249,0.8,1.3,0.6 250,1,1.1,0.6 251,0.2,-0.2,0.2 252,1,-0.2,0 253,0.6,-1.7,0.2 254,0.4,-2.1,0.6 255,0.6,1.3,0 256,1,-0.4,0.6 257,0.8,-2.3,0 258,0.2,-0.2,0 259,0.8,-2.3,0.6 260,0.6,0.9,0.6 261,0.8,1.3,0.4 262,0.8,-0.4,0.4 263,0.8,-2.1,0 264,0.2,-1.9,0.6 265,0.8,-2.1,0.6 266,0.6,-0.4,0 267,1,0.7,0.6 268,1,-0.8,0.4 269,0.6,-0.6,0.4 270,0.8,-1.9,0.6 271,1,0.7,0.4 272,0.4,-1.7,0 273,0.8,-1.7,0.6 274,0.8,0.9,0.4 275,0.8,-1.9,0 276,1,0.7,0 277,0,-1.9,0.2 278,0.4,-2.3,0.6 279,0.8,-1.7,0.4 280,1,0.9,0.4 281,0.2,-0.4,0.6 282,0.4,-1.9,0 283,0.2,-0.4,0.4 284,0.4,-1.7,0.6 285,0.2,-0.4,0.2 286,0.6,-0.8,0.6 287,0.4,-0.8,0.6 288,0.8,1.3,0 \*ELEMENT,TYPE=C3D8 1,81,3,155,39,120,22,151,64 2,23,90,135,277,229,175,124,152 3,164,56,169,68,243,182,216,232 4,236,103,256,138,98,160,159,262 5,98,160,159,262,26,83,213,51 6,26,83,213,51,100,252,178,69 7,243,182,216,232,282,55,197,5 8,41,264,99,58,173,165,246,10 9,40,170,187,57,42,172,189,59 10,42,172,189,59,44,174,191,61 11,44,174,191,61,46,176,193,63 12,173,165,246,10,277,135,156,237 13,277,135,156,237,152,124,117,171 14,57,187,204,74,59,189,206,76 15,59,189,206,76,61,191,208,78 16,61,191,208,78,63,193,210,80 17,169,212,195,129,216,201,184,181 18,216,201,184,181,197,263,257,222 19,74,204,221,91,76,206,223,93 20,76,206,223,93,78,208,225,95 21,78,208,225,95,80,210,227,97 22,58,99,240,36,10,246,132,186 23,10,246,132,186,237,156,128,130 24,11,244,281,18,17,248,283,32 25,17,248,283,32,70,251,285,67 26,70,251,285,67,96,258,1,6 27,237,156,128,130,171,117,161,198 28,254,219,211,278,68,169,129,202 29,18,281,92,12,32,283,104,13 30,32,283,104,13,67,285,147,14 31,67,285,147,14,6,1,166,16 32,68,169,129,202,232,216,181,77 33,138,256,20,234,262,159,29,116 34,170,28,72,187,172,33,79,189 35,172,33,79,189,174,38,84,191 36,174,38,84,191,176,50,89,193 37,262,159,29,116,51,213,3,81 38,51,213,3,81,69,178,22,120 39,187,72,137,204,189,79,149,206 40,189,79,149,206,191,84,154,208 41,191,84,154,208,193,89,158,210 42,116,29,268,4,81,3,155,39 43,232,216,181,77,5,197,222,136 44,204,137,185,221,206,149,188,223 45,206,149,188,223,208,154,200,225 46,208,154,200,225,210,158,205,227 47,217,270,265,219,56,218,212,169 48,273,133,127,270,279,82,65,218 49,12,92,220,34,13,104,233,43 50,13,104,233,43,14,147,238,60 51,14,147,238,60,16,166,242,73 52,279,82,65,218,224,71,54,235 53,56,218,212,169,182,235,201,216 54,122,236,138,153,157,98,262,144 55,157,98,262,144,163,26,51,45 56,163,26,51,45,25,100,69,266 57,105,153,75,180,107,144,269,192 58,109,284,228,264,141,140,164,165 59,28,110,245,72,33,139,19,79 60,33,139,19,79,38,145,94,84 61,38,145,94,84,50,255,179,89 62,141,140,164,165,90,239,243,135 63,90,239,243,135,175,272,282,124 64,72,245,260,137,79,19,30,149 65,79,19,30,149,84,94,118,154 66,84,94,118,154,89,179,146,158 67,107,144,269,192,111,45,148,196 68,111,45,148,196,123,266,167,199 69,137,260,35,185,149,30,112,188 70,149,30,112,188,154,118,203,200 71,154,118,203,200,158,146,49,205 72,264,228,254,99,165,164,68,246 73,165,164,68,246,135,243,232,156 74,153,138,234,75,144,262,116,269 75,144,262,116,269,45,51,81,148 76,45,51,81,148,266,69,120,167 77,135,243,232,156,124,282,5,117 78,235,54,37,201,275,125,119,263 79,244,190,105,281,248,209,107,283 80,248,209,107,283,251,230,111,285 81,251,230,111,285,258,24,123,1 82,182,235,201,216,55,275,263,197 83,99,254,278,240,246,68,202,132 84,110,249,15,245,139,261,177,19 85,139,261,177,19,145,88,53,94 86,145,88,53,94,255,288,207,179 87,246,68,202,132,156,232,77,128 88,156,232,77,128,117,5,136,161 89,245,15,168,260,19,177,274,30 90,19,177,274,30,94,53,52,118 91,94,53,52,118,179,207,183,146 92,180,75,286,287,192,269,126,7 93,265,121,115,259,212,48,31,195 94,260,168,2,35,30,274,62,112 95,30,274,62,112,118,52,214,203 96,118,52,214,203,146,183,247,49 97,190,122,153,105,209,157,144,107 98,209,157,144,107,230,163,45,111 99,281,105,180,92,283,107,192,104 100,283,107,192,104,285,111,196,147 101,285,111,196,147,1,123,199,166 102,230,163,45,111,24,25,266,123 103,224,71,54,235,241,131,125,275 104,92,180,287,220,104,192,7,233 105,104,192,7,233,147,196,9,238 106,147,196,9,238,166,199,66,242 107,212,48,31,195,201,37,21,184 108,192,269,126,7,196,148,87,9 109,249,47,250,15,261,143,142,177 110,261,143,142,177,88,102,162,53 111,88,102,162,53,288,231,150,207 112,270,127,121,265,218,65,48,212 113,218,65,48,212,235,54,37,201 114,15,250,85,168,177,142,280,274 115,177,142,280,274,53,162,194,52 116,53,162,194,52,207,150,86,183 117,219,265,259,211,169,212,195,129 118,101,273,270,217,215,279,218,56 119,168,85,267,2,274,280,271,62 120,274,280,271,62,52,194,134,214 121,52,194,134,214,183,86,276,247 122,284,101,217,228,140,215,56,164 123,140,215,56,164,239,253,182,243 124,75,234,226,286,269,116,4,126 125,269,116,4,126,148,81,39,87 126,148,81,39,87,167,120,64,108 127,196,148,87,9,199,167,108,66 128,201,37,21,184,263,119,113,257 129,239,253,182,243,272,114,55,282 130,215,279,218,56,253,224,235,182 131,234,20,106,226,116,29,268,4 132,253,224,235,182,114,241,275,55 133,228,217,219,254,164,56,169,68 134,8,109,264,41,27,141,165,173 135,27,141,165,173,23,90,135,277 \*ELSET,ELSET=UnilateralMazars 9 10 11 14 15 16 19 20 21 34 35 36 39 40 41 44 45 46 59 60 61 64 65 66 69 70 71 84 85 86 89 90 91 94 95 96 109 110 111 114 115 116 119 120 121 \*ELSET,ELSET=Mazars 1 4 5 6 24 25 26 29 30 31 33 37 38 42 49 50 51 54 55 56 57 67 68 74 75 76 79 80 81 92 97 98 99 100 101 102 104 105 106 108 124 125 126 127 131 \*ELSET,ELSET=EALL 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 \*ELSET,ELSET=Elastic 2 3 7 8 12 13 17 18 22 23 27 28 32 43 47 48 52 53 58 62 63 72 73 77 78 82 83 87 88 93 103 107 112 113 117 118 122 123 128 129 130 132 133 134 135 \*SURFACE,NAME=pressure\_faces 109,S4 114,S4 119,S4 120,S4 115,S4 110,S4 5,S4 4,S4 131,S4 113,S4 112,S4 107,S4 93,S4 52,S4 48,S4 42,S4 37,S4 33,S4 \*MATERIAL,NAME=Material \*ELASTIC,TYPE=ISOTROPIC 34180000000,0.2 \*DENSITY 1000 \*SOLID SECTION,ELSET=Elastic,MATERIAL=Material \*BOUNDARY 6,1,0 6,2,0 6,3,0 8,1,0 8,2,0 8,3,0 10,1,0 10,2,0 10,3,0 11,1,0 11,2,0 11,3,0 12,1,0 12,2,0 12,3,0 13,1,0 13,2,0 13,3,0 14,1,0 14,2,0 14,3,0 16,1,0 16,2,0 16,3,0 17,1,0 17,2,0 17,3,0 18,1,0 18,2,0 18,3,0 23,1,0 23,2,0 23,3,0 27,1,0 27,2,0 27,3,0 32,1,0 32,2,0 32,3,0 34,1,0 34,2,0 34,3,0 36,1,0 36,2,0 36,3,0 40,1,0 40,2,0 40,3,0 41,1,0 41,2,0 41,3,0 42,1,0 42,2,0 42,3,0 43,1,0 43,2,0 43,3,0 44,1,0 44,2,0 44,3,0 46,1,0 46,2,0 46,3,0 57,1,0 57,2,0 57,3,0 58,1,0 58,2,0 58,3,0 59,1,0 59,2,0 59,3,0 60,1,0 60,2,0 60,3,0 61,1,0 61,2,0 61,3,0 63,1,0 63,2,0 63,3,0 67,1,0 67,2,0 67,3,0 70,1,0 70,2,0 70,3,0 73,1,0 73,2,0 73,3,0 74,1,0 74,2,0 74,3,0 76,1,0 76,2,0 76,3,0 78,1,0 78,2,0 78,3,0 80,1,0 80,2,0 80,3,0 91,1,0 91,2,0 91,3,0 93,1,0 93,2,0 93,3,0 95,1,0 95,2,0 95,3,0 96,1,0 96,2,0 96,3,0 97,1,0 97,2,0 97,3,0 130,1,0 130,2,0 130,3,0 152,1,0 152,2,0 152,3,0 171,1,0 171,2,0 171,3,0 173,1,0 173,2,0 173,3,0 186,1,0 186,2,0 186,3,0 198,1,0 198,2,0 198,3,0 229,1,0 229,2,0 229,3,0 237,1,0 237,2,0 237,3,0 277,1,0 277,2,0 277,3,0 \*AMPLITUDE,NAME=A\_1 0,0 1,1000000 \*TIME POINTS,NAME=timepointsname,GENERATE 0,1,0.1 * * * \*\* MATERIALS: Mazars * * * \*\*USER MATERIAL (MFRONT) \*MATERIAL, NAME=@CALCULIXBEHAVIOUR\_MAZARS@Mazars \*DEPVAR 3 \*USER MATERIAL, constants=8 \*\*CALCULIX MAZARS \*\*, , , , , , , \*\* 23500e6, 0.2, 1.13, 0.7 , 1250e6,8000e6, 1, 0.000067 0.7 * * * \*SOLID SECTION, ELSET=Mazars, MATERIAL=@CALCULIXBEHAVIOUR\_MAZARS@Mazars * * * \*DENSITY 10000000000 * * * \*\* MATERIALS: UnilateralMazars * * * \*\*USER MATERIAL (MFRONT) \*MATERIAL, NAME=@CALCULIXBEHAVIOUR\_UNILATERALMAZARS@UnilateralMazars \*DEPVAR 6 \*USER MATERIAL, CONSTANTS=9 \*\*UNILATERALMAZARS \*\* YoungModulus\>, , , , , , \, \, \*\* 34180e6, 0.2, 0.9, 12200e6, 0.07, 101e6, 50e6, -300e6, 0.7 * * * \*SOLID SECTION, ELSET=UnilateralMazars, MATERIAL=@CALCULIXBEHAVIOUR\_UNILATERALMAZARS@UnilateralMazars * * * \*DENSITY 10000000000 \*STEP,NLGEOM=YES,INC=100,AMPLITUDE=STEP,SOLVER=PARDISO \*DYNAMIC 0.1,1,0,0.1 \*DSLOAD,AMPLITUDE=A\_1 pressure\_faces,P,1 \*NODE FILE,GLOBAL=YES U,RF \*EL FILE,OUTPUT=3D S,E,ENER \*ELPRINT,ELSET=EALL S,PEEQ \*END STEP ---
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [July 4, 2022, 12:37pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/10 "2022-07-04T12:37:21Z")
Hello @Disla, Can you have a go at the test model included at the top of this thread and see if it works? In my tests the material name was given by: > @CALCULIXBEHAVIOUR\_UNILATERALMAZARS or @CALCULIXBEHAVIOUR\_MAZARS But I think you have some extra text in your material definition: > @CALCULIXBEHAVIOUR\_UNILATERALMAZARS@UnilateralMazars or @CALCULIXBEHAVIOUR\_MAZARS@Mazars Best wishes, JW ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [July 4, 2022, 1:10pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/11 "2022-07-04T13:10:07Z")
Hello Williams, Thanks for giving attention to this. If I completely remove @… it runs but provides those same weird results. According to the manual a second @ caracter could be used to call the same behavior a second time but with different properties. I have used your material properties (numerical values) as reference too. ---
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [July 4, 2022, 1:24pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/12 "2022-07-04T13:24:12Z")
Hello @Disla, In my case when I got that unable to load message, I compiled both CalculiX and the MFRONT library on my machine to get everything to work 😅. I am not certain if anyone else has a much quicker fix. Best wishes, JW ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [July 4, 2022, 1:24pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/13 "2022-07-04T13:24:41Z")
Hi again, Maybe a silly comment but I have just noticed that ccx manual say that to call the material behavior one should write: @CALCULIXBEHAVIOURS\_ … With a final S Portable dll (libCALCULIXBEHAVIOUR.dll) and all the included examples .inp files in the zip refer to : @CALCULIXBEHAVIOUR\_… Without final S. ¿Could this be the reason why the library is not found or the result is nonsese ? I’m sorry but I do not have programing skills to inspect the code in this regards. ---
**Author:** ![JWilliams](https://avatars.discourse-cdn.com/v4/letter/j/c2a13f/32.png) [@JWilliams](https://calculix.discourse.group/u/JWilliams)\ **Post date:** [July 26, 2022, 2:40am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/14 "2022-07-26T02:40:07Z")
Hello @Disla and @xyont, I jumped to another machine today and had the same `unable to load` issue with the MFRONT material. In order to get things up and running again, I had to download the source files for MinGW-w64 from here [Releases · niXman/mingw-builds-binaries · GitHub](https://github.com/niXman/mingw-builds-binaries/releases) and add the bin folder to my path. In my case I downloaded [x86\_64-12.1.0-release-win32-sjlj-rt\_v10-rev3.7z](https://github.com/niXman/mingw-builds-binaries/releases/download/12.1.0-rt_v10-rev3/x86_64-12.1.0-release-win32-sjlj-rt_v10-rev3.7z). I hope this helps. Best wishes, JW ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [July 26, 2022, 6:27am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/15 "2022-07-26T06:27:56Z")
Thanks JWilliams. !! I will check if it helps to make my MFRONT materials work. ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [March 3, 2024, 10:55am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/16 "2024-03-03T10:55:37Z")
Hi JWilliams, I got the Unilateral Mazars material model up and running too following some of the posts recomendation. Yours in-between them . Thanks in advance. May I ask you how did you get VM and Strain from the model?. I can read demage and all related StateVariables but not VM or Principal strains. I would like to build the equivalent Strain for cracking and crushing but e1 e2 and e3 are not available or they haven’t been printed. Did you compute them from the input load or did you impose them by means of a displacement driven test? Now I feel like partially blind when postprocessing the results. Thanks ---
**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\ **Post date:** [March 4, 2024, 8:46pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/17 "2024-03-04T20:46:53Z")
can you show displacement results? it seems advanced material like damage ones is much sensitive to material parameter and boundary or loading scheme. ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [March 4, 2024, 10:55pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/18 "2024-03-04T22:55:55Z")
I’m more reading here and there now. Main driver Jacky Mazars • Franc¸ois Hamon • Stephane Grange : [A new 3D damage model for concrete under monotonic, cyclic and dynamic loadings](https://rdcu.be/dz2eG) and it’s references. I’m using material parameters of table 1. > [@xyont](#): > > can you show displacement results? Model is a 1m3 cube. Displacement x just before damage starts; 0.390474 mm Final Strain Sxx at the end of the load process (Damage =0.976132 ) : 0.0070 (Langrange Strain). I’m finding some weird things with the setting of the time step, there is no VM and just one Principal Strain. ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/7/70159f5099cdf350307c14679ec5eb4d4df8017c.png) ---
**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\ **Post date:** [March 4, 2024, 11:45pm UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/19 "2024-03-04T23:45:43Z")
> [@Disla](#): > > I’m using material parameters of table 1. is incomplete, what value of ‘k’ being used and derive from? These are related to shear, i have discussed to the author of TFEL/MFront in the past but can not remember exactly. > [@Disla](#): > > Displacement x just before damage starts; 0.390474 mm i mean displacement result plot, not strain. Animation maybe? ---
**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\ **Post date:** [March 5, 2024, 8:20am UTC](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987/20 "2024-03-05T08:20:44Z")
> [@xyont](#): > > is incomplete, what value of ‘k’ being used and derive from? You write like an auditor. 😂. > [@xyont](#): > > i have discussed to the author of TFEL/MFront in the past but can not remember exactly Pleaese, try to recover what was your question to the author and his response. It could help to progress on the model set up. My model is at an early stage now. I’m focus understanding material behavior , how the model responds to the different inputs and if those beautiful colors have any sense. The particular value and how to derive it from experimental data is not relevant to me now. That is something up to the user and the application. What I already know about k is that the original Mazars underestimate the resistance of the concrete in pure shearing. K parameter is introduced to fix that, increasing its resistance (it introduces an asymptote in the Stress Strain Curve in shearing). The value of k lower than 1 is very useful for modeling the effects of friction of reinforcements in reinforced concrete because it keeps a residual shear stress even when the model is completely damaged. That gave me a possible comparison parameter to check as for k=1, one should recovers the behavior of the original Mazars response under pure shear. The value advised for k is 0.7. [Next page](https://calculix.discourse.group/t/using-mfront-mazars-material-model/987.md?page=2)