-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCantileverToy_Hexahedra.m
executable file
·624 lines (593 loc) · 19.3 KB
/
CantileverToy_Hexahedra.m
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
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
%==================================================================
% General Data File
% Title: HEXAHEDRA
% Units: SI
% Dimensions: 3D
% Type of problem: Plane_Stress
% Type of Phisics: ELASTIC
% Micro/Macro: MACRO
%
%==================================================================
%% Data
Data_prb = {
'HEXAHEDRA';
'SI';
'3D';
'Plane_Stress';
'ELASTIC';
'MACRO';
};
%% Coordinates
% Node X Y Z
coord = [
1 12 0 2
2 12 0.5 1.5
3 12 -0.5 1.5
4 12 4.4409e-16 1
5 10.8 0 2
6 10.8 -0.5 1.5
7 10.8 0.5 1.5
8 12 -1 1
9 12 1 1
10 10.8 4.7184e-16 1
11 12 0.5 0.5
12 12 -0.5 0.5
13 10.8 1 1
14 10.8 -1 1
15 10.8 0.5 0.5
16 10.8 -0.5 0.5
17 12 0 0
18 12 1.5 0.5
19 12 -1.5 0.5
20 12 -1 0
21 12 1 0
22 10.8 0 0
23 9.6 0 2
24 10.8 1.5 0.5
25 10.8 -1.5 0.5
26 9.6 -0.5 1.5
27 9.6 0.5 1.5
28 10.8 -1 -2.7756e-17
29 10.8 1 0
30 12 0.5 -0.5
31 12 -0.5 -0.5
32 9.6 4.4409e-16 1
33 9.6 1 1
34 9.6 -1 1
35 10.8 0.5 -0.5
36 10.8 -0.5 -0.5
37 12 -2 0
38 12 2 0
39 9.6 0.5 0.5
40 9.6 -0.5 0.5
41 12 -1.5 -0.5
42 12 1.5 -0.5
43 12 4.4409e-16 -1
44 10.8 2 0
45 10.8 -2 0
46 9.6 0 0
47 10.8 1.5 -0.5
48 10.8 -1.5 -0.5
49 12 -1 -1
50 12 1 -1
51 9.6 1.5 0.5
52 9.6 -1.5 0.5
53 10.8 4.4409e-16 -1
54 9.6 -1 0
55 9.6 1 -2.7756e-17
56 10.8 -1 -1
57 10.8 1 -1
58 9.6 0.5 -0.5
59 9.6 -0.5 -0.5
60 12 -0.5 -1.5
61 12 0.5 -1.5
62 8.4 0 2
63 8.4 -0.5 1.5
64 8.4 0.5 1.5
65 9.6 -2 0
66 9.6 2 0
67 10.8 0.5 -1.5
68 10.8 -0.5 -1.5
69 8.4 4.4409e-16 1
70 9.6 -1.5 -0.5
71 9.6 1.5 -0.5
72 9.6 4.7184e-16 -1
73 8.4 -1 1
74 8.4 1 1
75 8.4 -0.5 0.5
76 8.4 0.5 0.5
77 9.6 -1 -1
78 9.6 1 -1
79 12 0 -2
80 8.4 0 0
81 10.8 0 -2
82 8.4 -1.5 0.5
83 8.4 1.5 0.5
84 8.4 -1 -2.7756e-17
85 8.4 1 0
86 9.6 -0.5 -1.5
87 9.6 0.5 -1.5
88 8.4 -0.5 -0.5
89 8.4 0.5 -0.5
90 8.4 2 0
91 8.4 -2 0
92 8.4 -1.5 -0.5
93 8.4 1.5 -0.5
94 9.6 0 -2
95 8.4 4.4409e-16 -1
96 8.4 1 -1
97 8.4 -1 -1
98 7.2 0 2
99 7.2 0.5 1.5
100 7.2 -0.5 1.5
101 7.2 4.4409e-16 1
102 7.2 1 1
103 7.2 -1 1
104 8.4 -0.5 -1.5
105 8.4 0.5 -1.5
106 7.2 -0.5 0.5
107 7.2 0.5 0.5
108 7.2 -5.5511e-17 5.5511e-17
109 7.2 -1.5 0.5
110 7.2 1.5 0.5
111 7.2 1 -5.5511e-17
112 7.2 -1 5.5511e-17
113 8.4 0 -2
114 7.2 -0.5 -0.5
115 7.2 0.5 -0.5
116 7.2 -2 0
117 7.2 2 0
118 7.2 1.5 -0.5
119 7.2 -1.5 -0.5
120 7.2 4.996e-16 -1
121 7.2 -1 -1
122 7.2 1 -1
123 7.2 -0.5 -1.5
124 7.2 0.5 -1.5
125 6 0 2
126 6 0.5 1.5
127 6 -0.5 1.5
128 6 4.4409e-16 1
129 6 -1 1
130 6 1 1
131 6 0.5 0.5
132 6 -0.5 0.5
133 7.2 0 -2
134 6 0 0
135 6 1.5 0.5
136 6 -1.5 0.5
137 6 1 0
138 6 -1 2.7756e-17
139 6 0.5 -0.5
140 6 -0.5 -0.5
141 6 2 0
142 6 -2 0
143 6 1.5 -0.5
144 6 -1.5 -0.5
145 6 4.4409e-16 -1
146 6 -1 -1
147 6 1 -1
148 6 -0.5 -1.5
149 6 0.5 -1.5
150 4.8 0 2
151 6 0 -2
152 4.8 0.5 1.5
153 4.8 -0.5 1.5
154 4.8 4.4409e-16 1
155 4.8 1 1
156 4.8 -1 1
157 4.8 -0.5 0.5
158 4.8 0.5 0.5
159 4.8 -2.7756e-17 2.7756e-17
160 4.8 1.5 0.5
161 4.8 -1.5 0.5
162 4.8 1 0
163 4.8 -1 0
164 4.8 0.5 -0.5
165 4.8 -0.5 -0.5
166 4.8 -2 0
167 4.8 2 0
168 4.8 1.5 -0.5
169 4.8 -1.5 -0.5
170 4.8 4.4409e-16 -1
171 4.8 -1 -1
172 4.8 1 -1
173 4.8 0.5 -1.5
174 4.8 -0.5 -1.5
175 4.8 0 -2
176 3.6 0 2
177 3.6 -0.5 1.5
178 3.6 0.5 1.5
179 3.6 4.4409e-16 1
180 3.6 -1 1
181 3.6 1 1
182 3.6 -0.5 0.5
183 3.6 0.5 0.5
184 3.6 -1.3878e-17 1.3878e-17
185 3.6 1.5 0.5
186 3.6 -1.5 0.5
187 3.6 1 0
188 3.6 -1 0
189 3.6 0.5 -0.5
190 3.6 -0.5 -0.5
191 3.6 2 0
192 3.6 -2 0
193 3.6 -1.5 -0.5
194 3.6 1.5 -0.5
195 3.6 4.996e-16 -1
196 3.6 1 -1
197 3.6 -1 -1
198 3.6 -0.5 -1.5
199 3.6 0.5 -1.5
200 3.6 0 -2
201 2.4 0 2
202 2.4 -0.5 1.5
203 2.4 0.5 1.5
204 2.4 4.7184e-16 1
205 2.4 -1 1
206 2.4 1 1
207 2.4 -0.5 0.5
208 2.4 0.5 0.5
209 2.4 -6.9389e-18 6.9389e-18
210 2.4 1.5 0.5
211 2.4 -1.5 0.5
212 2.4 1 -2.7756e-17
213 2.4 -1 -2.7756e-17
214 2.4 0.5 -0.5
215 2.4 -0.5 -0.5
216 2.4 -2 0
217 2.4 2 0
218 2.4 1.5 -0.5
219 2.4 -1.5 -0.5
220 2.4 4.1633e-16 -1
221 2.4 1 -1
222 2.4 -1 -1
223 2.4 -0.5 -1.5
224 2.4 0.5 -1.5
225 2.4 0 -2
226 1.2 0 2
227 1.2 0.5 1.5
228 1.2 -0.5 1.5
229 1.2 4.7184e-16 1
230 1.2 -1 1
231 1.2 1 1
232 1.2 0.5 0.5
233 1.2 -0.5 0.5
234 1.2 -1.7347e-18 1.7347e-18
235 1.2 1.5 0.5
236 1.2 -1.5 0.5
237 1.2 -1 -2.7756e-17
238 1.2 1 0
239 1.2 -0.5 -0.5
240 1.2 0.5 -0.5
241 1.2 -2 0
242 1.2 2 0
243 1.2 1.5 -0.5
244 1.2 -1.5 -0.5
245 1.2 4.4409e-16 -1
246 1.2 1 -1
247 1.2 -1 -1
248 1.2 0.5 -1.5
249 1.2 -0.5 -1.5
250 1.2 0 -2
251 0 0 2
252 0 -0.5 1.5
253 0 0.5 1.5
254 0 4.4409e-16 1
255 0 1 1
256 0 -1 1
257 0 -0.5 0.5
258 0 0.5 0.5
259 0 0 0
260 0 -1.5 0.5
261 0 1.5 0.5
262 0 1 0
263 0 -1 0
264 0 -0.5 -0.5
265 0 0.5 -0.5
266 0 2 0
267 0 -2 0
268 0 -1.5 -0.5
269 0 1.5 -0.5
270 0 4.4409e-16 -1
271 0 1 -1
272 0 -1 -1
273 0 -0.5 -1.5
274 0 0.5 -1.5
275 0 0 -2
];
%% Conectivities
% Element Node(1) Node(2) Node(3) Node(4) Node(5) Node(6) Node(7) Node(8) Material
connec = [
1 237 236 260 263 244 241 267 268 0
2 239 237 263 264 247 244 268 272 0
3 245 239 264 270 249 247 272 273 0
4 248 245 270 274 250 249 273 275 0
5 213 211 236 237 219 216 241 244 0
6 215 213 237 239 222 219 244 247 0
7 220 215 239 245 223 222 247 249 0
8 224 220 245 248 225 223 249 250 0
9 188 186 211 213 193 192 216 219 0
10 190 188 213 215 197 193 219 222 0
11 195 190 215 220 198 197 222 223 0
12 199 195 220 224 200 198 223 225 0
13 163 161 186 188 169 166 192 193 0
14 165 163 188 190 171 169 193 197 0
15 170 165 190 195 174 171 197 198 0
16 173 170 195 199 175 174 198 200 0
17 138 136 161 163 144 142 166 169 0
18 140 138 163 165 146 144 169 171 0
19 145 140 165 170 148 146 171 174 0
20 149 145 170 173 151 148 174 175 0
21 112 109 136 138 119 116 142 144 0
22 114 112 138 140 121 119 144 146 0
23 120 114 140 145 123 121 146 148 0
24 124 120 145 149 133 123 148 151 0
25 84 82 109 112 92 91 116 119 0
26 88 84 112 114 97 92 119 121 0
27 95 88 114 120 104 97 121 123 0
28 105 95 120 124 113 104 123 133 0
29 54 52 82 84 70 65 91 92 0
30 59 54 84 88 77 70 92 97 0
31 72 59 88 95 86 77 97 104 0
32 87 72 95 105 94 86 104 113 0
33 28 25 52 54 48 45 65 70 0
34 36 28 54 59 56 48 70 77 0
35 53 36 59 72 68 56 77 86 0
36 67 53 72 87 81 68 86 94 0
37 20 19 25 28 41 37 45 48 0
38 31 20 28 36 49 41 48 56 0
39 43 31 36 53 60 49 56 68 0
40 61 43 53 67 79 60 68 81 0
41 233 230 256 257 237 236 260 263 0
42 234 233 257 259 239 237 263 264 0
43 240 234 259 265 245 239 264 270 0
44 246 240 265 271 248 245 270 274 0
45 207 205 230 233 213 211 236 237 0
46 209 207 233 234 215 213 237 239 0
47 214 209 234 240 220 215 239 245 0
48 221 214 240 246 224 220 245 248 0
49 182 180 205 207 188 186 211 213 0
50 184 182 207 209 190 188 213 215 0
51 189 184 209 214 195 190 215 220 0
52 196 189 214 221 199 195 220 224 0
53 157 156 180 182 163 161 186 188 0
54 159 157 182 184 165 163 188 190 0
55 164 159 184 189 170 165 190 195 0
56 172 164 189 196 173 170 195 199 0
57 132 129 156 157 138 136 161 163 0
58 134 132 157 159 140 138 163 165 0
59 139 134 159 164 145 140 165 170 0
60 147 139 164 172 149 145 170 173 0
61 106 103 129 132 112 109 136 138 0
62 108 106 132 134 114 112 138 140 0
63 115 108 134 139 120 114 140 145 0
64 122 115 139 147 124 120 145 149 0
65 75 73 103 106 84 82 109 112 0
66 80 75 106 108 88 84 112 114 0
67 89 80 108 115 95 88 114 120 0
68 96 89 115 122 105 95 120 124 0
69 40 34 73 75 54 52 82 84 0
70 46 40 75 80 59 54 84 88 0
71 58 46 80 89 72 59 88 95 0
72 78 58 89 96 87 72 95 105 0
73 16 14 34 40 28 25 52 54 0
74 22 16 40 46 36 28 54 59 0
75 35 22 46 58 53 36 59 72 0
76 57 35 58 78 67 53 72 87 0
77 12 8 14 16 20 19 25 28 0
78 17 12 16 22 31 20 28 36 0
79 30 17 22 35 43 31 36 53 0
80 50 30 35 57 61 43 53 67 0
81 229 228 252 254 233 230 256 257 0
82 232 229 254 258 234 233 257 259 0
83 238 232 258 262 240 234 259 265 0
84 243 238 262 269 246 240 265 271 0
85 204 202 228 229 207 205 230 233 0
86 208 204 229 232 209 207 233 234 0
87 212 208 232 238 214 209 234 240 0
88 218 212 238 243 221 214 240 246 0
89 179 177 202 204 182 180 205 207 0
90 183 179 204 208 184 182 207 209 0
91 187 183 208 212 189 184 209 214 0
92 194 187 212 218 196 189 214 221 0
93 154 153 177 179 157 156 180 182 0
94 158 154 179 183 159 157 182 184 0
95 162 158 183 187 164 159 184 189 0
96 168 162 187 194 172 164 189 196 0
97 128 127 153 154 132 129 156 157 0
98 131 128 154 158 134 132 157 159 0
99 137 131 158 162 139 134 159 164 0
100 143 137 162 168 147 139 164 172 0
101 101 100 127 128 106 103 129 132 0
102 107 101 128 131 108 106 132 134 0
103 111 107 131 137 115 108 134 139 0
104 118 111 137 143 122 115 139 147 0
105 69 63 100 101 75 73 103 106 0
106 76 69 101 107 80 75 106 108 0
107 85 76 107 111 89 80 108 115 0
108 93 85 111 118 96 89 115 122 0
109 32 26 63 69 40 34 73 75 0
110 39 32 69 76 46 40 75 80 0
111 55 39 76 85 58 46 80 89 0
112 71 55 85 93 78 58 89 96 0
113 10 6 26 32 16 14 34 40 0
114 15 10 32 39 22 16 40 46 0
115 29 15 39 55 35 22 46 58 0
116 47 29 55 71 57 35 58 78 0
117 4 3 6 10 12 8 14 16 0
118 11 4 10 15 17 12 16 22 0
119 21 11 15 29 30 17 22 35 0
120 42 21 29 47 50 30 35 57 0
121 227 226 251 253 229 228 252 254 0
122 231 227 253 255 232 229 254 258 0
123 235 231 255 261 238 232 258 262 0
124 242 235 261 266 243 238 262 269 0
125 203 201 226 227 204 202 228 229 0
126 206 203 227 231 208 204 229 232 0
127 210 206 231 235 212 208 232 238 0
128 217 210 235 242 218 212 238 243 0
129 178 176 201 203 179 177 202 204 0
130 181 178 203 206 183 179 204 208 0
131 185 181 206 210 187 183 208 212 0
132 191 185 210 217 194 187 212 218 0
133 152 150 176 178 154 153 177 179 0
134 155 152 178 181 158 154 179 183 0
135 160 155 181 185 162 158 183 187 0
136 167 160 185 191 168 162 187 194 0
137 126 125 150 152 128 127 153 154 0
138 130 126 152 155 131 128 154 158 0
139 135 130 155 160 137 131 158 162 0
140 141 135 160 167 143 137 162 168 0
141 99 98 125 126 101 100 127 128 0
142 102 99 126 130 107 101 128 131 0
143 110 102 130 135 111 107 131 137 0
144 117 110 135 141 118 111 137 143 0
145 64 62 98 99 69 63 100 101 0
146 74 64 99 102 76 69 101 107 0
147 83 74 102 110 85 76 107 111 0
148 90 83 110 117 93 85 111 118 0
149 27 23 62 64 32 26 63 69 0
150 33 27 64 74 39 32 69 76 0
151 51 33 74 83 55 39 76 85 0
152 66 51 83 90 71 55 85 93 0
153 7 5 23 27 10 6 26 32 0
154 13 7 27 33 15 10 32 39 0
155 24 13 33 51 29 15 39 55 0
156 44 24 51 66 47 29 55 71 0
157 2 1 5 7 4 3 6 10 0
158 9 2 7 13 11 4 10 15 0
159 18 9 13 24 21 11 15 29 0
160 38 18 24 44 42 21 29 47 0
];
%% Variable Prescribed
% Node Dimension Value
dirichlet_data = [
251 1 0
251 2 0
251 3 0
252 1 0
252 2 0
252 3 0
253 1 0
253 2 0
253 3 0
254 1 0
254 2 0
254 3 0
255 1 0
255 2 0
255 3 0
256 1 0
256 2 0
256 3 0
257 1 0
257 2 0
257 3 0
258 1 0
258 2 0
258 3 0
259 1 0
259 2 0
259 3 0
260 1 0
260 2 0
260 3 0
261 1 0
261 2 0
261 3 0
262 1 0
262 2 0
262 3 0
263 1 0
263 2 0
263 3 0
264 1 0
264 2 0
264 3 0
265 1 0
265 2 0
265 3 0
266 1 0
266 2 0
266 3 0
267 1 0
267 2 0
267 3 0
268 1 0
268 2 0
268 3 0
269 1 0
269 2 0
269 3 0
270 1 0
270 2 0
270 3 0
271 1 0
271 2 0
271 3 0
272 1 0
272 2 0
272 3 0
273 1 0
273 2 0
273 3 0
274 1 0
274 2 0
274 3 0
275 1 0
275 2 0
275 3 0
];
%% Force Prescribed
% Node Dimension Value
pointload_complete = [
17 2 -1
];
%% Volumetric Force
% Element Dim Force_Dim
Vol_force = [
];
%% Group Elements
% Element Group_num
Group = [
];
%% Initial Holes
% Elements that are considered holes initially
% Element
Initial_holes = [
];
%% Boundary Elements
% Elements that can not be removed
% Element
Boundary_elements = [
];
%% Micro gauss post
%
% Element
Micro_gauss_post = [
];
%% Micro Slave-Master
% Nodes that are Slaves
% Nodes Value (1-Slave,0-Master)
Micro_slave = [
];
%% Nodes solid
% Nodes that must remain
% Nodes
nodesolid = unique(pointload_complete(:,1));
%% External border Elements
% Detect the elements that define the edge of the domain
% Element Node(1) Node(2)
External_border_elements = [
];
%% External border Nodes
% Detect the nodes that define the edge of the domain
% Node
External_border_nodes = [
];
%% Materials
% Materials that have been used
% Material_Num Mat_density Young_Modulus Poisson
Materials = [
];