I have observed an unusual phenomenon: when the Mohr-Coulomb constitutive model is combined with linear surface-to-surface contact in CalculiX, the model fails to converge. After replacing the linear contact definition with the *Tie constraint, the model runs successfully, yet no plastic zones are developed in the calculation results.
Today I ran a test and found that after activating the *SURFACE INTERACTION keyword (even without defining any corresponding contact pairs), the analysis appears to completely bypass the Mohr-Coulomb constitutive model, yielding zero equivalent plastic strain. The Mohr-Coulomb plasticity only works and plastic strain emerges after removing all *SURFACE INTERATION keywords. I would like to clarify whether this is a software bug or a user setup error on my part.
My test files are stored in Google Drive. Pay close attention to the *INCLUDE keyword—this section needs modification. Whenever I add a “Surface interaction” definition, the MC constitutive model gets ignored. I suspect this is a software bug.
With your initial stress you seem to inject a hydrostatic stress state to the entire structure which results in a Mises stress of 0. If the Mises stress is then the yield criteria for the Plastic strain you should consequently expect nothing else than a Plastik strain of 0.
The Mohr-Coulomb and von Mises yield criteria differ fundamentally.Calculation results reveal that simply deleting an irrelevant keyword leads to entirely different computational outputs.
I can reproduce the issue.
It could be related to the use of gapuni elements. The same thing happens if you replace the vertical displacement constraint with a compression-only support in a model that previously worked.
I am not entirely certain about the exact cause, yet my test results demonstrate that the Mohr-Coulomb constitutive model is disregarded and the material behaves as linear elasticity once the *SURFACE INTERACTION keyword is specified, regardless of whether contact elements exist in the model. This occurs even though I have not defined any contact pairs and no contact elements are created accordingly. Therefore, I suspect there may be some interference with the material definition.
plastic
SUMMARY OF C0NVERGENCE INFORMATION
STEP INC ATT ITER CONT. RESID. CORR. RESID. CORR.
EL. FORCE DISP FLUX TEMP.
(#) (%) (%) (%) (%)
1 1 1 1 0 0.1360E+01 0.1000E+03 0.0000E+00 0.0000E+00
1 1 1 2 0 0.3714E+00 0.3288E+01 0.0000E+00 0.0000E+00
1 1 1 3 0 0.4336E-01 0.4085E+00 0.0000E+00 0.0000E+00
elastic
SUMMARY OF C0NVERGENCE INFORMATION
STEP INC ATT ITER CONT. RESID. CORR. RESID. CORR.
EL. FORCE DISP FLUX TEMP.
(#) (%) (%) (%) (%)
1 1 1 1 0 0.0000E+00 0.1000E+03 0.0000E+00 0.0000E+00
1 1 1 2 0 0.0000E+00 0.5228E-10 0.0000E+00 0.0000E+00
I’ve been carrying out debugging work lately, but WSL2 suddenly malfunctioned, drastically slowing down all my operations. 
I have analyzed the calinput.f file and compared the keyword reading sequence for the plastic constitutive model.
For the plastic model, the keywords are read in the order: *material → *plastic → *surfaceinteraction.
In contrast, the Mohr-Coulomb constitutive model follows a different reading sequence: *material → *surfaceinteraction → *mohr-coulomb.
I modified the readinpt.c source file, which appears to resolve the read conflict between the *surfaceinteraction and *mohr-coulomb keywords. The Mohr-Coulomb constitutive model is now properly activated when contact interfaces are defined.
else if(strcmp1(&buff[0],"*PLASTIC")==0){
FORTRAN(keystart,(&ifreeinp,ipoinp,inp,"MATERIAL",
nline,&ikey));
}
else if(strcmp1(&buff[0],"*MOHRCOULOMB")==0){
FORTRAN(keystart,(&ifreeinp,ipoinp,inp,"MATERIAL",
nline,&ikey));
1 Like
Great and Thanks
¿Could you please check if your fix also solves the Mohr-Coulomb+ Only Compression issue?
Regards
Sure, could you share the test file with me? I’ll run a test to verify it.
Updated file: More stable. Ends up with ZERo Plastic Strain ?¿?
*NODE
1,0.125,0,0
2,0,0.125,0
3,0,0,-0.125
4,0,0.125,-0.125
5,0,0,0
6,0.125,0.125,0
7,0.125,0.125,-0.125
8,0.125,0,-0.125
9,0,0,0
10,0,0,-0.09375
11,0,0,-0.0625
12,0,0,-0.03125
13,0,0.03125,-0.125
14,0,0.03125,-0.09375
15,0,0.03125,-0.0625
16,0,0.03125,-0.03125
17,0,0.03125,0
18,0,0.0625,-0.125
19,0,0.0625,-0.09375
20,0,0.0625,-0.0625
21,0,0.0625,-0.03125
22,0,0.0625,0
23,0,0.09375,-0.125
24,0,0.09375,-0.09375
25,0,0.09375,-0.0625
26,0,0.09375,-0.03125
27,0,0.09375,0
28,0,0.125,-0.09375
29,0,0.125,-0.0625
30,0,0.125,-0.03125
31,0.03125,0,-0.125
32,0.03125,0,-0.09375
33,0.03125,0,-0.0625
34,0.03125,0,-0.03125
35,0.03125,0,0
36,0.03125,0.03125,-0.125
37,0.03125,0.03125,-0.09375
38,0.03125,0.03125,-0.0625
39,0.03125,0.03125,-0.03125
40,0.03125,0.03125,0
41,0.03125,0.0625,-0.125
42,0.03125,0.0625,-0.09375
43,0.03125,0.0625,-0.0625
44,0.03125,0.0625,-0.03125
45,0.03125,0.0625,0
46,0.03125,0.09375,-0.125
47,0.03125,0.09375,-0.09375
48,0.03125,0.09375,-0.0625
49,0.03125,0.09375,-0.03125
50,0.03125,0.09375,0
51,0.03125,0.125,-0.125
52,0.03125,0.125,-0.09375
53,0.03125,0.125,-0.0625
54,0.03125,0.125,-0.03125
55,0.03125,0.125,0
56,0.0625,0,-0.125
57,0.0625,0,-0.09375
58,0.0625,0,-0.0625
59,0.0625,0,-0.03125
60,0.0625,0,0
61,0.0625,0.03125,-0.125
62,0.0625,0.03125,-0.09375
63,0.0625,0.03125,-0.0625
64,0.0625,0.03125,-0.03125
65,0.0625,0.03125,0
66,0.0625,0.0625,-0.125
67,0.0625,0.0625,-0.09375
68,0.0625,0.0625,-0.0625
69,0.0625,0.0625,-0.03125
70,0.0625,0.0625,0
71,0.0625,0.09375,-0.125
72,0.0625,0.09375,-0.09375
73,0.0625,0.09375,-0.0625
74,0.0625,0.09375,-0.03125
75,0.0625,0.09375,0
76,0.0625,0.125,-0.125
77,0.0625,0.125,-0.09375
78,0.0625,0.125,-0.0625
79,0.0625,0.125,-0.03125
80,0.0625,0.125,0
81,0.09375,0,-0.125
82,0.09375,0,-0.09375
83,0.09375,0,-0.0625
84,0.09375,0,-0.03125
85,0.09375,0,0
86,0.09375,0.03125,-0.125
87,0.09375,0.03125,-0.09375
88,0.09375,0.03125,-0.0625
89,0.09375,0.03125,-0.03125
90,0.09375,0.03125,0
91,0.09375,0.0625,-0.125
92,0.09375,0.0625,-0.09375
93,0.09375,0.0625,-0.0625
94,0.09375,0.0625,-0.03125
95,0.09375,0.0625,0
96,0.09375,0.09375,-0.125
97,0.09375,0.09375,-0.09375
98,0.09375,0.09375,-0.0625
99,0.09375,0.09375,-0.03125
100,0.09375,0.09375,0
101,0.09375,0.125,-0.125
102,0.09375,0.125,-0.09375
103,0.09375,0.125,-0.0625
104,0.09375,0.125,-0.03125
105,0.09375,0.125,0
106,0.125,0,-0.09375
107,0.125,0,-0.0625
108,0.125,0,-0.03125
109,0.125,0.03125,-0.125
110,0.125,0.03125,-0.09375
111,0.125,0.03125,-0.0625
112,0.125,0.03125,-0.03125
113,0.125,0.03125,0
114,0.125,0.0625,-0.125
115,0.125,0.0625,-0.09375
116,0.125,0.0625,-0.0625
117,0.125,0.0625,-0.03125
118,0.125,0.0625,0
119,0.125,0.09375,-0.125
120,0.125,0.09375,-0.09375
121,0.125,0.09375,-0.0625
122,0.125,0.09375,-0.03125
123,0.125,0.09375,0
124,0.125,0.125,-0.09375
125,0.125,0.125,-0.0625
126,0.125,0.125,-0.03125
127,0,0,-0.125
128,0.03125,0,-0.125
129,0.03125,0,-0.09375
130,0,0,-0.09375
131,0,0,-0.09375
132,0.03125,0,-0.09375
133,0.03125,0,-0.0625
134,0,0,-0.0625
135,0,0,-0.0625
136,0.03125,0,-0.0625
137,0.03125,0,-0.03125
138,0,0,-0.03125
139,0,0,-0.03125
140,0.03125,0,-0.03125
141,0.03125,0,0
142,0,0,0
143,0.03125,0,-0.125
144,0.0625,0,-0.125
145,0.0625,0,-0.09375
146,0.03125,0,-0.09375
147,0.03125,0,-0.09375
148,0.0625,0,-0.09375
149,0.0625,0,-0.0625
150,0.03125,0,-0.0625
151,0.03125,0,-0.0625
152,0.0625,0,-0.0625
153,0.0625,0,-0.03125
154,0.03125,0,-0.03125
155,0.03125,0,-0.03125
156,0.0625,0,-0.03125
157,0.0625,0,0
158,0.03125,0,0
159,0.0625,0,-0.125
160,0.09375,0,-0.125
161,0.09375,0,-0.09375
162,0.0625,0,-0.09375
163,0.0625,0,-0.09375
164,0.09375,0,-0.09375
165,0.09375,0,-0.0625
166,0.0625,0,-0.0625
167,0.0625,0,-0.0625
168,0.09375,0,-0.0625
169,0.09375,0,-0.03125
170,0.0625,0,-0.03125
171,0.0625,0,-0.03125
172,0.09375,0,-0.03125
173,0.09375,0,0
174,0.0625,0,0
175,0.09375,0,-0.125
176,0.125,0,-0.125
177,0.125,0,-0.09375
178,0.09375,0,-0.09375
179,0.09375,0,-0.09375
180,0.125,0,-0.09375
181,0.125,0,-0.0625
182,0.09375,0,-0.0625
183,0.09375,0,-0.0625
184,0.125,0,-0.0625
185,0.125,0,-0.03125
186,0.09375,0,-0.03125
187,0.09375,0,-0.03125
188,0.125,0,-0.03125
189,0.125,0,0
190,0.09375,0,0
*ELEMENT,TYPE=C3D8I
1,99,122,126,104,100,123,6,105
2,3,31,36,13,10,32,37,14
3,10,32,37,14,11,33,38,15
4,11,33,38,15,12,34,39,16
5,12,34,39,16,5,35,40,17
6,13,36,41,18,14,37,42,19
7,14,37,42,19,15,38,43,20
8,15,38,43,20,16,39,44,21
9,16,39,44,21,17,40,45,22
10,18,41,46,23,19,42,47,24
11,19,42,47,24,20,43,48,25
12,20,43,48,25,21,44,49,26
13,21,44,49,26,22,45,50,27
14,23,46,51,4,24,47,52,28
15,24,47,52,28,25,48,53,29
16,25,48,53,29,26,49,54,30
17,26,49,54,30,27,50,55,2
18,31,56,61,36,32,57,62,37
19,32,57,62,37,33,58,63,38
20,33,58,63,38,34,59,64,39
21,34,59,64,39,35,60,65,40
22,36,61,66,41,37,62,67,42
23,37,62,67,42,38,63,68,43
24,38,63,68,43,39,64,69,44
25,39,64,69,44,40,65,70,45
26,41,66,71,46,42,67,72,47
27,42,67,72,47,43,68,73,48
28,43,68,73,48,44,69,74,49
29,44,69,74,49,45,70,75,50
30,46,71,76,51,47,72,77,52
31,47,72,77,52,48,73,78,53
32,48,73,78,53,49,74,79,54
33,49,74,79,54,50,75,80,55
34,56,81,86,61,57,82,87,62
35,57,82,87,62,58,83,88,63
36,58,83,88,63,59,84,89,64
37,59,84,89,64,60,85,90,65
38,61,86,91,66,62,87,92,67
39,62,87,92,67,63,88,93,68
40,63,88,93,68,64,89,94,69
41,64,89,94,69,65,90,95,70
42,66,91,96,71,67,92,97,72
43,67,92,97,72,68,93,98,73
44,68,93,98,73,69,94,99,74
45,69,94,99,74,70,95,100,75
46,71,96,101,76,72,97,102,77
47,72,97,102,77,73,98,103,78
48,73,98,103,78,74,99,104,79
49,74,99,104,79,75,100,105,80
50,81,8,109,86,82,106,110,87
51,82,106,110,87,83,107,111,88
52,83,107,111,88,84,108,112,89
53,84,108,112,89,85,1,113,90
54,86,109,114,91,87,110,115,92
55,87,110,115,92,88,111,116,93
56,88,111,116,93,89,112,117,94
57,89,112,117,94,90,113,118,95
58,91,114,119,96,92,115,120,97
59,92,115,120,97,93,116,121,98
60,93,116,121,98,94,117,122,99
61,94,117,122,99,95,118,123,100
62,96,119,7,101,97,120,124,102
63,97,120,124,102,98,121,125,103
64,98,121,125,103,99,122,126,104
*ELEMENT,TYPE=GAPUNI
65,3,127
66,31,128
67,32,129
68,10,130
69,10,131
70,32,132
71,33,133
72,11,134
73,11,135
74,33,136
75,34,137
76,12,138
77,12,139
78,34,140
79,35,141
80,5,142
81,31,143
82,56,144
83,57,145
84,32,146
85,32,147
86,57,148
87,58,149
88,33,150
89,33,151
90,58,152
91,59,153
92,34,154
93,34,155
94,59,156
95,60,157
96,35,158
97,56,159
98,81,160
99,82,161
100,57,162
101,57,163
102,82,164
103,83,165
104,58,166
105,58,167
106,83,168
107,84,169
108,59,170
109,59,171
110,84,172
111,85,173
112,60,174
113,81,175
114,8,176
115,106,177
116,82,178
117,82,179
118,106,180
119,107,181
120,83,182
121,83,183
122,107,184
123,108,185
124,84,186
125,84,187
126,108,188
127,1,189
128,85,190
*NSET,NSET=COMPRESSIONONLY
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
*NSET,NSET=TEMPERATURE_ELEMENTS
1
2
3
4
5
6
7
8
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
*NSET,NSET=NALL
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
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
*ELSET,ELSET=MOHR_TEST
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
*ELSET,ELSET=COMPRESSIONONLY_1
65
*ELSET,ELSET=COMPRESSIONONLY_2
66
*ELSET,ELSET=COMPRESSIONONLY_3
67
*ELSET,ELSET=COMPRESSIONONLY_4
68
*ELSET,ELSET=COMPRESSIONONLY_5
69
*ELSET,ELSET=COMPRESSIONONLY_6
70
*ELSET,ELSET=COMPRESSIONONLY_7
71
*ELSET,ELSET=COMPRESSIONONLY_8
72
*ELSET,ELSET=COMPRESSIONONLY_9
73
*ELSET,ELSET=COMPRESSIONONLY_A
74
*ELSET,ELSET=COMPRESSIONONLY_B
75
*ELSET,ELSET=COMPRESSIONONLY_C
76
*ELSET,ELSET=COMPRESSIONONLY_D
77
*ELSET,ELSET=COMPRESSIONONLY_E
78
*ELSET,ELSET=COMPRESSIONONLY_F
79
*ELSET,ELSET=COMPRESSIONONLY_10
80
*ELSET,ELSET=COMPRESSIONONLY_11
81
*ELSET,ELSET=COMPRESSIONONLY_12
82
*ELSET,ELSET=COMPRESSIONONLY_13
83
*ELSET,ELSET=COMPRESSIONONLY_14
84
*ELSET,ELSET=COMPRESSIONONLY_15
85
*ELSET,ELSET=COMPRESSIONONLY_16
86
*ELSET,ELSET=COMPRESSIONONLY_17
87
*ELSET,ELSET=COMPRESSIONONLY_18
88
*ELSET,ELSET=COMPRESSIONONLY_19
89
*ELSET,ELSET=COMPRESSIONONLY_1A
90
*ELSET,ELSET=COMPRESSIONONLY_1B
91
*ELSET,ELSET=COMPRESSIONONLY_1C
92
*ELSET,ELSET=COMPRESSIONONLY_1D
93
*ELSET,ELSET=COMPRESSIONONLY_1E
94
*ELSET,ELSET=COMPRESSIONONLY_1F
95
*ELSET,ELSET=COMPRESSIONONLY_20
96
*ELSET,ELSET=COMPRESSIONONLY_21
97
*ELSET,ELSET=COMPRESSIONONLY_22
98
*ELSET,ELSET=COMPRESSIONONLY_23
99
*ELSET,ELSET=COMPRESSIONONLY_24
100
*ELSET,ELSET=COMPRESSIONONLY_25
101
*ELSET,ELSET=COMPRESSIONONLY_26
102
*ELSET,ELSET=COMPRESSIONONLY_27
103
*ELSET,ELSET=COMPRESSIONONLY_28
104
*ELSET,ELSET=COMPRESSIONONLY_29
105
*ELSET,ELSET=COMPRESSIONONLY_2A
106
*ELSET,ELSET=COMPRESSIONONLY_2B
107
*ELSET,ELSET=COMPRESSIONONLY_2C
108
*ELSET,ELSET=COMPRESSIONONLY_2D
109
*ELSET,ELSET=COMPRESSIONONLY_2E
110
*ELSET,ELSET=COMPRESSIONONLY_2F
111
*ELSET,ELSET=COMPRESSIONONLY_30
112
*ELSET,ELSET=COMPRESSIONONLY_31
113
*ELSET,ELSET=COMPRESSIONONLY_32
114
*ELSET,ELSET=COMPRESSIONONLY_33
115
*ELSET,ELSET=COMPRESSIONONLY_34
116
*ELSET,ELSET=COMPRESSIONONLY_35
117
*ELSET,ELSET=COMPRESSIONONLY_36
118
*ELSET,ELSET=COMPRESSIONONLY_37
119
*ELSET,ELSET=COMPRESSIONONLY_38
120
*ELSET,ELSET=COMPRESSIONONLY_39
121
*ELSET,ELSET=COMPRESSIONONLY_3A
122
*ELSET,ELSET=COMPRESSIONONLY_3B
123
*ELSET,ELSET=COMPRESSIONONLY_3C
124
*ELSET,ELSET=COMPRESSIONONLY_3D
125
*ELSET,ELSET=COMPRESSIONONLY_3E
126
*ELSET,ELSET=COMPRESSIONONLY_3F
127
*ELSET,ELSET=COMPRESSIONONLY_40
128
*SURFACE,NAME=REACTION
2,S3
3,S3
4,S3
5,S3
18,S3
19,S3
20,S3
21,S3
34,S3
35,S3
36,S3
37,S3
50,S3
51,S3
52,S3
53,S3
*SURFACE,NAME=MASTER
14,S5
15,S5
16,S5
17,S5
30,S5
31,S5
32,S5
33,S5
46,S5
47,S5
48,S5
49,S5
62,S5
63,S5
64,S5
1,S5
*SURFACE,NAME=SIDE_FACES
6,S1
2,S1
10,S1
14,S1
18,S1
22,S1
26,S1
30,S1
34,S1
38,S1
42,S1
46,S1
50,S1
54,S1
58,S1
62,S1
5,S2
9,S2
13,S2
17,S2
21,S2
25,S2
29,S2
33,S2
37,S2
41,S2
45,S2
49,S2
53,S2
57,S2
61,S2
1,S2
*BOUNDARY
COMPRESSIONONLY,1,,0
COMPRESSIONONLY,2,,0
COMPRESSIONONLY,3,,0
*BOUNDARY
3,1,,0
5,1,,0
5,3,,0
10,1,,0
11,1,,0
12,1,,0
*GAP,ELSET=COMPRESSIONONLY_1
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_4
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_5
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_6
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_7
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_8
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_9
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_A
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_B
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_C
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_D
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_E
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_F
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_10
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_11
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_12
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_13
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_14
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_15
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_16
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_17
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_18
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_19
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_1A
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_1B
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_1C
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_1D
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_1E
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_1F
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_20
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_21
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_22
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_23
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_24
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_25
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_26
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_27
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_28
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_29
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2A
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2B
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2C
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2D
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2E
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_2F
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_30
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_31
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_32
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_33
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_34
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_35
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_36
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_37
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_38
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_39
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3A
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3B
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3C
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3D
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3E
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_3F
0,0,-1,0,,5126953125,2441.40625
*GAP,ELSET=COMPRESSIONONLY_40
0,0,-1,0,,5126953125,2441.40625
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,300
*AMPLITUDE,NAME=A_41
0,0
1,50000000
*AMPLITUDE,NAME=A_42
0,0
1,-10000000
*MATERIAL,NAME=MOHR-MB
*ELASTIC
210000.E6,0.26,300.
*DENSITY
7850
*MOHR COULOMB
35.,35.,300.
*MOHR COULOMB HARDENING
2.5E6,0.001,300.
60.0E6,0.1,300.
*SOLID SECTION,ELSET=Mohr_Test,material=MOHR-MB
*STEP,NLGEOM=YES,INC=200,AMPLITUDE=STEP
*STATIC
0.01,1,0,0.01
*TEMPERATURE
TEMPERATURE_ELEMENTS,300
*TEMPERATURE
9,300
127,300
128,300
129,300
130,300
131,300
132,300
133,300
134,300
135,300
136,300
137,300
138,300
139,300
140,300
141,300
142,300
143,300
144,300
145,300
146,300
147,300
148,300
149,300
150,300
151,300
152,300
153,300
154,300
155,300
156,300
157,300
158,300
159,300
160,300
161,300
162,300
163,300
164,300
165,300
166,300
167,300
168,300
169,300
170,300
171,300
172,300
173,300
174,300
175,300
176,300
177,300
178,300
179,300
180,300
181,300
182,300
183,300
184,300
185,300
186,300
187,300
188,300
189,300
190,300
*DSLOAD,AMPLITUDE=A_41
MASTER,P,1
*DSLOAD,AMPLITUDE=A_42
SIDE_FACES,P,1
*NODE FILE,GLOBAL=YES
U,RF
*EL FILE
S,NOE,E,PEEQ,ENER
*END STEP
My results are the same.
After some consideration, I added a few extra revisions to readinput.c as a precaution.
else if(strcmp1(&buff[0],"*PLASTIC")==0){
FORTRAN(keystart,(&ifreeinp,ipoinp,inp,"MATERIAL",
nline,&ikey));
}
else if(strcmp1(&buff[0],"*MOHRCOULOMB")==0){
FORTRAN(keystart,(&ifreeinp,ipoinp,inp,"MATERIAL",
nline,&ikey));
}
else if(strcmp1(&buff[0],"*MOHRCOULOMBHARDENING")==0){
FORTRAN(keystart,(&ifreeinp,ipoinp,inp,"MATERIAL",
nline,&ikey));
}
Thank you very much 