-
Notifications
You must be signed in to change notification settings - Fork 0
/
Bridge_Hexahedra_Linear_Structured.m
executable file
·673 lines (642 loc) · 22.6 KB
/
Bridge_Hexahedra_Linear_Structured.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
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
%==================================================================
% 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
gidcoord = [
1 0 20 20
2 5 20 20
3 0 15 20
4 0 20 15
5 5 20 15
6 5 15 20
7 0 15 15
8 5 15 15
9 0 20 10
10 10 20 20
11 0 10 20
12 10 20 15
13 5 20 10
14 0 15 10
15 0 10 15
16 10 15 20
17 5 10 20
18 10 15 15
19 5 10 15
20 5 15 10
21 10 20 10
22 0 10 10
23 10 10 20
24 10 10 15
25 0 20 5
26 15 20 20
27 10 15 10
28 0 5 20
29 5 10 10
30 15 15 20
31 5 5 20
32 5 20 5
33 0 15 5
34 0 5 15
35 15 20 15
36 15 15 15
37 5 5 15
38 5 15 5
39 10 10 10
40 10 20 5
41 0 5 10
42 10 5 20
43 15 10 20
44 0 10 5
45 15 20 10
46 15 15 10
47 5 10 5
48 15 10 15
49 10 5 15
50 10 15 5
51 5 5 10
52 20 20 20
53 0 20 0
54 0 0 20
55 10 5 10
56 10 10 5
57 15 10 10
58 20 20 15
59 5 0 20
60 0 0 15
61 5 20 0
62 0 15 0
63 20 15 20
64 15 5 20
65 0 5 5
66 5 0 15
67 5 15 0
68 20 15 15
69 15 20 5
70 5 5 5
71 15 15 5
72 15 5 15
73 20 20 10
74 0 0 10
75 10 0 20
76 0 10 0
77 10 20 0
78 20 10 20
79 10 0 15
80 5 0 10
81 20 15 10
82 5 10 0
83 10 15 0
84 20 10 15
85 10 5 5
86 15 5 10
87 15 10 5
88 10 10 0
89 10 0 10
90 20 10 10
91 20 20 5
92 25 20 20
93 15 0 20
94 0 0 5
95 0 5 0
96 15 20 0
97 20 5 20
98 15 15 0
99 5 5 0
100 20 15 5
101 5 0 5
102 25 20 15
103 15 0 15
104 20 5 15
105 25 15 20
106 15 5 5
107 25 15 15
108 25 10 20
109 10 0 5
110 15 10 0
111 10 5 0
112 20 5 10
113 20 10 5
114 25 20 10
115 15 0 10
116 25 10 15
117 25 15 10
118 20 0 20
119 20 20 0
120 0 0 0
121 25 10 10
122 5 0 0
123 20 15 0
124 20 0 15
125 25 5 20
126 25 20 5
127 15 5 0
128 20 5 5
129 15 0 5
130 25 5 15
131 25 15 5
132 30 20 20
133 10 0 0
134 20 10 0
135 20 0 10
136 30 15 20
137 30 20 15
138 25 5 10
139 25 10 5
140 30 15 15
141 30 10 20
142 30 20 10
143 30 15 10
144 30 10 15
145 25 0 20
146 15 0 0
147 20 0 5
148 25 20 0
149 20 5 0
150 25 0 15
151 25 15 0
152 25 5 5
153 30 10 10
154 30 5 20
155 25 10 0
156 30 20 5
157 25 0 10
158 30 5 15
159 30 15 5
160 20 0 0
161 30 10 5
162 30 5 10
163 35 20 20
164 25 5 0
165 35 15 20
166 25 0 5
167 35 20 15
168 35 15 15
169 30 20 0
170 30 20 0
171 30 0 20
172 30 16.667 0
173 35 10 20
174 30 15 0
175 35 20 10
176 30 0 15
177 30 13.333 0
178 35 10 15
179 30 5 5
180 35 15 10
181 30 10 0
182 30 0 10
183 30 10 0
184 35 10 10
185 25 0 0
186 35 5 20
187 35 20 5
188 35 5 15
189 35 15 5
190 30 6.6667 0
191 30 5 0
192 30 0 5
193 35 5 10
194 35 10 5
195 30 3.3333 0
196 40 20 20
197 40 15 20
198 40 20 15
199 35 0 20
200 35 20 0
201 40 15 15
202 35 15 0
203 35 0 15
204 35 5 5
205 40 20 10
206 40 10 20
207 30 0 0
208 30 0 0
209 40 10 15
210 40 15 10
211 35 10 0
212 35 0 10
213 40 10 10
214 40 5 20
215 40 20 5
216 40 5 15
217 40 15 5
218 35 5 0
219 35 0 5
220 40 10 5
221 40 5 10
222 40 20 0
223 40 0 20
224 40 15 0
225 35 0 0
226 40 0 15
227 45 20 20
228 45 15 20
229 40 5 5
230 45 20 15
231 45 15 15
232 40 10 0
233 40 0 10
234 45 10 20
235 45 20 10
236 45 15 10
237 45 10 15
238 40 5 0
239 45 10 10
240 40 0 5
241 45 5 20
242 45 20 5
243 45 15 5
244 45 5 15
245 45 10 5
246 45 5 10
247 40 0 0
248 45 20 0
249 45 0 20
250 45 15 0
251 45 0 15
252 45 5 5
253 50 20 20
254 50 15 20
255 45 10 0
256 50 20 15
257 45 0 10
258 50 15 15
259 50 20 10
260 50 10 20
261 50 15 10
262 50 10 15
263 45 5 0
264 45 0 5
265 50 10 10
266 50 20 5
267 50 5 20
268 50 5 15
269 50 15 5
270 50 10 5
271 50 5 10
272 45 0 0
273 50 0 20
274 50 20 0
275 50 15 0
276 50 0 15
277 50 5 5
278 50 10 0
279 50 0 10
280 55 20 20
281 55 20 15
282 55 15 20
283 55 15 15
284 50 5 0
285 55 10 20
286 55 20 10
287 50 0 5
288 55 10 15
289 55 15 10
290 55 10 10
291 55 20 5
292 55 5 20
293 55 5 15
294 55 15 5
295 50 0 0
296 55 10 5
297 55 5 10
298 55 20 0
299 55 0 20
300 55 0 15
301 55 15 0
302 55 5 5
303 55 10 0
304 55 0 10
305 60 20 20
306 60 15 20
307 60 20 15
308 60 15 15
309 55 0 5
310 55 5 0
311 60 20 10
312 60 10 20
313 60 15 10
314 60 10 15
315 60 10 10
316 60 5 20
317 55 0 0
318 60 20 5
319 60 5 15
320 60 15 5
321 60 5 10
322 60 10 5
323 60 0 20
324 60 20 0
325 60 0 15
326 60 15 0
327 60 5 5
328 60 0 10
329 60 10 0
330 60 5 0
331 60 0 5
332 60 0 0
];
%% Conectivities
% Element Node(1) Node(2) Node(3) Node(4) Node(5) Node(6) Node(7) Node(8) Material
gidlnods = [
1 70 65 95 99 101 94 120 122 0
2 85 70 99 111 109 101 122 133 0
3 106 85 111 127 129 109 133 146 0
4 128 106 127 149 147 129 146 160 0
5 152 128 149 164 166 147 160 185 0
6 179 152 164 191 192 166 185 208 0
7 204 179 191 218 219 192 208 225 0
8 229 204 218 238 240 219 225 247 0
9 252 229 238 263 264 240 247 272 0
10 277 252 263 284 287 264 272 295 0
11 302 277 284 310 309 287 295 317 0
12 327 302 310 330 331 309 317 332 0
13 51 41 65 70 80 74 94 101 0
14 55 51 70 85 89 80 101 109 0
15 86 55 85 106 115 89 109 129 0
16 112 86 106 128 135 115 129 147 0
17 138 112 128 152 157 135 147 166 0
18 162 138 152 179 182 157 166 192 0
19 193 162 179 204 212 182 192 219 0
20 221 193 204 229 233 212 219 240 0
21 246 221 229 252 257 233 240 264 0
22 271 246 252 277 279 257 264 287 0
23 297 271 277 302 304 279 287 309 0
24 321 297 302 327 328 304 309 331 0
25 37 34 41 51 66 60 74 80 0
26 49 37 51 55 79 66 80 89 0
27 72 49 55 86 103 79 89 115 0
28 104 72 86 112 124 103 115 135 0
29 130 104 112 138 150 124 135 157 0
30 158 130 138 162 176 150 157 182 0
31 188 158 162 193 203 176 182 212 0
32 216 188 193 221 226 203 212 233 0
33 244 216 221 246 251 226 233 257 0
34 268 244 246 271 276 251 257 279 0
35 293 268 271 297 300 276 279 304 0
36 319 293 297 321 325 300 304 328 0
37 31 28 34 37 59 54 60 66 0
38 42 31 37 49 75 59 66 79 0
39 64 42 49 72 93 75 79 103 0
40 97 64 72 104 118 93 103 124 0
41 125 97 104 130 145 118 124 150 0
42 154 125 130 158 171 145 150 176 0
43 186 154 158 188 199 171 176 203 0
44 214 186 188 216 223 199 203 226 0
45 241 214 216 244 249 223 226 251 0
46 267 241 244 268 273 249 251 276 0
47 292 267 268 293 299 273 276 300 0
48 316 292 293 319 323 299 300 325 0
49 47 44 76 82 70 65 95 99 0
50 56 47 82 88 85 70 99 111 0
51 87 56 88 110 106 85 111 127 0
52 113 87 110 134 128 106 127 149 0
53 139 113 134 155 152 128 149 164 0
54 161 139 155 181 179 152 164 191 0
55 194 161 181 211 204 179 191 218 0
56 220 194 211 232 229 204 218 238 0
57 245 220 232 255 252 229 238 263 0
58 270 245 255 278 277 252 263 284 0
59 296 270 278 303 302 277 284 310 0
60 322 296 303 329 327 302 310 330 0
61 29 22 44 47 51 41 65 70 0
62 39 29 47 56 55 51 70 85 0
63 57 39 56 87 86 55 85 106 0
64 90 57 87 113 112 86 106 128 0
65 121 90 113 139 138 112 128 152 0
66 153 121 139 161 162 138 152 179 0
67 184 153 161 194 193 162 179 204 0
68 213 184 194 220 221 193 204 229 0
69 239 213 220 245 246 221 229 252 0
70 265 239 245 270 271 246 252 277 0
71 290 265 270 296 297 271 277 302 0
72 315 290 296 322 321 297 302 327 0
73 19 15 22 29 37 34 41 51 0
74 24 19 29 39 49 37 51 55 0
75 48 24 39 57 72 49 55 86 0
76 84 48 57 90 104 72 86 112 0
77 116 84 90 121 130 104 112 138 0
78 144 116 121 153 158 130 138 162 0
79 178 144 153 184 188 158 162 193 0
80 209 178 184 213 216 188 193 221 0
81 237 209 213 239 244 216 221 246 0
82 262 237 239 265 268 244 246 271 0
83 288 262 265 290 293 268 271 297 0
84 314 288 290 315 319 293 297 321 0
85 17 11 15 19 31 28 34 37 0
86 23 17 19 24 42 31 37 49 0
87 43 23 24 48 64 42 49 72 0
88 78 43 48 84 97 64 72 104 0
89 108 78 84 116 125 97 104 130 0
90 141 108 116 144 154 125 130 158 0
91 173 141 144 178 186 154 158 188 0
92 206 173 178 209 214 186 188 216 0
93 234 206 209 237 241 214 216 244 0
94 260 234 237 262 267 241 244 268 0
95 285 260 262 288 292 267 268 293 0
96 312 285 288 314 316 292 293 319 0
97 38 33 62 67 47 44 76 82 0
98 50 38 67 83 56 47 82 88 0
99 71 50 83 98 87 56 88 110 0
100 100 71 98 123 113 87 110 134 0
101 131 100 123 151 139 113 134 155 0
102 159 131 151 174 161 139 155 181 0
103 189 159 174 202 194 161 181 211 0
104 217 189 202 224 220 194 211 232 0
105 243 217 224 250 245 220 232 255 0
106 269 243 250 275 270 245 255 278 0
107 294 269 275 301 296 270 278 303 0
108 320 294 301 326 322 296 303 329 0
109 20 14 33 38 29 22 44 47 0
110 27 20 38 50 39 29 47 56 0
111 46 27 50 71 57 39 56 87 0
112 81 46 71 100 90 57 87 113 0
113 117 81 100 131 121 90 113 139 0
114 143 117 131 159 153 121 139 161 0
115 180 143 159 189 184 153 161 194 0
116 210 180 189 217 213 184 194 220 0
117 236 210 217 243 239 213 220 245 0
118 261 236 243 269 265 239 245 270 0
119 289 261 269 294 290 265 270 296 0
120 313 289 294 320 315 290 296 322 0
121 8 7 14 20 19 15 22 29 0
122 18 8 20 27 24 19 29 39 0
123 36 18 27 46 48 24 39 57 0
124 68 36 46 81 84 48 57 90 0
125 107 68 81 117 116 84 90 121 0
126 140 107 117 143 144 116 121 153 0
127 168 140 143 180 178 144 153 184 0
128 201 168 180 210 209 178 184 213 0
129 231 201 210 236 237 209 213 239 0
130 258 231 236 261 262 237 239 265 0
131 283 258 261 289 288 262 265 290 0
132 308 283 289 313 314 288 290 315 0
133 6 3 7 8 17 11 15 19 0
134 16 6 8 18 23 17 19 24 0
135 30 16 18 36 43 23 24 48 0
136 63 30 36 68 78 43 48 84 0
137 105 63 68 107 108 78 84 116 0
138 136 105 107 140 141 108 116 144 0
139 165 136 140 168 173 141 144 178 0
140 197 165 168 201 206 173 178 209 0
141 228 197 201 231 234 206 209 237 0
142 254 228 231 258 260 234 237 262 0
143 282 254 258 283 285 260 262 288 0
144 306 282 283 308 312 285 288 314 0
145 32 25 53 61 38 33 62 67 0
146 40 32 61 77 50 38 67 83 0
147 69 40 77 96 71 50 83 98 0
148 91 69 96 119 100 71 98 123 0
149 126 91 119 148 131 100 123 151 0
150 156 126 148 169 159 131 151 174 0
151 187 156 169 200 189 159 174 202 0
152 215 187 200 222 217 189 202 224 0
153 242 215 222 248 243 217 224 250 0
154 266 242 248 274 269 243 250 275 0
155 291 266 274 298 294 269 275 301 0
156 318 291 298 324 320 294 301 326 0
157 13 9 25 32 20 14 33 38 0
158 21 13 32 40 27 20 38 50 0
159 45 21 40 69 46 27 50 71 0
160 73 45 69 91 81 46 71 100 0
161 114 73 91 126 117 81 100 131 0
162 142 114 126 156 143 117 131 159 0
163 175 142 156 187 180 143 159 189 0
164 205 175 187 215 210 180 189 217 0
165 235 205 215 242 236 210 217 243 0
166 259 235 242 266 261 236 243 269 0
167 286 259 266 291 289 261 269 294 0
168 311 286 291 318 313 289 294 320 0
169 5 4 9 13 8 7 14 20 0
170 12 5 13 21 18 8 20 27 0
171 35 12 21 45 36 18 27 46 0
172 58 35 45 73 68 36 46 81 0
173 102 58 73 114 107 68 81 117 0
174 137 102 114 142 140 107 117 143 0
175 167 137 142 175 168 140 143 180 0
176 198 167 175 205 201 168 180 210 0
177 230 198 205 235 231 201 210 236 0
178 256 230 235 259 258 231 236 261 0
179 281 256 259 286 283 258 261 289 0
180 307 281 286 311 308 283 289 313 0
181 2 1 4 5 6 3 7 8 0
182 10 2 5 12 16 6 8 18 0
183 26 10 12 35 30 16 18 36 0
184 52 26 35 58 63 30 36 68 0
185 92 52 58 102 105 63 68 107 0
186 132 92 102 137 136 105 107 140 0
187 163 132 137 167 165 136 140 168 0
188 196 163 167 198 197 165 168 201 0
189 227 196 198 230 228 197 201 231 0
190 253 227 230 256 254 228 231 258 0
191 280 253 256 281 282 254 258 283 0
192 305 280 281 307 306 282 283 308 0
193 207 195 0
194 195 190 0
195 190 183 0
196 183 177 0
197 177 172 0
198 172 170 0
];
%% Variable Prescribed
% Node Dimension Value
lnodes = [
53 1 0
53 2 0
53 3 0
62 1 0
62 2 0
62 3 0
76 1 0
76 2 0
76 3 0
95 1 0
95 2 0
95 3 0
120 1 0
120 2 0
120 3 0
324 1 0
324 2 0
324 3 0
326 1 0
326 2 0
326 3 0
329 1 0
329 2 0
329 3 0
330 1 0
330 2 0
330 3 0
332 1 0
332 2 0
332 3 0
];
%% Force Prescribed
% Node Dimension Value
pointload_complete = [
];
%% 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 = [
];