-
Notifications
You must be signed in to change notification settings - Fork 0
/
CantileverToy_Triangular.m
executable file
·151 lines (120 loc) · 2.65 KB
/
CantileverToy_Triangular.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
%==================================================================
% General Data File
% Title: TRIANGLE
% Units: SI
% Dimensions: 2D
% Type of problem: Plane_Stress
% Type of Phisics: ELASTIC
% Micro/Macro: MACRO
%
%==================================================================
%% Data
Data_prb = {
'TRIANGLE';
'SI';
'2D';
'Plane_Stress';
'ELASTIC';
'MACRO';
};
%% Coordinates
% Node X Y Z
coord = [
1 2 0 0
2 2 -0.2 0
3 2 -0.4 0
4 1.5 0 0
5 1.5 -0.2 0
6 1.5 -0.4 0
7 1 0 0
8 1 -0.2 0
9 1 -0.4 0
10 0.5 0 0
11 0.5 -0.2 0
12 0.5 -0.4 0
13 0 0 0
14 0 -0.2 0
15 0 -0.4 0
];
%% Conectivities
% Element Node(1) Node(2) Node(3) Material
connec = [
1 14 11 13 0
2 11 10 13 0
3 15 12 14 0
4 12 11 14 0
5 11 8 10 0
6 8 7 10 0
7 12 9 11 0
8 9 8 11 0
9 8 5 7 0
10 5 4 7 0
11 9 6 8 0
12 6 5 8 0
13 5 2 4 0
14 2 1 4 0
15 6 3 5 0
16 3 2 5 0
];
%% Variable Prescribed
% Node Dimension Value
dirichlet_data = [
13 1 0
13 2 0
14 1 0
14 2 0
15 1 0
15 2 0
];
%% Force Prescribed
% Node Dimension Value
pointload_complete = [
2 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 = [
];