Benoit44
Mechanical
- Jan 13, 2014
- 1
Dear All,
I would like to compute an anisotropic linear elastic law for composite.
I wish to output stress and strain in order to calculate 2D and 3D failure criteria (Hashin 3D, Puck...)
I am a beginner in UMAT computing.
When I run a calculation without UMAT, I do not find the same thing when I do the calculation with the UMAT I computed.
For instance, please find below my the .inp without the UMAT:
*Heading
** Job name: sans_umat_90degre_2D Model name: sans_umat_orthotrope_2D
** Generated by: Abaqus/CAE 6.12-3
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
1, 55., 25., 0.
2, 51., 25., 0.
3, 47., 25., 0.
4, 43., 25., 0.
5, 39., 25., 0.
6, 35., 25., 0.
7, 31., 25., 0.
8, 27., 25., 0.
9, 23., 25., 0.
10, 19., 25., 0.
11, 15., 25., 0.
12, 11., 25., 0.
13, 7., 25., 0.
14, 3., 25., 0.
15, -1., 25., 0.
16, -5., 25., 0.
17, -9., 25., 0.
18, -13., 25., 0.
19, -17., 25., 0.
20, -21., 25., 0.
21, -25., 25., 0.
22, -29., 25., 0.
23, -33., 25., 0.
24, -37., 25., 0.
25, -41., 25., 0.
26, -45., 25., 0.
27, 55., 20.833334, 0.
28, 51., 20.833334, 0.
29, 47., 20.833334, 0.
30, 43., 20.833334, 0.
31, 39., 20.833334, 0.
32, 35., 20.833334, 0.
33, 31., 20.833334, 0.
34, 27., 20.833334, 0.
35, 23., 20.833334, 0.
36, 19., 20.833334, 0.
37, 15., 20.833334, 0.
38, 11., 20.833334, 0.
39, 7., 20.833334, 0.
40, 3., 20.833334, 0.
41, -1., 20.833334, 0.
42, -5., 20.833334, 0.
43, -9., 20.833334, 0.
44, -13., 20.833334, 0.
45, -17., 20.833334, 0.
46, -21., 20.833334, 0.
47, -25., 20.833334, 0.
48, -29., 20.833334, 0.
49, -33., 20.833334, 0.
50, -37., 20.833334, 0.
51, -41., 20.833334, 0.
52, -45., 20.833334, 0.
53, 55., 16.666666, 0.
54, 51., 16.666666, 0.
55, 47., 16.666666, 0.
56, 43., 16.666666, 0.
57, 39., 16.666666, 0.
58, 35., 16.666666, 0.
59, 31., 16.666666, 0.
60, 27., 16.666666, 0.
61, 23., 16.666666, 0.
62, 19., 16.666666, 0.
63, 15., 16.666666, 0.
64, 11., 16.666666, 0.
65, 7., 16.666666, 0.
66, 3., 16.666666, 0.
67, -1., 16.666666, 0.
68, -5., 16.666666, 0.
69, -9., 16.666666, 0.
70, -13., 16.666666, 0.
71, -17., 16.666666, 0.
72, -21., 16.666666, 0.
73, -25., 16.666666, 0.
74, -29., 16.666666, 0.
75, -33., 16.666666, 0.
76, -37., 16.666666, 0.
77, -41., 16.666666, 0.
78, -45., 16.666666, 0.
79, 55., 12.5, 0.
80, 51., 12.5, 0.
81, 47., 12.5, 0.
82, 43., 12.5, 0.
83, 39., 12.5, 0.
84, 35., 12.5, 0.
85, 31., 12.5, 0.
86, 27., 12.5, 0.
87, 23., 12.5, 0.
88, 19., 12.5, 0.
89, 15., 12.5, 0.
90, 11., 12.5, 0.
91, 7., 12.5, 0.
92, 3., 12.5, 0.
93, -1., 12.5, 0.
94, -5., 12.5, 0.
95, -9., 12.5, 0.
96, -13., 12.5, 0.
97, -17., 12.5, 0.
98, -21., 12.5, 0.
99, -25., 12.5, 0.
100, -29., 12.5, 0.
101, -33., 12.5, 0.
102, -37., 12.5, 0.
103, -41., 12.5, 0.
104, -45., 12.5, 0.
105, 55., 8.33333302, 0.
106, 51., 8.33333302, 0.
107, 47., 8.33333302, 0.
108, 43., 8.33333302, 0.
109, 39., 8.33333302, 0.
110, 35., 8.33333302, 0.
111, 31., 8.33333302, 0.
112, 27., 8.33333302, 0.
113, 23., 8.33333302, 0.
114, 19., 8.33333302, 0.
115, 15., 8.33333302, 0.
116, 11., 8.33333302, 0.
117, 7., 8.33333302, 0.
118, 3., 8.33333302, 0.
119, -1., 8.33333302, 0.
120, -5., 8.33333302, 0.
121, -9., 8.33333302, 0.
122, -13., 8.33333302, 0.
123, -17., 8.33333302, 0.
124, -21., 8.33333302, 0.
125, -25., 8.33333302, 0.
126, -29., 8.33333302, 0.
127, -33., 8.33333302, 0.
128, -37., 8.33333302, 0.
129, -41., 8.33333302, 0.
130, -45., 8.33333302, 0.
131, 55., 4.16666651, 0.
132, 51., 4.16666651, 0.
133, 47., 4.16666651, 0.
134, 43., 4.16666651, 0.
135, 39., 4.16666651, 0.
136, 35., 4.16666651, 0.
137, 31., 4.16666651, 0.
138, 27., 4.16666651, 0.
139, 23., 4.16666651, 0.
140, 19., 4.16666651, 0.
141, 15., 4.16666651, 0.
142, 11., 4.16666651, 0.
143, 7., 4.16666651, 0.
144, 3., 4.16666651, 0.
145, -1., 4.16666651, 0.
146, -5., 4.16666651, 0.
147, -9., 4.16666651, 0.
148, -13., 4.16666651, 0.
149, -17., 4.16666651, 0.
150, -21., 4.16666651, 0.
151, -25., 4.16666651, 0.
152, -29., 4.16666651, 0.
153, -33., 4.16666651, 0.
154, -37., 4.16666651, 0.
155, -41., 4.16666651, 0.
156, -45., 4.16666651, 0.
157, 55., 0., 0.
158, 51., 0., 0.
159, 47., 0., 0.
160, 43., 0., 0.
161, 39., 0., 0.
162, 35., 0., 0.
163, 31., 0., 0.
164, 27., 0., 0.
165, 23., 0., 0.
166, 19., 0., 0.
167, 15., 0., 0.
168, 11., 0., 0.
169, 7., 0., 0.
170, 3., 0., 0.
171, -1., 0., 0.
172, -5., 0., 0.
173, -9., 0., 0.
174, -13., 0., 0.
175, -17., 0., 0.
176, -21., 0., 0.
177, -25., 0., 0.
178, -29., 0., 0.
179, -33., 0., 0.
180, -37., 0., 0.
181, -41., 0., 0.
182, -45., 0., 0.
*Element, type=S4R
1, 1, 2, 28, 27
2, 2, 3, 29, 28
3, 3, 4, 30, 29
4, 4, 5, 31, 30
5, 5, 6, 32, 31
6, 6, 7, 33, 32
7, 7, 8, 34, 33
8, 8, 9, 35, 34
9, 9, 10, 36, 35
10, 10, 11, 37, 36
11, 11, 12, 38, 37
12, 12, 13, 39, 38
13, 13, 14, 40, 39
14, 14, 15, 41, 40
15, 15, 16, 42, 41
16, 16, 17, 43, 42
17, 17, 18, 44, 43
18, 18, 19, 45, 44
19, 19, 20, 46, 45
20, 20, 21, 47, 46
21, 21, 22, 48, 47
22, 22, 23, 49, 48
23, 23, 24, 50, 49
24, 24, 25, 51, 50
25, 25, 26, 52, 51
26, 27, 28, 54, 53
27, 28, 29, 55, 54
28, 29, 30, 56, 55
29, 30, 31, 57, 56
30, 31, 32, 58, 57
31, 32, 33, 59, 58
32, 33, 34, 60, 59
33, 34, 35, 61, 60
34, 35, 36, 62, 61
35, 36, 37, 63, 62
36, 37, 38, 64, 63
37, 38, 39, 65, 64
38, 39, 40, 66, 65
39, 40, 41, 67, 66
40, 41, 42, 68, 67
41, 42, 43, 69, 68
42, 43, 44, 70, 69
43, 44, 45, 71, 70
44, 45, 46, 72, 71
45, 46, 47, 73, 72
46, 47, 48, 74, 73
47, 48, 49, 75, 74
48, 49, 50, 76, 75
49, 50, 51, 77, 76
50, 51, 52, 78, 77
51, 53, 54, 80, 79
52, 54, 55, 81, 80
53, 55, 56, 82, 81
54, 56, 57, 83, 82
55, 57, 58, 84, 83
56, 58, 59, 85, 84
57, 59, 60, 86, 85
58, 60, 61, 87, 86
59, 61, 62, 88, 87
60, 62, 63, 89, 88
61, 63, 64, 90, 89
62, 64, 65, 91, 90
63, 65, 66, 92, 91
64, 66, 67, 93, 92
65, 67, 68, 94, 93
66, 68, 69, 95, 94
67, 69, 70, 96, 95
68, 70, 71, 97, 96
69, 71, 72, 98, 97
70, 72, 73, 99, 98
71, 73, 74, 100, 99
72, 74, 75, 101, 100
73, 75, 76, 102, 101
74, 76, 77, 103, 102
75, 77, 78, 104, 103
76, 79, 80, 106, 105
77, 80, 81, 107, 106
78, 81, 82, 108, 107
79, 82, 83, 109, 108
80, 83, 84, 110, 109
81, 84, 85, 111, 110
82, 85, 86, 112, 111
83, 86, 87, 113, 112
84, 87, 88, 114, 113
85, 88, 89, 115, 114
86, 89, 90, 116, 115
87, 90, 91, 117, 116
88, 91, 92, 118, 117
89, 92, 93, 119, 118
90, 93, 94, 120, 119
91, 94, 95, 121, 120
92, 95, 96, 122, 121
93, 96, 97, 123, 122
94, 97, 98, 124, 123
95, 98, 99, 125, 124
96, 99, 100, 126, 125
97, 100, 101, 127, 126
98, 101, 102, 128, 127
99, 102, 103, 129, 128
100, 103, 104, 130, 129
101, 105, 106, 132, 131
102, 106, 107, 133, 132
103, 107, 108, 134, 133
104, 108, 109, 135, 134
105, 109, 110, 136, 135
106, 110, 111, 137, 136
107, 111, 112, 138, 137
108, 112, 113, 139, 138
109, 113, 114, 140, 139
110, 114, 115, 141, 140
111, 115, 116, 142, 141
112, 116, 117, 143, 142
113, 117, 118, 144, 143
114, 118, 119, 145, 144
115, 119, 120, 146, 145
116, 120, 121, 147, 146
117, 121, 122, 148, 147
118, 122, 123, 149, 148
119, 123, 124, 150, 149
120, 124, 125, 151, 150
121, 125, 126, 152, 151
122, 126, 127, 153, 152
123, 127, 128, 154, 153
124, 128, 129, 155, 154
125, 129, 130, 156, 155
126, 131, 132, 158, 157
127, 132, 133, 159, 158
128, 133, 134, 160, 159
129, 134, 135, 161, 160
130, 135, 136, 162, 161
131, 136, 137, 163, 162
132, 137, 138, 164, 163
133, 138, 139, 165, 164
134, 139, 140, 166, 165
135, 140, 141, 167, 166
136, 141, 142, 168, 167
137, 142, 143, 169, 168
138, 143, 144, 170, 169
139, 144, 145, 171, 170
140, 145, 146, 172, 171
141, 146, 147, 173, 172
142, 147, 148, 174, 173
143, 148, 149, 175, 174
144, 149, 150, 176, 175
145, 150, 151, 177, 176
146, 151, 152, 178, 177
147, 152, 153, 179, 178
148, 153, 154, 180, 179
149, 154, 155, 181, 180
150, 155, 156, 182, 181
*Nset, nset=Set-1, generate
1, 182, 1
*Elset, elset=Set-1, generate
1, 150, 1
*Nset, nset=_PickedSet5, internal, generate
1, 182, 1
*Elset, elset=_PickedSet5, internal, generate
1, 150, 1
** Region: (CompositeLayup-1-1: Generated From Layup)
*Elset, elset=CompositeLayup-1-1, generate
1, 150, 1
** Section: CompositeLayup-1-1
*Shell Section, elset=CompositeLayup-1-1, composite, layup=CompositeLayup-1
1., 3, Carbone, 90., Ply-1
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*End Instance
**
*Nset, nset=Set-1, instance=Part-1-1
182,
*Nset, nset=Set-2, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-2, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-3, instance=Part-1-1, generate
157, 182, 1
*Elset, elset=Set-3, instance=Part-1-1, generate
126, 150, 1
*Nset, nset=Set-4, instance=Part-1-1, generate
1, 26, 1
*Elset, elset=Set-4, instance=Part-1-1, generate
1, 25, 1
*Nset, nset=Set-5, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-5, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-6, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-6, instance=Part-1-1, generate
1, 126, 25
*Nset, nset=Set-7, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-7, instance=Part-1-1, generate
1, 126, 25
*Elset, elset=_Surf-1_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_E4, E4
*Elset, elset=_Surf-2_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-2
_Surf-2_E4, E4
*Elset, elset=_Surf-3_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-3
_Surf-3_E4, E4
*Elset, elset=_Surf-4_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-4
_Surf-4_E4, E4
*End Assembly
**
** MATERIALS
**
*Material, name=Carbone
*Damage Initiation, criterion=HASHIN
2500.,1000., 50., 65., 30., 30.
*Elastic, type=ENGINEERING CONSTANTS
130000.,9000.,9000., 0.36, 0.3, 0.3,3600.,3600.
3100.,
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES
*Static
1., 1., 1e-05, 1.
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
Set-1, ENCASTRE
** Name: BC-2 Type: Displacement/Rotation
*Boundary
Set-2, 1, 1
Set-2, 2, 2
Set-2, 3, 3
Set-2, 4, 4
Set-2, 5, 5
Set-2, 6, 6
**
** LOADS
**
** Name: Load-1 Type: Shell edge load
*Dsload
Surf-2, EDNOR, -10.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
U,
*Element Output, directions=YES
E, LE, NE, S
**
** FIELD OUTPUT: F-Output-2
**
*Element Output, elset=Part-1-1.CompositeLayup-1-1, directions=YES
1, 2, 3
E, LE, NE, S
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history
*Energy Output
ALLAE, ALLCD, ALLDMD, ALLEE, ALLFD, ALLIE, ALLJD, ALLKE, ALLKL, ALLPD, ALLQB, ALLSD, ALLSE, ALLVD, ALLWK, ETOTAL
*Contact Output
BDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC
*End Step
And the .inp with the UMAT (you have to copy/paste in a .inp file to run in Abaqus and launch the calculation):
*Heading
** Job name: avec_umat_90degre_2D Model name: avec_umat_orthotrope_2D
** Generated by: Abaqus/CAE 6.12-3
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
1, 55., 25., 0.
2, 51., 25., 0.
3, 47., 25., 0.
4, 43., 25., 0.
5, 39., 25., 0.
6, 35., 25., 0.
7, 31., 25., 0.
8, 27., 25., 0.
9, 23., 25., 0.
10, 19., 25., 0.
11, 15., 25., 0.
12, 11., 25., 0.
13, 7., 25., 0.
14, 3., 25., 0.
15, -1., 25., 0.
16, -5., 25., 0.
17, -9., 25., 0.
18, -13., 25., 0.
19, -17., 25., 0.
20, -21., 25., 0.
21, -25., 25., 0.
22, -29., 25., 0.
23, -33., 25., 0.
24, -37., 25., 0.
25, -41., 25., 0.
26, -45., 25., 0.
27, 55., 20.833334, 0.
28, 51., 20.833334, 0.
29, 47., 20.833334, 0.
30, 43., 20.833334, 0.
31, 39., 20.833334, 0.
32, 35., 20.833334, 0.
33, 31., 20.833334, 0.
34, 27., 20.833334, 0.
35, 23., 20.833334, 0.
36, 19., 20.833334, 0.
37, 15., 20.833334, 0.
38, 11., 20.833334, 0.
39, 7., 20.833334, 0.
40, 3., 20.833334, 0.
41, -1., 20.833334, 0.
42, -5., 20.833334, 0.
43, -9., 20.833334, 0.
44, -13., 20.833334, 0.
45, -17., 20.833334, 0.
46, -21., 20.833334, 0.
47, -25., 20.833334, 0.
48, -29., 20.833334, 0.
49, -33., 20.833334, 0.
50, -37., 20.833334, 0.
51, -41., 20.833334, 0.
52, -45., 20.833334, 0.
53, 55., 16.666666, 0.
54, 51., 16.666666, 0.
55, 47., 16.666666, 0.
56, 43., 16.666666, 0.
57, 39., 16.666666, 0.
58, 35., 16.666666, 0.
59, 31., 16.666666, 0.
60, 27., 16.666666, 0.
61, 23., 16.666666, 0.
62, 19., 16.666666, 0.
63, 15., 16.666666, 0.
64, 11., 16.666666, 0.
65, 7., 16.666666, 0.
66, 3., 16.666666, 0.
67, -1., 16.666666, 0.
68, -5., 16.666666, 0.
69, -9., 16.666666, 0.
70, -13., 16.666666, 0.
71, -17., 16.666666, 0.
72, -21., 16.666666, 0.
73, -25., 16.666666, 0.
74, -29., 16.666666, 0.
75, -33., 16.666666, 0.
76, -37., 16.666666, 0.
77, -41., 16.666666, 0.
78, -45., 16.666666, 0.
79, 55., 12.5, 0.
80, 51., 12.5, 0.
81, 47., 12.5, 0.
82, 43., 12.5, 0.
83, 39., 12.5, 0.
84, 35., 12.5, 0.
85, 31., 12.5, 0.
86, 27., 12.5, 0.
87, 23., 12.5, 0.
88, 19., 12.5, 0.
89, 15., 12.5, 0.
90, 11., 12.5, 0.
91, 7., 12.5, 0.
92, 3., 12.5, 0.
93, -1., 12.5, 0.
94, -5., 12.5, 0.
95, -9., 12.5, 0.
96, -13., 12.5, 0.
97, -17., 12.5, 0.
98, -21., 12.5, 0.
99, -25., 12.5, 0.
100, -29., 12.5, 0.
101, -33., 12.5, 0.
102, -37., 12.5, 0.
103, -41., 12.5, 0.
104, -45., 12.5, 0.
105, 55., 8.33333302, 0.
106, 51., 8.33333302, 0.
107, 47., 8.33333302, 0.
108, 43., 8.33333302, 0.
109, 39., 8.33333302, 0.
110, 35., 8.33333302, 0.
111, 31., 8.33333302, 0.
112, 27., 8.33333302, 0.
113, 23., 8.33333302, 0.
114, 19., 8.33333302, 0.
115, 15., 8.33333302, 0.
116, 11., 8.33333302, 0.
117, 7., 8.33333302, 0.
118, 3., 8.33333302, 0.
119, -1., 8.33333302, 0.
120, -5., 8.33333302, 0.
121, -9., 8.33333302, 0.
122, -13., 8.33333302, 0.
123, -17., 8.33333302, 0.
124, -21., 8.33333302, 0.
125, -25., 8.33333302, 0.
126, -29., 8.33333302, 0.
127, -33., 8.33333302, 0.
128, -37., 8.33333302, 0.
129, -41., 8.33333302, 0.
130, -45., 8.33333302, 0.
131, 55., 4.16666651, 0.
132, 51., 4.16666651, 0.
133, 47., 4.16666651, 0.
134, 43., 4.16666651, 0.
135, 39., 4.16666651, 0.
136, 35., 4.16666651, 0.
137, 31., 4.16666651, 0.
138, 27., 4.16666651, 0.
139, 23., 4.16666651, 0.
140, 19., 4.16666651, 0.
141, 15., 4.16666651, 0.
142, 11., 4.16666651, 0.
143, 7., 4.16666651, 0.
144, 3., 4.16666651, 0.
145, -1., 4.16666651, 0.
146, -5., 4.16666651, 0.
147, -9., 4.16666651, 0.
148, -13., 4.16666651, 0.
149, -17., 4.16666651, 0.
150, -21., 4.16666651, 0.
151, -25., 4.16666651, 0.
152, -29., 4.16666651, 0.
153, -33., 4.16666651, 0.
154, -37., 4.16666651, 0.
155, -41., 4.16666651, 0.
156, -45., 4.16666651, 0.
157, 55., 0., 0.
158, 51., 0., 0.
159, 47., 0., 0.
160, 43., 0., 0.
161, 39., 0., 0.
162, 35., 0., 0.
163, 31., 0., 0.
164, 27., 0., 0.
165, 23., 0., 0.
166, 19., 0., 0.
167, 15., 0., 0.
168, 11., 0., 0.
169, 7., 0., 0.
170, 3., 0., 0.
171, -1., 0., 0.
172, -5., 0., 0.
173, -9., 0., 0.
174, -13., 0., 0.
175, -17., 0., 0.
176, -21., 0., 0.
177, -25., 0., 0.
178, -29., 0., 0.
179, -33., 0., 0.
180, -37., 0., 0.
181, -41., 0., 0.
182, -45., 0., 0.
*Element, type=S4R
1, 1, 2, 28, 27
2, 2, 3, 29, 28
3, 3, 4, 30, 29
4, 4, 5, 31, 30
5, 5, 6, 32, 31
6, 6, 7, 33, 32
7, 7, 8, 34, 33
8, 8, 9, 35, 34
9, 9, 10, 36, 35
10, 10, 11, 37, 36
11, 11, 12, 38, 37
12, 12, 13, 39, 38
13, 13, 14, 40, 39
14, 14, 15, 41, 40
15, 15, 16, 42, 41
16, 16, 17, 43, 42
17, 17, 18, 44, 43
18, 18, 19, 45, 44
19, 19, 20, 46, 45
20, 20, 21, 47, 46
21, 21, 22, 48, 47
22, 22, 23, 49, 48
23, 23, 24, 50, 49
24, 24, 25, 51, 50
25, 25, 26, 52, 51
26, 27, 28, 54, 53
27, 28, 29, 55, 54
28, 29, 30, 56, 55
29, 30, 31, 57, 56
30, 31, 32, 58, 57
31, 32, 33, 59, 58
32, 33, 34, 60, 59
33, 34, 35, 61, 60
34, 35, 36, 62, 61
35, 36, 37, 63, 62
36, 37, 38, 64, 63
37, 38, 39, 65, 64
38, 39, 40, 66, 65
39, 40, 41, 67, 66
40, 41, 42, 68, 67
41, 42, 43, 69, 68
42, 43, 44, 70, 69
43, 44, 45, 71, 70
44, 45, 46, 72, 71
45, 46, 47, 73, 72
46, 47, 48, 74, 73
47, 48, 49, 75, 74
48, 49, 50, 76, 75
49, 50, 51, 77, 76
50, 51, 52, 78, 77
51, 53, 54, 80, 79
52, 54, 55, 81, 80
53, 55, 56, 82, 81
54, 56, 57, 83, 82
55, 57, 58, 84, 83
56, 58, 59, 85, 84
57, 59, 60, 86, 85
58, 60, 61, 87, 86
59, 61, 62, 88, 87
60, 62, 63, 89, 88
61, 63, 64, 90, 89
62, 64, 65, 91, 90
63, 65, 66, 92, 91
64, 66, 67, 93, 92
65, 67, 68, 94, 93
66, 68, 69, 95, 94
67, 69, 70, 96, 95
68, 70, 71, 97, 96
69, 71, 72, 98, 97
70, 72, 73, 99, 98
71, 73, 74, 100, 99
72, 74, 75, 101, 100
73, 75, 76, 102, 101
74, 76, 77, 103, 102
75, 77, 78, 104, 103
76, 79, 80, 106, 105
77, 80, 81, 107, 106
78, 81, 82, 108, 107
79, 82, 83, 109, 108
80, 83, 84, 110, 109
81, 84, 85, 111, 110
82, 85, 86, 112, 111
83, 86, 87, 113, 112
84, 87, 88, 114, 113
85, 88, 89, 115, 114
86, 89, 90, 116, 115
87, 90, 91, 117, 116
88, 91, 92, 118, 117
89, 92, 93, 119, 118
90, 93, 94, 120, 119
91, 94, 95, 121, 120
92, 95, 96, 122, 121
93, 96, 97, 123, 122
94, 97, 98, 124, 123
95, 98, 99, 125, 124
96, 99, 100, 126, 125
97, 100, 101, 127, 126
98, 101, 102, 128, 127
99, 102, 103, 129, 128
100, 103, 104, 130, 129
101, 105, 106, 132, 131
102, 106, 107, 133, 132
103, 107, 108, 134, 133
104, 108, 109, 135, 134
105, 109, 110, 136, 135
106, 110, 111, 137, 136
107, 111, 112, 138, 137
108, 112, 113, 139, 138
109, 113, 114, 140, 139
110, 114, 115, 141, 140
111, 115, 116, 142, 141
112, 116, 117, 143, 142
113, 117, 118, 144, 143
114, 118, 119, 145, 144
115, 119, 120, 146, 145
116, 120, 121, 147, 146
117, 121, 122, 148, 147
118, 122, 123, 149, 148
119, 123, 124, 150, 149
120, 124, 125, 151, 150
121, 125, 126, 152, 151
122, 126, 127, 153, 152
123, 127, 128, 154, 153
124, 128, 129, 155, 154
125, 129, 130, 156, 155
126, 131, 132, 158, 157
127, 132, 133, 159, 158
128, 133, 134, 160, 159
129, 134, 135, 161, 160
130, 135, 136, 162, 161
131, 136, 137, 163, 162
132, 137, 138, 164, 163
133, 138, 139, 165, 164
134, 139, 140, 166, 165
135, 140, 141, 167, 166
136, 141, 142, 168, 167
137, 142, 143, 169, 168
138, 143, 144, 170, 169
139, 144, 145, 171, 170
140, 145, 146, 172, 171
141, 146, 147, 173, 172
142, 147, 148, 174, 173
143, 148, 149, 175, 174
144, 149, 150, 176, 175
145, 150, 151, 177, 176
146, 151, 152, 178, 177
147, 152, 153, 179, 178
148, 153, 154, 180, 179
149, 154, 155, 181, 180
150, 155, 156, 182, 181
*Nset, nset=Set-1, generate
1, 182, 1
*Elset, elset=Set-1, generate
1, 150, 1
*Nset, nset=_PickedSet5, internal, generate
1, 182, 1
*Elset, elset=_PickedSet5, internal, generate
1, 150, 1
*Nset, nset=_PickedSet8, internal, generate
1, 182, 1
*Elset, elset=_PickedSet8, internal, generate
1, 150, 1
** Region: (CompositeLayup-1-1: Generated From Layup)
*Elset, elset=CompositeLayup-1-1, generate
1, 150, 1
** Section: CompositeLayup-1-1
*Shell Section, elset=CompositeLayup-1-1, composite, layup=CompositeLayup-1
1., 3, Carbone, 90., Ply-1
*Transverse Shear
3000., 3000., 0.
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*End Instance
**
*Nset, nset=Set-1, instance=Part-1-1
182,
*Nset, nset=Set-2, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-2, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-3, instance=Part-1-1, generate
157, 182, 1
*Elset, elset=Set-3, instance=Part-1-1, generate
126, 150, 1
*Nset, nset=Set-4, instance=Part-1-1, generate
1, 26, 1
*Elset, elset=Set-4, instance=Part-1-1, generate
1, 25, 1
*Nset, nset=Set-5, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-5, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-6, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-6, instance=Part-1-1, generate
1, 126, 25
*Nset, nset=Set-7, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-7, instance=Part-1-1, generate
1, 126, 25
*Elset, elset=_Surf-1_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_E4, E4
*Elset, elset=_Surf-2_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-2
_Surf-2_E4, E4
*Elset, elset=_Surf-3_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-3
_Surf-3_E4, E4
*Elset, elset=_Surf-4_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-4
_Surf-4_E4, E4
*End Assembly
**
** MATERIALS
**
*Material, name=Carbone
*Depvar
23,
*User Material, constants=15
130000.,9000.,9000., 0.36, 0.3, 0.3,3600.,3600.
3100.,2500.,1000., 50., 65., 30., 30.
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES
*Static
1., 1., 1e-05, 1.
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
Set-1, ENCASTRE
** Name: BC-2 Type: Displacement/Rotation
*Boundary
Set-2, 1, 1
Set-2, 2, 2
Set-2, 3, 3
Set-2, 4, 4
Set-2, 5, 5
Set-2, 6, 6
**
** LOADS
**
** Name: Load-1 Type: Shell edge load
*Dsload
Surf-2, EDNOR, -10.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
PHILSM, PSILSM, U
*Element Output, directions=YES
CFAILURE, DAMAGEC, DAMAGEFC, DAMAGEFT, DAMAGEMC, DAMAGEMT, DAMAGESHR, DAMAGET, DMICRT, E, ERPRATIO, HSNFCCRT, HSNFTCRT, HSNMCCRT, HSNMTCRT, JK
S, SDEG, SHRRATIO
*Contact Output
BDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history
*Energy Output
ALLAE, ALLCD, ALLDMD, ALLEE, ALLFD, ALLIE, ALLJD, ALLKE, ALLKL, ALLPD, ALLQB, ALLSD, ALLSE, ALLVD, ALLWK, ETOTAL
*Contact Output
BDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC
*End Step
And the routine FORTRAN:
C ---------------------------------------------------------------------
C
C LIGNES STANDARD SUBROUTINE UMAT
C
C
SUBROUTINE UMAT(STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL,
1 DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, TIME, DTIME, TEMP, DTEMP,
2 PREDEF, DPRED, CMNAME, NDI, NSHR, NTENS, NSTATV, PROPS, NPROPS,
3 COORDS, DROT, PNEWDT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, LAYER,
4 KSPT, KSTEP, KINC)
C
INCLUDE 'ABA_PARAM.INC'
CHARACTER *80 CMNAME
DIMENSION STRESS(NTENS), STATEV(NSTATV), DDSDDE(NTENS, NTENS),
1 DDSDDT(NTENS), DRPLDE(NTENS), STRAN(NTENS), DSTRAN(NTENS),
2 PREDEF(1), DPRED(1), PROPS(NPROPS), COORDS(3), DROT(3, 3),
3 DFGRD0(3, 3), DFGRD1(3, 3), EELAS(6), EPLAS(6), FLOW(6)
C Attention !!, ne pas faire commencer mes variables par i, j, k, l, m, n sinon la valeur sortie est automatiquement la partie entière
C Material properties
E1 = PROPS(1)
E2 = PROPS(2)
E3 = PROPS(3)
PU12 = PROPS(4)
PU13 = PROPS(5)
PU23 = PROPS(6)
G12 = PROPS(7)
G13 = PROPS(8)
G23 = PROPS(9)
XT = PROPS(10)
XC = PROPS(11)
YT = PROPS(12)
YC = PROPS(13)
SL = PROPS(14)
ST = PROPS(15)
C calcul du NU21
PU21=((E2 * PU12) / E1)
C calcul du NU31
PU31=(E3*PU13)/E1
C calcul du NU32
PU32=(E3*PU23)/E2
C calcul du paramètre gamma
GAMMA = 1/(1-(PU12*PU21)-(PU23*PU32)-(PU31*PU13)-(2*PU21*PU32*PU13))
C ================================================================
C ELASTIC STIFFNESS MATRIX (in terms of E, v)
C-----------------------------------------------------------------------------------------------------
C [ D11 D21 D13 0 0 0 ]
C [ D12 D22 D23 0 0 0 ]
C [ D13 D32 D33 0 0 0 ]
C [ 0 0 0 D44 0 0 ]
C [ 0 0 0 0 D55 0 ]
C [ 0 0 0 0 0 D66]
C calcul des constantes de l'ingénieur
C DDSDDE(1,1)=E1*(1-(NU23*NU32))*GAMMA
C DDSDDE(2,2)=E2*(1-(NU13*NU31))*GAMMA
C DDSDDE(3,3)=E3*(1-(NU12*NU21))*GAMMA
C DDSDDE(1,2)=E2*(NU12+(NU32*NU13))*GAMMA
C DDSDDE(2,1)=DDSDDE(1,2)
C DDSDDE(1,3)=E3*(NU13+(NU12*NU23))*GAMMA
C DDSDDE(3,1)=DDSDDE(1,3)
C DDSDDE(2,3)=E3*(NU23+(NU21*NU13))*GAMMA
C DDSDDE(3,2)=DDSDDE(2,3)
C DDSDDE(4,4)=G12
C DDSDDE(5,5)=G13
C DDSDDE(6,6)=G23
DO K1=1,NDI
IF (K1 .EQ. 1 .AND. NDI .EQ. 2) THEN
DDSDDE(1,K1)=E1*(1-(PU23*PU32))*GAMMA
DDSDDE(2,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(3,K1)=0
ELSE
IF (K1 .EQ. 1 .AND. NDI .EQ. 3) THEN
DDSDDE(1,K1)=E1*(1-(PU23*PU32))*GAMMA
DDSDDE(2,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(3,K1)=E3*(PU13+(PU12*PU23))*GAMMA
ELSE
IF (K1 .EQ. 2 .AND. NDI .EQ. 2) THEN
DDSDDE(1,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(2,K1)=E2*(1-(PU13*PU31))*GAMMA
DDSDDE(3,K1)=0
ELSE
IF (K1 .EQ. 2 .AND. NDI .EQ. 3) THEN
DDSDDE(1,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(2,K1)=E2*(1-(PU13*PU31))*GAMMA
DDSDDE(3,K1)=E3*(PU23+(PU21*PU13))*GAMMA
ELSE
IF (K1 .EQ. 3 .AND. NDI .EQ. 3) THEN
DDSDDE(1,K1)=E3*(PU13+(PU12*PU23))*GAMMA
DDSDDE(2,K1)=E3*(PU23+(PU21*PU13))*GAMMA
DDSDDE(3,K1)=E3*(1-(PU12*PU21))*GAMMA
END IF
END IF
END IF
END IF
END IF
END DO
DO K1=NDI+1,NTENS
IF (K1 .EQ. 3 .AND. NTENS .EQ. 6) THEN
DDSDDE(K1,K1)=E3*(1-(PU12*PU21))*GAMMA
ELSE
IF (K1 .EQ. 4) THEN
DDSDDE(K1,K1)=G12
ELSE
IF (K1 .EQ. 5) THEN
DDSDDE(K1,K1)=G13
ELSE
IF (K1 .EQ. NTENS .AND. NDI .EQ. 2 .OR. NDI .EQ. 3) THEN
DDSDDE(K1,K1)=G23
END IF
END IF
END IF
END IF
END DO
C ================================================================
C STRESS CALCULATION
C-----------------------------------------------------------------------------------------------------
C [Sig11] [ D11 D12 D13 0 0 0 ] [eps11]11
C [Sig22] [ D21 D22 D23 0 0 0 ] [eps22]
C [Sig33] [ D31 D32 D33 0 0 0 ] [eps33]
C [Tau12] [ 0 0 0 D44 0 0 ] [gam12]
C [Tau13] [ 0 0 0 0 D55 0 ] [gam13]
C [Tau23] [ 0 0 0 0 0 D66] [gam23]
C-----------------------------------------------------------------------------------------------------
C DO K1=1,NTENS
C DO K2=1,NTENS
C STRESS(K2)=STRESS(K2)+(DDSDDE(K2,K1)*DSTRAN(K1))
C END DO
C END DO
IF (NTENS .EQ. 3) THEN
STRESS(1)=(DDSDDE(1,1)*DSTRAN(1))+(DDSDDE(1,2)*DSTRAN(2))
STRESS(2)=(DDSDDE(2,1)*DSTRAN(1))+(DDSDDE(2,2)*DSTRAN(2))
STRESS(3)=(DDSDDE(3,1)*DSTRAN(1))+(DDSDDE(3,2)*DSTRAN(2))
1 +(DDSDDE(3,3)*DSTRAN(3))
ELSE
IF (NTENS .EQ. 6) THEN
STRESS(1)=(DDSDDE(1,1)*DSTRAN(1))+(DDSDDE(1,2)*DSTRAN(2))
STRESS(2)=(DDSDDE(2,1)*DSTRAN(1))+(DDSDDE(2,2)*DSTRAN(2))
STRESS(3)=(DDSDDE(3,1)*DSTRAN(1))+(DDSDDE(3,2)*DSTRAN(2))
1 +(DDSDDE(3,3)*DSTRAN(3))
STRESS(4)=(DDSDDE(4,4)*DSTRAN(4))
STRESS(5)=(DDSDDE(5,5)*DSTRAN(5))
STRESS(6)=(DDSDDE(6,6)*DSTRAN(6))
END IF
END IF
C HASHIN FIBRE
IF (STRESS(1) .GE. 0) THEN
uFTF=(STRESS(1)/XT)**2
END IF
IF (STRESS(1) .LT. 0) THEN
uFCF=(STRESS(1)/XC)**2
END IF
C HASHIN MATRICE
IF (STRESS(2) .GE. 0 .AND. NTENS .EQ. 3) THEN
uFTM=((STRESS(2)/YT)**2)+((STRESS(3)/SL)**2)
IF (STRESS(2) .GE. 0 .AND. NTENS .EQ. 6 ) THEN
uFTM=((STRESS(2)/YT)**2)+((STRESS(6)/SL)**2)
END IF
END IF
IF (STRESS(2) .LT. 0) THEN
uFCM=((STRESS(2)/(2*ST))**2)+((((YC/(2*ST))**2)-1)*(STRESS(2)/YC))+
* (STRESS(3)/SL)**2
END IF
STATEV(1) = STRESS(1)
STATEV(2) = STRESS(2)
STATEV(3) = STRESS(3)
STATEV(4) = STRESS(5)
STATEV(5) = STRESS(6)
STATEV(6) = DDSDDE(1,1)
STATEV(7) = DDSDDE(1,2)
STATEV(8) = DDSDDE(1,3)
STATEV(9) = DDSDDE(2,1)
STATEV(10) = DDSDDE(2,2)
STATEV(11) = DDSDDE(2,3)
STATEV(12) = DDSDDE(3,1)
STATEV(13) = DDSDDE(3,2)
STATEV(14) = DDSDDE(3,3)
STATEV(15) = DDSDDE(4,4)
STATEV(16) = DDSDDE(5,5)
STATEV(17) = DDSDDE(6,6)
STATEV(18) = DSTRAN(1)
STATEV(19) = DSTRAN(2)
STATEV(20) = DSTRAN(3)
STATEV(21) = DSTRAN(4)
STATEV(22) = DSTRAN(5)
STATEV(23) = DSTRAN(6)
C-----------------------------------------------------------------------------------------------------
RETURN
END
The STRESSES (S11,S22,S12) and the STRAINS (E11,E22,E12) are not the same in the both model, why ?
How could you explain that ?
Please help me !
I would like to compute an anisotropic linear elastic law for composite.
I wish to output stress and strain in order to calculate 2D and 3D failure criteria (Hashin 3D, Puck...)
I am a beginner in UMAT computing.
When I run a calculation without UMAT, I do not find the same thing when I do the calculation with the UMAT I computed.
For instance, please find below my the .inp without the UMAT:
*Heading
** Job name: sans_umat_90degre_2D Model name: sans_umat_orthotrope_2D
** Generated by: Abaqus/CAE 6.12-3
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
1, 55., 25., 0.
2, 51., 25., 0.
3, 47., 25., 0.
4, 43., 25., 0.
5, 39., 25., 0.
6, 35., 25., 0.
7, 31., 25., 0.
8, 27., 25., 0.
9, 23., 25., 0.
10, 19., 25., 0.
11, 15., 25., 0.
12, 11., 25., 0.
13, 7., 25., 0.
14, 3., 25., 0.
15, -1., 25., 0.
16, -5., 25., 0.
17, -9., 25., 0.
18, -13., 25., 0.
19, -17., 25., 0.
20, -21., 25., 0.
21, -25., 25., 0.
22, -29., 25., 0.
23, -33., 25., 0.
24, -37., 25., 0.
25, -41., 25., 0.
26, -45., 25., 0.
27, 55., 20.833334, 0.
28, 51., 20.833334, 0.
29, 47., 20.833334, 0.
30, 43., 20.833334, 0.
31, 39., 20.833334, 0.
32, 35., 20.833334, 0.
33, 31., 20.833334, 0.
34, 27., 20.833334, 0.
35, 23., 20.833334, 0.
36, 19., 20.833334, 0.
37, 15., 20.833334, 0.
38, 11., 20.833334, 0.
39, 7., 20.833334, 0.
40, 3., 20.833334, 0.
41, -1., 20.833334, 0.
42, -5., 20.833334, 0.
43, -9., 20.833334, 0.
44, -13., 20.833334, 0.
45, -17., 20.833334, 0.
46, -21., 20.833334, 0.
47, -25., 20.833334, 0.
48, -29., 20.833334, 0.
49, -33., 20.833334, 0.
50, -37., 20.833334, 0.
51, -41., 20.833334, 0.
52, -45., 20.833334, 0.
53, 55., 16.666666, 0.
54, 51., 16.666666, 0.
55, 47., 16.666666, 0.
56, 43., 16.666666, 0.
57, 39., 16.666666, 0.
58, 35., 16.666666, 0.
59, 31., 16.666666, 0.
60, 27., 16.666666, 0.
61, 23., 16.666666, 0.
62, 19., 16.666666, 0.
63, 15., 16.666666, 0.
64, 11., 16.666666, 0.
65, 7., 16.666666, 0.
66, 3., 16.666666, 0.
67, -1., 16.666666, 0.
68, -5., 16.666666, 0.
69, -9., 16.666666, 0.
70, -13., 16.666666, 0.
71, -17., 16.666666, 0.
72, -21., 16.666666, 0.
73, -25., 16.666666, 0.
74, -29., 16.666666, 0.
75, -33., 16.666666, 0.
76, -37., 16.666666, 0.
77, -41., 16.666666, 0.
78, -45., 16.666666, 0.
79, 55., 12.5, 0.
80, 51., 12.5, 0.
81, 47., 12.5, 0.
82, 43., 12.5, 0.
83, 39., 12.5, 0.
84, 35., 12.5, 0.
85, 31., 12.5, 0.
86, 27., 12.5, 0.
87, 23., 12.5, 0.
88, 19., 12.5, 0.
89, 15., 12.5, 0.
90, 11., 12.5, 0.
91, 7., 12.5, 0.
92, 3., 12.5, 0.
93, -1., 12.5, 0.
94, -5., 12.5, 0.
95, -9., 12.5, 0.
96, -13., 12.5, 0.
97, -17., 12.5, 0.
98, -21., 12.5, 0.
99, -25., 12.5, 0.
100, -29., 12.5, 0.
101, -33., 12.5, 0.
102, -37., 12.5, 0.
103, -41., 12.5, 0.
104, -45., 12.5, 0.
105, 55., 8.33333302, 0.
106, 51., 8.33333302, 0.
107, 47., 8.33333302, 0.
108, 43., 8.33333302, 0.
109, 39., 8.33333302, 0.
110, 35., 8.33333302, 0.
111, 31., 8.33333302, 0.
112, 27., 8.33333302, 0.
113, 23., 8.33333302, 0.
114, 19., 8.33333302, 0.
115, 15., 8.33333302, 0.
116, 11., 8.33333302, 0.
117, 7., 8.33333302, 0.
118, 3., 8.33333302, 0.
119, -1., 8.33333302, 0.
120, -5., 8.33333302, 0.
121, -9., 8.33333302, 0.
122, -13., 8.33333302, 0.
123, -17., 8.33333302, 0.
124, -21., 8.33333302, 0.
125, -25., 8.33333302, 0.
126, -29., 8.33333302, 0.
127, -33., 8.33333302, 0.
128, -37., 8.33333302, 0.
129, -41., 8.33333302, 0.
130, -45., 8.33333302, 0.
131, 55., 4.16666651, 0.
132, 51., 4.16666651, 0.
133, 47., 4.16666651, 0.
134, 43., 4.16666651, 0.
135, 39., 4.16666651, 0.
136, 35., 4.16666651, 0.
137, 31., 4.16666651, 0.
138, 27., 4.16666651, 0.
139, 23., 4.16666651, 0.
140, 19., 4.16666651, 0.
141, 15., 4.16666651, 0.
142, 11., 4.16666651, 0.
143, 7., 4.16666651, 0.
144, 3., 4.16666651, 0.
145, -1., 4.16666651, 0.
146, -5., 4.16666651, 0.
147, -9., 4.16666651, 0.
148, -13., 4.16666651, 0.
149, -17., 4.16666651, 0.
150, -21., 4.16666651, 0.
151, -25., 4.16666651, 0.
152, -29., 4.16666651, 0.
153, -33., 4.16666651, 0.
154, -37., 4.16666651, 0.
155, -41., 4.16666651, 0.
156, -45., 4.16666651, 0.
157, 55., 0., 0.
158, 51., 0., 0.
159, 47., 0., 0.
160, 43., 0., 0.
161, 39., 0., 0.
162, 35., 0., 0.
163, 31., 0., 0.
164, 27., 0., 0.
165, 23., 0., 0.
166, 19., 0., 0.
167, 15., 0., 0.
168, 11., 0., 0.
169, 7., 0., 0.
170, 3., 0., 0.
171, -1., 0., 0.
172, -5., 0., 0.
173, -9., 0., 0.
174, -13., 0., 0.
175, -17., 0., 0.
176, -21., 0., 0.
177, -25., 0., 0.
178, -29., 0., 0.
179, -33., 0., 0.
180, -37., 0., 0.
181, -41., 0., 0.
182, -45., 0., 0.
*Element, type=S4R
1, 1, 2, 28, 27
2, 2, 3, 29, 28
3, 3, 4, 30, 29
4, 4, 5, 31, 30
5, 5, 6, 32, 31
6, 6, 7, 33, 32
7, 7, 8, 34, 33
8, 8, 9, 35, 34
9, 9, 10, 36, 35
10, 10, 11, 37, 36
11, 11, 12, 38, 37
12, 12, 13, 39, 38
13, 13, 14, 40, 39
14, 14, 15, 41, 40
15, 15, 16, 42, 41
16, 16, 17, 43, 42
17, 17, 18, 44, 43
18, 18, 19, 45, 44
19, 19, 20, 46, 45
20, 20, 21, 47, 46
21, 21, 22, 48, 47
22, 22, 23, 49, 48
23, 23, 24, 50, 49
24, 24, 25, 51, 50
25, 25, 26, 52, 51
26, 27, 28, 54, 53
27, 28, 29, 55, 54
28, 29, 30, 56, 55
29, 30, 31, 57, 56
30, 31, 32, 58, 57
31, 32, 33, 59, 58
32, 33, 34, 60, 59
33, 34, 35, 61, 60
34, 35, 36, 62, 61
35, 36, 37, 63, 62
36, 37, 38, 64, 63
37, 38, 39, 65, 64
38, 39, 40, 66, 65
39, 40, 41, 67, 66
40, 41, 42, 68, 67
41, 42, 43, 69, 68
42, 43, 44, 70, 69
43, 44, 45, 71, 70
44, 45, 46, 72, 71
45, 46, 47, 73, 72
46, 47, 48, 74, 73
47, 48, 49, 75, 74
48, 49, 50, 76, 75
49, 50, 51, 77, 76
50, 51, 52, 78, 77
51, 53, 54, 80, 79
52, 54, 55, 81, 80
53, 55, 56, 82, 81
54, 56, 57, 83, 82
55, 57, 58, 84, 83
56, 58, 59, 85, 84
57, 59, 60, 86, 85
58, 60, 61, 87, 86
59, 61, 62, 88, 87
60, 62, 63, 89, 88
61, 63, 64, 90, 89
62, 64, 65, 91, 90
63, 65, 66, 92, 91
64, 66, 67, 93, 92
65, 67, 68, 94, 93
66, 68, 69, 95, 94
67, 69, 70, 96, 95
68, 70, 71, 97, 96
69, 71, 72, 98, 97
70, 72, 73, 99, 98
71, 73, 74, 100, 99
72, 74, 75, 101, 100
73, 75, 76, 102, 101
74, 76, 77, 103, 102
75, 77, 78, 104, 103
76, 79, 80, 106, 105
77, 80, 81, 107, 106
78, 81, 82, 108, 107
79, 82, 83, 109, 108
80, 83, 84, 110, 109
81, 84, 85, 111, 110
82, 85, 86, 112, 111
83, 86, 87, 113, 112
84, 87, 88, 114, 113
85, 88, 89, 115, 114
86, 89, 90, 116, 115
87, 90, 91, 117, 116
88, 91, 92, 118, 117
89, 92, 93, 119, 118
90, 93, 94, 120, 119
91, 94, 95, 121, 120
92, 95, 96, 122, 121
93, 96, 97, 123, 122
94, 97, 98, 124, 123
95, 98, 99, 125, 124
96, 99, 100, 126, 125
97, 100, 101, 127, 126
98, 101, 102, 128, 127
99, 102, 103, 129, 128
100, 103, 104, 130, 129
101, 105, 106, 132, 131
102, 106, 107, 133, 132
103, 107, 108, 134, 133
104, 108, 109, 135, 134
105, 109, 110, 136, 135
106, 110, 111, 137, 136
107, 111, 112, 138, 137
108, 112, 113, 139, 138
109, 113, 114, 140, 139
110, 114, 115, 141, 140
111, 115, 116, 142, 141
112, 116, 117, 143, 142
113, 117, 118, 144, 143
114, 118, 119, 145, 144
115, 119, 120, 146, 145
116, 120, 121, 147, 146
117, 121, 122, 148, 147
118, 122, 123, 149, 148
119, 123, 124, 150, 149
120, 124, 125, 151, 150
121, 125, 126, 152, 151
122, 126, 127, 153, 152
123, 127, 128, 154, 153
124, 128, 129, 155, 154
125, 129, 130, 156, 155
126, 131, 132, 158, 157
127, 132, 133, 159, 158
128, 133, 134, 160, 159
129, 134, 135, 161, 160
130, 135, 136, 162, 161
131, 136, 137, 163, 162
132, 137, 138, 164, 163
133, 138, 139, 165, 164
134, 139, 140, 166, 165
135, 140, 141, 167, 166
136, 141, 142, 168, 167
137, 142, 143, 169, 168
138, 143, 144, 170, 169
139, 144, 145, 171, 170
140, 145, 146, 172, 171
141, 146, 147, 173, 172
142, 147, 148, 174, 173
143, 148, 149, 175, 174
144, 149, 150, 176, 175
145, 150, 151, 177, 176
146, 151, 152, 178, 177
147, 152, 153, 179, 178
148, 153, 154, 180, 179
149, 154, 155, 181, 180
150, 155, 156, 182, 181
*Nset, nset=Set-1, generate
1, 182, 1
*Elset, elset=Set-1, generate
1, 150, 1
*Nset, nset=_PickedSet5, internal, generate
1, 182, 1
*Elset, elset=_PickedSet5, internal, generate
1, 150, 1
** Region: (CompositeLayup-1-1: Generated From Layup)
*Elset, elset=CompositeLayup-1-1, generate
1, 150, 1
** Section: CompositeLayup-1-1
*Shell Section, elset=CompositeLayup-1-1, composite, layup=CompositeLayup-1
1., 3, Carbone, 90., Ply-1
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*End Instance
**
*Nset, nset=Set-1, instance=Part-1-1
182,
*Nset, nset=Set-2, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-2, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-3, instance=Part-1-1, generate
157, 182, 1
*Elset, elset=Set-3, instance=Part-1-1, generate
126, 150, 1
*Nset, nset=Set-4, instance=Part-1-1, generate
1, 26, 1
*Elset, elset=Set-4, instance=Part-1-1, generate
1, 25, 1
*Nset, nset=Set-5, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-5, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-6, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-6, instance=Part-1-1, generate
1, 126, 25
*Nset, nset=Set-7, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-7, instance=Part-1-1, generate
1, 126, 25
*Elset, elset=_Surf-1_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_E4, E4
*Elset, elset=_Surf-2_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-2
_Surf-2_E4, E4
*Elset, elset=_Surf-3_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-3
_Surf-3_E4, E4
*Elset, elset=_Surf-4_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-4
_Surf-4_E4, E4
*End Assembly
**
** MATERIALS
**
*Material, name=Carbone
*Damage Initiation, criterion=HASHIN
2500.,1000., 50., 65., 30., 30.
*Elastic, type=ENGINEERING CONSTANTS
130000.,9000.,9000., 0.36, 0.3, 0.3,3600.,3600.
3100.,
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES
*Static
1., 1., 1e-05, 1.
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
Set-1, ENCASTRE
** Name: BC-2 Type: Displacement/Rotation
*Boundary
Set-2, 1, 1
Set-2, 2, 2
Set-2, 3, 3
Set-2, 4, 4
Set-2, 5, 5
Set-2, 6, 6
**
** LOADS
**
** Name: Load-1 Type: Shell edge load
*Dsload
Surf-2, EDNOR, -10.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
U,
*Element Output, directions=YES
E, LE, NE, S
**
** FIELD OUTPUT: F-Output-2
**
*Element Output, elset=Part-1-1.CompositeLayup-1-1, directions=YES
1, 2, 3
E, LE, NE, S
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history
*Energy Output
ALLAE, ALLCD, ALLDMD, ALLEE, ALLFD, ALLIE, ALLJD, ALLKE, ALLKL, ALLPD, ALLQB, ALLSD, ALLSE, ALLVD, ALLWK, ETOTAL
*Contact Output
BDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC
*End Step
And the .inp with the UMAT (you have to copy/paste in a .inp file to run in Abaqus and launch the calculation):
*Heading
** Job name: avec_umat_90degre_2D Model name: avec_umat_orthotrope_2D
** Generated by: Abaqus/CAE 6.12-3
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
1, 55., 25., 0.
2, 51., 25., 0.
3, 47., 25., 0.
4, 43., 25., 0.
5, 39., 25., 0.
6, 35., 25., 0.
7, 31., 25., 0.
8, 27., 25., 0.
9, 23., 25., 0.
10, 19., 25., 0.
11, 15., 25., 0.
12, 11., 25., 0.
13, 7., 25., 0.
14, 3., 25., 0.
15, -1., 25., 0.
16, -5., 25., 0.
17, -9., 25., 0.
18, -13., 25., 0.
19, -17., 25., 0.
20, -21., 25., 0.
21, -25., 25., 0.
22, -29., 25., 0.
23, -33., 25., 0.
24, -37., 25., 0.
25, -41., 25., 0.
26, -45., 25., 0.
27, 55., 20.833334, 0.
28, 51., 20.833334, 0.
29, 47., 20.833334, 0.
30, 43., 20.833334, 0.
31, 39., 20.833334, 0.
32, 35., 20.833334, 0.
33, 31., 20.833334, 0.
34, 27., 20.833334, 0.
35, 23., 20.833334, 0.
36, 19., 20.833334, 0.
37, 15., 20.833334, 0.
38, 11., 20.833334, 0.
39, 7., 20.833334, 0.
40, 3., 20.833334, 0.
41, -1., 20.833334, 0.
42, -5., 20.833334, 0.
43, -9., 20.833334, 0.
44, -13., 20.833334, 0.
45, -17., 20.833334, 0.
46, -21., 20.833334, 0.
47, -25., 20.833334, 0.
48, -29., 20.833334, 0.
49, -33., 20.833334, 0.
50, -37., 20.833334, 0.
51, -41., 20.833334, 0.
52, -45., 20.833334, 0.
53, 55., 16.666666, 0.
54, 51., 16.666666, 0.
55, 47., 16.666666, 0.
56, 43., 16.666666, 0.
57, 39., 16.666666, 0.
58, 35., 16.666666, 0.
59, 31., 16.666666, 0.
60, 27., 16.666666, 0.
61, 23., 16.666666, 0.
62, 19., 16.666666, 0.
63, 15., 16.666666, 0.
64, 11., 16.666666, 0.
65, 7., 16.666666, 0.
66, 3., 16.666666, 0.
67, -1., 16.666666, 0.
68, -5., 16.666666, 0.
69, -9., 16.666666, 0.
70, -13., 16.666666, 0.
71, -17., 16.666666, 0.
72, -21., 16.666666, 0.
73, -25., 16.666666, 0.
74, -29., 16.666666, 0.
75, -33., 16.666666, 0.
76, -37., 16.666666, 0.
77, -41., 16.666666, 0.
78, -45., 16.666666, 0.
79, 55., 12.5, 0.
80, 51., 12.5, 0.
81, 47., 12.5, 0.
82, 43., 12.5, 0.
83, 39., 12.5, 0.
84, 35., 12.5, 0.
85, 31., 12.5, 0.
86, 27., 12.5, 0.
87, 23., 12.5, 0.
88, 19., 12.5, 0.
89, 15., 12.5, 0.
90, 11., 12.5, 0.
91, 7., 12.5, 0.
92, 3., 12.5, 0.
93, -1., 12.5, 0.
94, -5., 12.5, 0.
95, -9., 12.5, 0.
96, -13., 12.5, 0.
97, -17., 12.5, 0.
98, -21., 12.5, 0.
99, -25., 12.5, 0.
100, -29., 12.5, 0.
101, -33., 12.5, 0.
102, -37., 12.5, 0.
103, -41., 12.5, 0.
104, -45., 12.5, 0.
105, 55., 8.33333302, 0.
106, 51., 8.33333302, 0.
107, 47., 8.33333302, 0.
108, 43., 8.33333302, 0.
109, 39., 8.33333302, 0.
110, 35., 8.33333302, 0.
111, 31., 8.33333302, 0.
112, 27., 8.33333302, 0.
113, 23., 8.33333302, 0.
114, 19., 8.33333302, 0.
115, 15., 8.33333302, 0.
116, 11., 8.33333302, 0.
117, 7., 8.33333302, 0.
118, 3., 8.33333302, 0.
119, -1., 8.33333302, 0.
120, -5., 8.33333302, 0.
121, -9., 8.33333302, 0.
122, -13., 8.33333302, 0.
123, -17., 8.33333302, 0.
124, -21., 8.33333302, 0.
125, -25., 8.33333302, 0.
126, -29., 8.33333302, 0.
127, -33., 8.33333302, 0.
128, -37., 8.33333302, 0.
129, -41., 8.33333302, 0.
130, -45., 8.33333302, 0.
131, 55., 4.16666651, 0.
132, 51., 4.16666651, 0.
133, 47., 4.16666651, 0.
134, 43., 4.16666651, 0.
135, 39., 4.16666651, 0.
136, 35., 4.16666651, 0.
137, 31., 4.16666651, 0.
138, 27., 4.16666651, 0.
139, 23., 4.16666651, 0.
140, 19., 4.16666651, 0.
141, 15., 4.16666651, 0.
142, 11., 4.16666651, 0.
143, 7., 4.16666651, 0.
144, 3., 4.16666651, 0.
145, -1., 4.16666651, 0.
146, -5., 4.16666651, 0.
147, -9., 4.16666651, 0.
148, -13., 4.16666651, 0.
149, -17., 4.16666651, 0.
150, -21., 4.16666651, 0.
151, -25., 4.16666651, 0.
152, -29., 4.16666651, 0.
153, -33., 4.16666651, 0.
154, -37., 4.16666651, 0.
155, -41., 4.16666651, 0.
156, -45., 4.16666651, 0.
157, 55., 0., 0.
158, 51., 0., 0.
159, 47., 0., 0.
160, 43., 0., 0.
161, 39., 0., 0.
162, 35., 0., 0.
163, 31., 0., 0.
164, 27., 0., 0.
165, 23., 0., 0.
166, 19., 0., 0.
167, 15., 0., 0.
168, 11., 0., 0.
169, 7., 0., 0.
170, 3., 0., 0.
171, -1., 0., 0.
172, -5., 0., 0.
173, -9., 0., 0.
174, -13., 0., 0.
175, -17., 0., 0.
176, -21., 0., 0.
177, -25., 0., 0.
178, -29., 0., 0.
179, -33., 0., 0.
180, -37., 0., 0.
181, -41., 0., 0.
182, -45., 0., 0.
*Element, type=S4R
1, 1, 2, 28, 27
2, 2, 3, 29, 28
3, 3, 4, 30, 29
4, 4, 5, 31, 30
5, 5, 6, 32, 31
6, 6, 7, 33, 32
7, 7, 8, 34, 33
8, 8, 9, 35, 34
9, 9, 10, 36, 35
10, 10, 11, 37, 36
11, 11, 12, 38, 37
12, 12, 13, 39, 38
13, 13, 14, 40, 39
14, 14, 15, 41, 40
15, 15, 16, 42, 41
16, 16, 17, 43, 42
17, 17, 18, 44, 43
18, 18, 19, 45, 44
19, 19, 20, 46, 45
20, 20, 21, 47, 46
21, 21, 22, 48, 47
22, 22, 23, 49, 48
23, 23, 24, 50, 49
24, 24, 25, 51, 50
25, 25, 26, 52, 51
26, 27, 28, 54, 53
27, 28, 29, 55, 54
28, 29, 30, 56, 55
29, 30, 31, 57, 56
30, 31, 32, 58, 57
31, 32, 33, 59, 58
32, 33, 34, 60, 59
33, 34, 35, 61, 60
34, 35, 36, 62, 61
35, 36, 37, 63, 62
36, 37, 38, 64, 63
37, 38, 39, 65, 64
38, 39, 40, 66, 65
39, 40, 41, 67, 66
40, 41, 42, 68, 67
41, 42, 43, 69, 68
42, 43, 44, 70, 69
43, 44, 45, 71, 70
44, 45, 46, 72, 71
45, 46, 47, 73, 72
46, 47, 48, 74, 73
47, 48, 49, 75, 74
48, 49, 50, 76, 75
49, 50, 51, 77, 76
50, 51, 52, 78, 77
51, 53, 54, 80, 79
52, 54, 55, 81, 80
53, 55, 56, 82, 81
54, 56, 57, 83, 82
55, 57, 58, 84, 83
56, 58, 59, 85, 84
57, 59, 60, 86, 85
58, 60, 61, 87, 86
59, 61, 62, 88, 87
60, 62, 63, 89, 88
61, 63, 64, 90, 89
62, 64, 65, 91, 90
63, 65, 66, 92, 91
64, 66, 67, 93, 92
65, 67, 68, 94, 93
66, 68, 69, 95, 94
67, 69, 70, 96, 95
68, 70, 71, 97, 96
69, 71, 72, 98, 97
70, 72, 73, 99, 98
71, 73, 74, 100, 99
72, 74, 75, 101, 100
73, 75, 76, 102, 101
74, 76, 77, 103, 102
75, 77, 78, 104, 103
76, 79, 80, 106, 105
77, 80, 81, 107, 106
78, 81, 82, 108, 107
79, 82, 83, 109, 108
80, 83, 84, 110, 109
81, 84, 85, 111, 110
82, 85, 86, 112, 111
83, 86, 87, 113, 112
84, 87, 88, 114, 113
85, 88, 89, 115, 114
86, 89, 90, 116, 115
87, 90, 91, 117, 116
88, 91, 92, 118, 117
89, 92, 93, 119, 118
90, 93, 94, 120, 119
91, 94, 95, 121, 120
92, 95, 96, 122, 121
93, 96, 97, 123, 122
94, 97, 98, 124, 123
95, 98, 99, 125, 124
96, 99, 100, 126, 125
97, 100, 101, 127, 126
98, 101, 102, 128, 127
99, 102, 103, 129, 128
100, 103, 104, 130, 129
101, 105, 106, 132, 131
102, 106, 107, 133, 132
103, 107, 108, 134, 133
104, 108, 109, 135, 134
105, 109, 110, 136, 135
106, 110, 111, 137, 136
107, 111, 112, 138, 137
108, 112, 113, 139, 138
109, 113, 114, 140, 139
110, 114, 115, 141, 140
111, 115, 116, 142, 141
112, 116, 117, 143, 142
113, 117, 118, 144, 143
114, 118, 119, 145, 144
115, 119, 120, 146, 145
116, 120, 121, 147, 146
117, 121, 122, 148, 147
118, 122, 123, 149, 148
119, 123, 124, 150, 149
120, 124, 125, 151, 150
121, 125, 126, 152, 151
122, 126, 127, 153, 152
123, 127, 128, 154, 153
124, 128, 129, 155, 154
125, 129, 130, 156, 155
126, 131, 132, 158, 157
127, 132, 133, 159, 158
128, 133, 134, 160, 159
129, 134, 135, 161, 160
130, 135, 136, 162, 161
131, 136, 137, 163, 162
132, 137, 138, 164, 163
133, 138, 139, 165, 164
134, 139, 140, 166, 165
135, 140, 141, 167, 166
136, 141, 142, 168, 167
137, 142, 143, 169, 168
138, 143, 144, 170, 169
139, 144, 145, 171, 170
140, 145, 146, 172, 171
141, 146, 147, 173, 172
142, 147, 148, 174, 173
143, 148, 149, 175, 174
144, 149, 150, 176, 175
145, 150, 151, 177, 176
146, 151, 152, 178, 177
147, 152, 153, 179, 178
148, 153, 154, 180, 179
149, 154, 155, 181, 180
150, 155, 156, 182, 181
*Nset, nset=Set-1, generate
1, 182, 1
*Elset, elset=Set-1, generate
1, 150, 1
*Nset, nset=_PickedSet5, internal, generate
1, 182, 1
*Elset, elset=_PickedSet5, internal, generate
1, 150, 1
*Nset, nset=_PickedSet8, internal, generate
1, 182, 1
*Elset, elset=_PickedSet8, internal, generate
1, 150, 1
** Region: (CompositeLayup-1-1: Generated From Layup)
*Elset, elset=CompositeLayup-1-1, generate
1, 150, 1
** Section: CompositeLayup-1-1
*Shell Section, elset=CompositeLayup-1-1, composite, layup=CompositeLayup-1
1., 3, Carbone, 90., Ply-1
*Transverse Shear
3000., 3000., 0.
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*End Instance
**
*Nset, nset=Set-1, instance=Part-1-1
182,
*Nset, nset=Set-2, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-2, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-3, instance=Part-1-1, generate
157, 182, 1
*Elset, elset=Set-3, instance=Part-1-1, generate
126, 150, 1
*Nset, nset=Set-4, instance=Part-1-1, generate
1, 26, 1
*Elset, elset=Set-4, instance=Part-1-1, generate
1, 25, 1
*Nset, nset=Set-5, instance=Part-1-1, generate
26, 182, 26
*Elset, elset=Set-5, instance=Part-1-1, generate
25, 150, 25
*Nset, nset=Set-6, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-6, instance=Part-1-1, generate
1, 126, 25
*Nset, nset=Set-7, instance=Part-1-1, generate
1, 157, 26
*Elset, elset=Set-7, instance=Part-1-1, generate
1, 126, 25
*Elset, elset=_Surf-1_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_E4, E4
*Elset, elset=_Surf-2_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-2
_Surf-2_E4, E4
*Elset, elset=_Surf-3_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-3
_Surf-3_E4, E4
*Elset, elset=_Surf-4_E4, internal, instance=Part-1-1, generate
1, 126, 25
*Surface, type=ELEMENT, name=Surf-4
_Surf-4_E4, E4
*End Assembly
**
** MATERIALS
**
*Material, name=Carbone
*Depvar
23,
*User Material, constants=15
130000.,9000.,9000., 0.36, 0.3, 0.3,3600.,3600.
3100.,2500.,1000., 50., 65., 30., 30.
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES
*Static
1., 1., 1e-05, 1.
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
Set-1, ENCASTRE
** Name: BC-2 Type: Displacement/Rotation
*Boundary
Set-2, 1, 1
Set-2, 2, 2
Set-2, 3, 3
Set-2, 4, 4
Set-2, 5, 5
Set-2, 6, 6
**
** LOADS
**
** Name: Load-1 Type: Shell edge load
*Dsload
Surf-2, EDNOR, -10.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
PHILSM, PSILSM, U
*Element Output, directions=YES
CFAILURE, DAMAGEC, DAMAGEFC, DAMAGEFT, DAMAGEMC, DAMAGEMT, DAMAGESHR, DAMAGET, DMICRT, E, ERPRATIO, HSNFCCRT, HSNFTCRT, HSNMCCRT, HSNMTCRT, JK
S, SDEG, SHRRATIO
*Contact Output
BDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history
*Energy Output
ALLAE, ALLCD, ALLDMD, ALLEE, ALLFD, ALLIE, ALLJD, ALLKE, ALLKL, ALLPD, ALLQB, ALLSD, ALLSE, ALLVD, ALLWK, ETOTAL
*Contact Output
BDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC
*End Step
And the routine FORTRAN:
C ---------------------------------------------------------------------
C
C LIGNES STANDARD SUBROUTINE UMAT
C
C
SUBROUTINE UMAT(STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL,
1 DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, TIME, DTIME, TEMP, DTEMP,
2 PREDEF, DPRED, CMNAME, NDI, NSHR, NTENS, NSTATV, PROPS, NPROPS,
3 COORDS, DROT, PNEWDT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, LAYER,
4 KSPT, KSTEP, KINC)
C
INCLUDE 'ABA_PARAM.INC'
CHARACTER *80 CMNAME
DIMENSION STRESS(NTENS), STATEV(NSTATV), DDSDDE(NTENS, NTENS),
1 DDSDDT(NTENS), DRPLDE(NTENS), STRAN(NTENS), DSTRAN(NTENS),
2 PREDEF(1), DPRED(1), PROPS(NPROPS), COORDS(3), DROT(3, 3),
3 DFGRD0(3, 3), DFGRD1(3, 3), EELAS(6), EPLAS(6), FLOW(6)
C Attention !!, ne pas faire commencer mes variables par i, j, k, l, m, n sinon la valeur sortie est automatiquement la partie entière
C Material properties
E1 = PROPS(1)
E2 = PROPS(2)
E3 = PROPS(3)
PU12 = PROPS(4)
PU13 = PROPS(5)
PU23 = PROPS(6)
G12 = PROPS(7)
G13 = PROPS(8)
G23 = PROPS(9)
XT = PROPS(10)
XC = PROPS(11)
YT = PROPS(12)
YC = PROPS(13)
SL = PROPS(14)
ST = PROPS(15)
C calcul du NU21
PU21=((E2 * PU12) / E1)
C calcul du NU31
PU31=(E3*PU13)/E1
C calcul du NU32
PU32=(E3*PU23)/E2
C calcul du paramètre gamma
GAMMA = 1/(1-(PU12*PU21)-(PU23*PU32)-(PU31*PU13)-(2*PU21*PU32*PU13))
C ================================================================
C ELASTIC STIFFNESS MATRIX (in terms of E, v)
C-----------------------------------------------------------------------------------------------------
C [ D11 D21 D13 0 0 0 ]
C [ D12 D22 D23 0 0 0 ]
C [ D13 D32 D33 0 0 0 ]
C [ 0 0 0 D44 0 0 ]
C [ 0 0 0 0 D55 0 ]
C [ 0 0 0 0 0 D66]
C calcul des constantes de l'ingénieur
C DDSDDE(1,1)=E1*(1-(NU23*NU32))*GAMMA
C DDSDDE(2,2)=E2*(1-(NU13*NU31))*GAMMA
C DDSDDE(3,3)=E3*(1-(NU12*NU21))*GAMMA
C DDSDDE(1,2)=E2*(NU12+(NU32*NU13))*GAMMA
C DDSDDE(2,1)=DDSDDE(1,2)
C DDSDDE(1,3)=E3*(NU13+(NU12*NU23))*GAMMA
C DDSDDE(3,1)=DDSDDE(1,3)
C DDSDDE(2,3)=E3*(NU23+(NU21*NU13))*GAMMA
C DDSDDE(3,2)=DDSDDE(2,3)
C DDSDDE(4,4)=G12
C DDSDDE(5,5)=G13
C DDSDDE(6,6)=G23
DO K1=1,NDI
IF (K1 .EQ. 1 .AND. NDI .EQ. 2) THEN
DDSDDE(1,K1)=E1*(1-(PU23*PU32))*GAMMA
DDSDDE(2,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(3,K1)=0
ELSE
IF (K1 .EQ. 1 .AND. NDI .EQ. 3) THEN
DDSDDE(1,K1)=E1*(1-(PU23*PU32))*GAMMA
DDSDDE(2,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(3,K1)=E3*(PU13+(PU12*PU23))*GAMMA
ELSE
IF (K1 .EQ. 2 .AND. NDI .EQ. 2) THEN
DDSDDE(1,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(2,K1)=E2*(1-(PU13*PU31))*GAMMA
DDSDDE(3,K1)=0
ELSE
IF (K1 .EQ. 2 .AND. NDI .EQ. 3) THEN
DDSDDE(1,K1)=E2*(PU12+(PU32*PU13))*GAMMA
DDSDDE(2,K1)=E2*(1-(PU13*PU31))*GAMMA
DDSDDE(3,K1)=E3*(PU23+(PU21*PU13))*GAMMA
ELSE
IF (K1 .EQ. 3 .AND. NDI .EQ. 3) THEN
DDSDDE(1,K1)=E3*(PU13+(PU12*PU23))*GAMMA
DDSDDE(2,K1)=E3*(PU23+(PU21*PU13))*GAMMA
DDSDDE(3,K1)=E3*(1-(PU12*PU21))*GAMMA
END IF
END IF
END IF
END IF
END IF
END DO
DO K1=NDI+1,NTENS
IF (K1 .EQ. 3 .AND. NTENS .EQ. 6) THEN
DDSDDE(K1,K1)=E3*(1-(PU12*PU21))*GAMMA
ELSE
IF (K1 .EQ. 4) THEN
DDSDDE(K1,K1)=G12
ELSE
IF (K1 .EQ. 5) THEN
DDSDDE(K1,K1)=G13
ELSE
IF (K1 .EQ. NTENS .AND. NDI .EQ. 2 .OR. NDI .EQ. 3) THEN
DDSDDE(K1,K1)=G23
END IF
END IF
END IF
END IF
END DO
C ================================================================
C STRESS CALCULATION
C-----------------------------------------------------------------------------------------------------
C [Sig11] [ D11 D12 D13 0 0 0 ] [eps11]11
C [Sig22] [ D21 D22 D23 0 0 0 ] [eps22]
C [Sig33] [ D31 D32 D33 0 0 0 ] [eps33]
C [Tau12] [ 0 0 0 D44 0 0 ] [gam12]
C [Tau13] [ 0 0 0 0 D55 0 ] [gam13]
C [Tau23] [ 0 0 0 0 0 D66] [gam23]
C-----------------------------------------------------------------------------------------------------
C DO K1=1,NTENS
C DO K2=1,NTENS
C STRESS(K2)=STRESS(K2)+(DDSDDE(K2,K1)*DSTRAN(K1))
C END DO
C END DO
IF (NTENS .EQ. 3) THEN
STRESS(1)=(DDSDDE(1,1)*DSTRAN(1))+(DDSDDE(1,2)*DSTRAN(2))
STRESS(2)=(DDSDDE(2,1)*DSTRAN(1))+(DDSDDE(2,2)*DSTRAN(2))
STRESS(3)=(DDSDDE(3,1)*DSTRAN(1))+(DDSDDE(3,2)*DSTRAN(2))
1 +(DDSDDE(3,3)*DSTRAN(3))
ELSE
IF (NTENS .EQ. 6) THEN
STRESS(1)=(DDSDDE(1,1)*DSTRAN(1))+(DDSDDE(1,2)*DSTRAN(2))
STRESS(2)=(DDSDDE(2,1)*DSTRAN(1))+(DDSDDE(2,2)*DSTRAN(2))
STRESS(3)=(DDSDDE(3,1)*DSTRAN(1))+(DDSDDE(3,2)*DSTRAN(2))
1 +(DDSDDE(3,3)*DSTRAN(3))
STRESS(4)=(DDSDDE(4,4)*DSTRAN(4))
STRESS(5)=(DDSDDE(5,5)*DSTRAN(5))
STRESS(6)=(DDSDDE(6,6)*DSTRAN(6))
END IF
END IF
C HASHIN FIBRE
IF (STRESS(1) .GE. 0) THEN
uFTF=(STRESS(1)/XT)**2
END IF
IF (STRESS(1) .LT. 0) THEN
uFCF=(STRESS(1)/XC)**2
END IF
C HASHIN MATRICE
IF (STRESS(2) .GE. 0 .AND. NTENS .EQ. 3) THEN
uFTM=((STRESS(2)/YT)**2)+((STRESS(3)/SL)**2)
IF (STRESS(2) .GE. 0 .AND. NTENS .EQ. 6 ) THEN
uFTM=((STRESS(2)/YT)**2)+((STRESS(6)/SL)**2)
END IF
END IF
IF (STRESS(2) .LT. 0) THEN
uFCM=((STRESS(2)/(2*ST))**2)+((((YC/(2*ST))**2)-1)*(STRESS(2)/YC))+
* (STRESS(3)/SL)**2
END IF
STATEV(1) = STRESS(1)
STATEV(2) = STRESS(2)
STATEV(3) = STRESS(3)
STATEV(4) = STRESS(5)
STATEV(5) = STRESS(6)
STATEV(6) = DDSDDE(1,1)
STATEV(7) = DDSDDE(1,2)
STATEV(8) = DDSDDE(1,3)
STATEV(9) = DDSDDE(2,1)
STATEV(10) = DDSDDE(2,2)
STATEV(11) = DDSDDE(2,3)
STATEV(12) = DDSDDE(3,1)
STATEV(13) = DDSDDE(3,2)
STATEV(14) = DDSDDE(3,3)
STATEV(15) = DDSDDE(4,4)
STATEV(16) = DDSDDE(5,5)
STATEV(17) = DDSDDE(6,6)
STATEV(18) = DSTRAN(1)
STATEV(19) = DSTRAN(2)
STATEV(20) = DSTRAN(3)
STATEV(21) = DSTRAN(4)
STATEV(22) = DSTRAN(5)
STATEV(23) = DSTRAN(6)
C-----------------------------------------------------------------------------------------------------
RETURN
END
The STRESSES (S11,S22,S12) and the STRAINS (E11,E22,E12) are not the same in the both model, why ?
How could you explain that ?
Please help me !