Using MFront Mazars Material Model

JWilliams · February 7, 2022, 9:25pm

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

*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>, <PoissonRatio>, <At>, <Bt>, <Ac>, <Bc>, <st_max>, <sc_max>, 
** <k> 
**************************************
*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

xyont · February 7, 2022, 10:10pm

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,

JWilliams · February 7, 2022, 10:41pm

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

xyont · February 7, 2022, 11:13pm

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.

JWilliams · February 9, 2022, 5:13am

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

Best wishes,
JW

xyont · February 9, 2022, 6:17am

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,

JWilliams · February 9, 2022, 11:24am

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) 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

    } else if ((smt == SECANTOPERATOR) || (smt == CONSISTENTTANGENTOPERATOR)) {
      Dt = (1 - d) * Hooke;

I am not certain how to submit a patch though.

Best wishes,
JW

xyont · February 10, 2022, 2:36am

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

Disla · July 2, 2022, 6:24pm

Good afternoon,

We have managed to set up an analysis calling mechanical behaviors deﬁned 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.

there is no recover of displacement even if the pressure is much smaller than the tensile strength.
no stress is computed;
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

** 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>, , , , , , <st_max>, <sc_max>,
**

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

JWilliams · July 4, 2022, 12:37pm

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

Disla · July 4, 2022, 1:10pm

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.

JWilliams · July 4, 2022, 1:24pm

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

Disla · July 4, 2022, 1:24pm

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.

JWilliams · July 26, 2022, 2:40am

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 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. I hope this helps.

Best wishes,
JW

Disla · July 26, 2022, 6:27am

Thanks JWilliams. !!
I will check if it helps to make my MFRONT materials work.

Disla · March 3, 2024, 10:55am

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

xyont · March 4, 2024, 8:46pm

can you show displacement results? it seems advanced material like damage ones is much sensitive to material parameter and boundary or loading scheme.

Disla · March 4, 2024, 10:55pm

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 and it’s references.

I’m using material parameters of table 1.

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.

xyont · March 4, 2024, 11:45pm

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.

i mean displacement result plot, not strain. Animation maybe?

Disla · March 5, 2024, 8:20am

You write like an auditor. .

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.

Topic		Replies	Views
Problems with Mazars' model	1	399	July 11, 2022
Output specific variables from MFront to CalculiX frd output file	2	378	February 15, 2022
CalculiX & Mfront	80	4216	September 15, 2024
MFront Orthrotropic behaviour	0	349	July 17, 2023
Does someone try WARP3D?	2	370	June 13, 2023

Using MFront Mazars Material Model

Related Topics