CantileverToy_Quadrilater.m

Data_prb = {
'QUAD';
'SI';
'2D';
'Plane_Stress';
'ELASTIC';
'MACRO';
};
coord=[
 1              12               0               0
    2              12             0.2               0
    3              12             0.4               0
    4              12             0.6               0
    5              12             0.8               0
    6              12               1               0
    7            10.8               0               0
    8              12             1.2               0
    9            10.8             0.2               0
   10            10.8             0.4               0
   11            10.8             0.6               0
   12              12             1.4               0
   13            10.8             0.8               0
   14            10.8               1               0
   15              12             1.6               0
   16            10.8             1.2               0
   17              12             1.8               0
   18            10.8             1.4               0
   19              12               2               0
   20            10.8             1.6               0
   21            10.8             1.8               0
   22            10.8               2               0
   23             9.6               0               0
   24             9.6             0.2               0
   25             9.6             0.4               0
   26             9.6             0.6               0
   27             9.6             0.8               0
   28             9.6               1               0
   29             9.6             1.2               0
   30             9.6             1.4               0
   31             9.6             1.6               0
   32             9.6             1.8               0
   33             9.6               2               0
   34             8.4               0               0
   35             8.4             0.2               0
   36             8.4             0.4               0
   37             8.4             0.6               0
   38             8.4             0.8               0
   39             8.4               1               0
   40             8.4             1.2               0
   41             8.4             1.4               0
   42             8.4             1.6               0
   43             8.4             1.8               0
   44             8.4               2               0
   45             7.2               0               0
   46             7.2             0.2               0
   47             7.2             0.4               0
   48             7.2             0.6               0
   49             7.2             0.8               0
   50             7.2               1               0
   51             7.2             1.2               0
   52             7.2             1.4               0
   53             7.2             1.6               0
   54             7.2             1.8               0
   55             7.2               2               0
   56               6               0               0
   57               6             0.2               0
   58               6             0.4               0
   59               6             0.6               0
   60               6             0.8               0
   61               6               1               0
   62               6             1.2               0
   63               6             1.4               0
   64               6             1.6               0
   65               6             1.8               0
   66               6               2               0
   67             4.8               0               0
   68             4.8             0.2               0
   69             4.8             0.4               0
   70             4.8             0.6               0
   71             4.8             0.8               0
   72             4.8               1               0
   73             4.8             1.2               0
   74             4.8             1.4               0
   75             4.8             1.6               0
   76             4.8             1.8               0
   77             4.8               2               0
   78             3.6               0               0
   79             3.6             0.2               0
   80             3.6             0.4               0
   81             3.6             0.6               0
   82             3.6             0.8               0
   83             3.6               1               0
   84             3.6             1.2               0
   85             3.6             1.4               0
   86             3.6             1.6               0
   87             3.6             1.8               0
   88             3.6               2               0
   89             2.4               0               0
   90             2.4             0.2               0
   91             2.4             0.4               0
   92             2.4             0.6               0
   93             2.4             0.8               0
   94             2.4               1               0
   95             2.4             1.2               0
   96             2.4             1.4               0
   97             2.4             1.6               0
   98             2.4             1.8               0
   99             2.4               2               0
  100             1.2               0               0
  101             1.2             0.2               0
  102             1.2             0.4               0
  103             1.2             0.6               0
  104             1.2             0.8               0
  105             1.2               1               0
  106             1.2             1.2               0
  107             1.2             1.4               0
  108             1.2             1.6               0
  109             1.2             1.8               0
  110             1.2               2               0
  111               0               0               0
  112               0             0.2               0
  113               0             0.4               0
  114               0             0.6               0
  115               0             0.8               0
  116               0               1               0
  117               0             1.2               0
  118               0             1.4               0
  119               0             1.6               0
  120               0             1.8               0
  121               0               2               0];
connec=[
1 100 101 112 111
2 89 90 101 100
3 78 79 90 89
4 67 68 79 78
5 56 57 68 67
6 45 46 57 56
7 34 35 46 45
8 23 24 35 34
9 7 9 24 23
10 1 2 9 7
11 101 102 113 112
12 90 91 102 101
13 79 80 91 90
14 68 69 80 79
15 57 58 69 68
16 46 47 58 57
17 35 36 47 46
18 24 25 36 35
19 9 10 25 24
20 2 3 10 9
21 102 103 114 113
22 91 92 103 102
23 80 81 92 91
24 69 70 81 80
25 58 59 70 69
26 47 48 59 58
27 36 37 48 47
28 25 26 37 36
29 10 11 26 25
30 3 4 11 10
31 103 104 115 114
32 92 93 104 103
33 81 82 93 92
34 70 71 82 81
35 59 60 71 70
36 48 49 60 59
37 37 38 49 48
38 26 27 38 37
39 11 13 27 26
40 4 5 13 11
41 104 105 116 115
42 93 94 105 104
43 82 83 94 93
44 71 72 83 82
45 60 61 72 71
46 49 50 61 60
47 38 39 50 49
48 27 28 39 38
49 13 14 28 27
50 5 6 14 13
51 105 106 117 116
52 94 95 106 105
53 83 84 95 94
54 72 73 84 83
55 61 62 73 72
56 50 51 62 61
57 39 40 51 50
58 28 29 40 39
59 14 16 29 28
60 6 8 16 14
61 106 107 118 117
62 95 96 107 106
63 84 85 96 95
64 73 74 85 84
65 62 63 74 73
66 51 52 63 62
67 40 41 52 51
68 29 30 41 40
69 16 18 30 29
70 8 12 18 16
71 107 108 119 118
72 96 97 108 107
73 85 86 97 96
74 74 75 86 85
75 63 64 75 74
76 52 53 64 63
77 41 42 53 52
78 30 31 42 41
79 18 20 31 30
80 12 15 20 18
81 108 109 120 119
82 97 98 109 108
83 86 87 98 97
84 75 76 87 86
85 64 65 76 75
86 53 54 65 64
87 42 43 54 53
88 31 32 43 42
89 20 21 32 31
90 15 17 21 20
91 109 110 121 120
92 98 99 110 109
93 87 88 99 98
94 76 77 88 87
95 65 66 77 76
96 54 55 66 65
97 43 44 55 54
98 32 33 44 43
99 21 22 33 32
100 17 19 22 21];
pointload_complete=[6 2 -1];
dirichlet_data=[121 1 0
           121 2 0
           120 1 0
            120 2 0 
            119 1 0
            119 2 0
            118 1 0
            118 2 0
            117 1 0
            117 2 0
            116 1 0
            116 2 0
            115 1 0
            115 2 0
            114 1 0
            114 2 0
            113 1 0
            113 2 0
            112 1 0
            112 2 0
            111 1 0
            111 2 0];