Using MFront Mazars Material Model

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

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,

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

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.

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

1 Like

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,

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

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

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

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

1 Like

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

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.

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 :sweat_smile:. I am not certain if anyone else has a much quicker fix.

Best wishes,
JW

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.