ORIGIN 1

1. Coordenas de los nodos

x := ( 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 0 0 0 0 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.
y := ( 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 0 0 0 0 2.5 2.5 2.5 2.5 2.
z := ( 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 0 0 0 0 1.
x := ( 7 7 7 7 4 4 4 4 5.5 5.5 5.5 5.5 7 4 4 7 7 4 4 7 7 7 4 4 5.5 5.5 7 4 4 7
y := ( 10 5 5 10 10 5 5 10 10 5 5 10 7.5 7.5 7.5 7.5 10 10 5 5
z := ( 11 11 8 8 11 11 8 8 11 11 8 8 11 11 8 8 9.5 9.5 9.5 9.5 11 8 11 8 11 8 11
T

x := x

y := y

z := z

2. Conectividad de los Elementos

a := 5

b := 5

la := 2 ( a + b )

la = 20

c := 4

d := 3

lb := 2 ( c + d )

lb = 14

h := 3

 1
2
3
4

5
6

7
8
9

10
c :=
11
12

13
14

15
16
17

18
19
20

6 23 35

22

34

21 33 45 6 23 35 47

22

34

46

2 21 33 5

22 34 46 7 24 36 48 90

102 109 116 92

7 24 36 91

90

102 109 93

6 23 35 57 58 69 76

24

36

60

23 35 47 58 69 76 83

24

36

48

59

3 22 34 8

24 36 48 59 70 77 84 92

103 110 117 123

7 24 36 60 59 70 77 93

92

103 110 124

10 25 37 61 62 71 78 12

11

26

38

64

25 37 49 62 71 78 85 11

26

38

50

63

26 38 50 63 72 79 86 94

104 111 118 125

11 26 38 64 63 72 79 95

94

104 111 126

27 39 51 14 28 40 52 16

30

42

54

15

28 40 52 65 73 80 87 15

29

41

53

66

29 41 53 66 74 81 88 96

105 112 119 127

30 42 54 15 29 41 53 97

106 113 120 96

20 32 44 18 19 31 43 98

101 108 115 99

19 31 43 67 68 75 82 99

100 107 114 128

31 43 55 68 75 82 89 100 107 114 121 129

32 44 56 19 31 43 55 101 108 115 122 100

3. Definir Grados de libertad

i := 1 .. 20
e := 1 .. 14
G

e , 3 i

:= 3 c

i, e

e , 3 i 1

:= 3 c

3. Funciones de Interpolacion
Coordenadas de los nodos:

i, e

e , 3 i 2

:= 3 c

i, e

48
117

110

77
84

130

131
79

86
132

133

53
88

134

119
114

135
136

121
36

 1
1
1
1

1
1

1
1
1

1
:=
1
1

0
0

0
0
1

1
1
1

1
1
1
1

1
1

1
1
0

0
:=
0
0

1
1

1
1
1

1
1
1

1
1
1
1

1
1

1
1
1

1
:=
1
1

1
1

1
1
0

0
0
0

Nodos en los vertices:

i := 1 .. 8
N

i, j

:=

1
8

j := 1 .. 20

)(

)(

)(

1 + 1 + 1 + + + 2
i j

i j

i j

i j

i j

i j

Nodos con =0:

i := 9 .. 12
N

i, j

:=

1
4

( j) (1 + i j) (1 + i j)
2

Nodos con =0
i := 13 .. 16
N

i, j

:=

1
4

( j) (1 + i j) (1 + i j)
2

Nodos con =0

2
1

3
0

4
0

7
N= 8

10

11

12

13

14

15

16

i := 17 .. 20

i, j

:=

( j) (1 + i j) (1 + i j)
2

5. Coordenadas puntos de integracion de Gauss:

0.57735
0.57735

0.57735
0.57735
g :=
g :=
0.57735
0.57735
0.57735

0.57735

0.57735
0.57735

0.57735
0.57735

g :=
0.57735
0.57735
0.57735

0.57735

0.57735
0.57735

0.57735
0.57735

0.57735
0.57735
0.57735

0.57735

5. Matriz Elastica D:
:= 0.2
:=

E := 3000000
E

( 1 + ) ( 1 2)

:=

E
2 ( 1 + )


+ 2

+ 2

+ 2
D :=
0
0
0
0
0
0

0
0
0

0 0 0
0 0 0

0 0 0
0 0

0 0
0 0

5. Matriz B
Nodos en los vertices:
i := 1 .. 8

j := 1 .. 8

:=

:=

:=

i, j

i, j

i, j

)(

)(

)(

1 + g 1 + g 2 g + g + g 1
i

)(

)(

:=

)(

:=

1 g

:=

1 g

i, j

( j) i (1 + i gj)

)(

g 1 + g 1 + g
j

i, j

( j) (1 + i g j)i

Nodos con =0
i := 13 .. 16
N

i, j

i, j

1 + g 1 + g g + g + 2 g 1
i j
i
j
i j
i
j
i j
8 i

i := 9 .. 12
i, j

1 + g 1 + g g + 2 g + g 1

Nodos con =0:

:=

:=

1
4

1 g

1
2

( j) i (1 + i gj)

)(

g 1 + g 1 + g
j

j i

j i

)
)
)

i, j

:=

1 g

( j) (1 + i g j)i

Nodos con =0
i := 17 .. 20
:=

:=

:=

i, j

i, j

i, j

1 g

( j) i (1 + i g j)
2

1 g

( j) (1 + i g j)i

1
2

)(

g 1 + g 1 + g
j

6. Jacobiano:
i := 1 .. 20

J11

j,e

j := 1 .. 8

e := 1 .. 14

(Ni, jxc )

:=

J12

j,e

i,e

:=

(Ni, jyc )

J21

j,e

J22

j,e

i,e

:=

j,e

j,e

:= J11

:=

J32

j,e

i,e

J33

j,e

(Ni, jyc )

:=

J33

j,e

i,e

j,e

+ J12

J23

j,e

J31

j,e

j,e

+ J13

J21

j,e

J32

j,e

j,e

J13

(Ni, jzc )
i,e

:= J12

:= J13

:= J21

:= J22

:= J23

3 e 1 , 3 j 2

j,e

3 e 2 , 3 j 1

j,e

3 e 1 , 3 j 1

:= J31

j,e

3 e , 3 j 1

2
15

j,e
j,e

:= J32

j,e

3
0 447

1 332 10 15

3 e 2 , 3 j
3 e 1 , 3 j

3 e , 3 j

J22

j,e

3 e , 3 j 2

i,e

:= J11

(Ni, jzc )
i

3 e 2 , 3 j 2

j,e

J22

j,e

:=

J23

i,e

detJ

i,e

(Ni, jxc )

:=

(Ni, jzc )
i

(Ni, jyc )

J31

j,e

(Ni, jxc )

:=

:=

J13

i,e

J31

j,e

j,e

j,e
j,e

:= J33

j,e

5
15

0 258

J12

-1.33210-15

1.5

-0.447

1.5

0.258

6.091

8.72

3.591

1.5

1.77610-15

-0.258

1.5

-1.33210-15

1.5

2.5

2.5

1.5

1.77610-15

7
J= 8

1.5

-1.33210-15

1.5

2.5

2.5

1.5

1.77610-15

1.5

10

1.5

-1.33210-15

11

2.5

2.5
-3.55310-15

12

1.5

1.77610-15

13

-1.33210-15

14

2.5

2.5

15

1.5

-1.33210-15

...

16

6. Matriz B

e, j

:= submatrix( J , 3 e 2 , 3 e , 3 j 2 , 3 j )

Nx

20 ( e 1) + i , j

Ny

20 ( e 1) + i , j

Nz

20 ( e 1) + i , j

:=

J 1 N + J 1 N + J 1 N
( e , j ) 1 , 1 i , j ( e , j ) 1 , 2 i , j ( e , j ) 1 , 3 i , j

:=

J 1 N + J 1 N + J 1 N
( e , j ) 2 , 1 i , j ( e , j ) 2 , 2 i , j ( e , j ) 2 , 3 i , j

:=

J 1 N + J 1 N + J 1 N
( e , j ) 3 , 1 i , j ( e , j ) 3 , 2 i , j ( e , j ) 3 , 3 i , j

j := 1 .. 8
:= 0

j,i

48 ( e 1) + 6 j 5 , 3 i 2

48 ( e 1) + 3+ ( 6 j 5) , 3 i 2

:= Nx

20 ( e 1) + i , j

:= Ny

20 ( e 1) + i , j

48 ( e 1) + 1+ ( 6 j 5) , 3 i 1

48 ( e 1) + 3+ ( 6 j 5) , 3 i 1

48 ( e 1) + 4+ ( 6 j 5) , 3 i 1

:= Ny

:= Nx
:= Nz

20 ( e 1) + i

20 ( e 1) + i

20 ( e 1) + i ,

48 ( e 1) + 5+ ( 6 j 5) , 3 i 2

:= Nz

0.024

0.291

0.051

0.067

-0.114

0.051

0.024

0.067

0.291

-0.114

0.051

0.111

-0.114

0.024

0.111

8.3410-3

1.47810-3

-0.018

0.021

-0.116

0
8.3410-3

10

0.021

-0.018

1.47810-3

11

-0.116

0.021

-0.207

12

-0.116

-0.018

-0.207

13

0.303

0.032

14

-0.182

-0.139

15

1.025

16

-0.182

0.303

-0.139

...

K :=
e

7
B=

20 ( e 1) + i , j

submatrix[B , 48(e 1) + 6j 5 , 48(e 1) + 6j , 1 , 60] Dsubmatrix[B , 48(e

T

7. Ensamblaje de Matriz de Rigidez

m := 136
Kglobal

3m , 3m

i := 1 .. 60

Kglobal

1
:= 0

k := 1 .. 60

Ge , i , Ge , k

:= Kglobal

e := 1 .. 14

Ge , i , Ge , k

( e)i, k

+ K

Kglobal =

9.96

-6.23

3.40

7.23

2.26

1.08

5.45

-3.0

-5.68

10

2.76

11

-7.

12

9.10

13

4.44

14

8.04

15

6.41

16

5.6

8. Condicion de Borde
3
4

7
8
emp :=
11
12
15

16

e := 9 .. 14

Ge , 3 emp , Ge , 3 emp
i
i

Kglobal

Ge , 3 emp 1 , G e , 3 emp 1
i

Ge , 3 emp 2 , G e , 3 emp 2

:= 10
:= 10

Ge , 3 sim 2 , Ge , 3 sim 2

:= 10

10

:= 10

10

e := 13 .. 14
Ge , 3 sim 2 , Ge , 3 sim 2
i

9. Vector de Cargas
3m

:= 0

2( La+ bp+ 1)

10

i := 1 .. 8

Kglobal

10

e := 5 .. 8

Kglobal

10

5
6

7
8
sim :=
14
15
18

19

:= 10

Kglobal

Kglobal

i := 1 .. 8

:= 100

10. Vector Desplazamientos Nodales

desp := Kglobal

61

97

13

16

G= 7

67

103

10

10

11

64

12

100

13

22

14

106

.5 1.5 1.5 1.5 0 0 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 0 0 1.5 1.5 1.5 1.5 1.
.5 2.5 2.5 2.5 2.5 2.5 0 0 0 0 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 0 0 0 0 2.5 2.5 2.5 2.5 2.5
5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 0 0 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5

4 7 7 7 4 4 5.5 5.5 7 4 4 7 4 7 7 7 4 4 5.5 5.5 7 4 4 7 4 7 0 0 0 0 2 2 2 2

0 0 0 0 0 0 2.5 2.5 2.5 2.5 0 0 5 5 5 5 5 5 2.5 2.5 2.5 2.5 5 5 10 10 10 10 1

1 11 8 8 9.5 9.5 11 8 11 8 11 8 11 11 8 8 9.5 9.5 11 8 11 8 11 8 11 11 8 8 9.5 9.5

6
0

7
0

8
0

9
0

10
0

...

 J21

,e

J33

j,e

j,e

J11

J23

j,e

J32

j,e

j,e

,j

48 ( e 1) + 2+ ( 6 j 5) , 3 i

:= Nz

20 ( e 1) + i , j

,j

,j

48 ( e 1) + 4+ ( 6 j 5) , 3 i

:= Ny

20 ( e 1) + i , j

48 ( e 1) + 5+ ( 6 j 5) , 3 i

:= Nx

1) + 6 j 5 , 48 ( e 1 ) + 6 j , 1 , 60] detJ

20 ( e 1) + i , j

j , e

64106

-6.238105

3.404106

7.235106

2.261104

38105

8.616106

-2.289106

-4.683104

6.009106

04106

-2.289106

1.938107

6.746105

-1.319106

35106

-4.683104

6.746105

1.473107

3.198105

61104

6.009106

-1.319106

3.198105

9.897106

88105

-1.388106

1.31107

2.421106

-1.411106

55106

-6.543104

5.492104

7.536106

3.972105

07104

4.868106

-7.111105

3.972105

5.942106

88105

-7.111105

1.155107

2.408105

-7.905105

65106

-7.9104

2.067105

1.725106

3.487105

.9104

2.289106

-2.191105

3.14105

1.246106

01104

-1.496105

5.142106

4.159105

-4.255105

49106

1.499105

-4.022104

7.104106

3.39105

44104

2.173106

-3.785105

3.39105

3.736106

14105

-3.785105

2.84106

-1.157106

-4.158105

63106

-1.14104

4.246105

1.178107

...

10

10

61

62

63

64

65

66

62

63

97

98

99

100

101

102

64

98

99

133

134

135

136

137

138

100

14

15

16

17

18

19

20

21

22

17

18

67

68

69

70

71

72

19

68

69

103

104

105

106

107

108

70

104

105

139

140

141

142

143

144

106

11

12

268

269

270

271

64

65

66

304

305

306

268

65

66

100

101

102

325

326

327

304

101

102

136

137

138

346

347

348

325

23

24

19

20

21

274

275

276

277

107

108

142

143

144

349

350

351

...

5 1.5 1.5 1.5 1.5 1.5 0 0 1.5 1.5 1.5 1.5 1.5 1.5 2 2 2 2 0 0 0 0 2 2 2
2.5 2.5 2.5 0 0 0 0 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 0 0 2.5 2.5 2.5 2.5 2.5
0 0 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 0 0 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1

0 0 0 0 0 0 2 2 0 0 0 0 0 2 2 0 0 0 0 0 2 2 0 0 0 7 7 4 4 5.5 5.5 4 7 7 4 4
0 10 7.5 7.5 7.5 7.5 10 10 10 5 5 10 10 5 5 10 7.5 7.5 10 5 0 0 0 0 2.

11 11 8 8 11 11 8 8 11 8 9.5 9.5 11 8 11 8 11 8 9.5 11 8 11 8 11 8 9.5 11 8 11

:=

2 2 2 0 0 2 2 2 2 2 0 0 2 2 2 2 2 0 0 2 2 2 1.5 1.5 1.5 1.5 0

2.5 0 0 2.5 2.5 2.5 2.5 2.5 0 0 2.5 2.5 2.5 2.5 2.5 0 0 2.5 2.5 2.5 2.5 2.5 2.5 2.5 0
1.5 1.5 0 0 1.5 1.5 1.5 1.5 1.5 1.5 0 1.5 1.5 1.5 1.5 1.5 1.5 0 1.5 1.5 1.5 1.5 1.5

7 7 4 5.5 4 7 4 7 7 4 5.5 4 7 4 7 7 4 5.5 4 7 4 7 4 4 5.5 5.5 4 4 4 4 4 5.5 5.5

.5 2.5 0 5 5 5 5 2.5 2.5 5 10 10 10 10 7.5 7.5 10 5 10 5 10 5 10 7.5 7.5 1

8 11 8 9.5 0 0 0 0 0 0 0 0 4 4 4 4 0 0 0 0 0 4 4 0 0 0 0 0 4 4 0 0 0 0 0 4

0 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 0 1.5 1.5 1.5 1.5 1.5 1.5 0 1.5 1.5 1.5 1.5 1.5
0 2.5 2.5 2.5 2.5 2.5 2.5 2.5 0 0 2.5 2.5 2.5 2.5 2.5 0 0 2.5 2.5 2.5 2.5 2.5 0 0 2.5 2.
5 1.5 0 4 4 4 4 4 4 4 4 0 0 0 0 4 4 4 4 4 0 0 4 4 4 4 4 0

5 4 4 4)
10 10 5 5 0 0 0 2.5 2.5 0 0 5 5 5 2.5 2.5 5 5 10 10 10 7.5 7.5 10 10 5 10 5

4 0 0 0 0 0 4 4 0 0 0 0 0 4 4)

1.5 0 1.5 1.5 1.5 1.5 2 2 0 0 2 2 2 2 2 0 0 2 2 2 )

5 2.5 2.5 2.5 2.5 0 2.5 2.5 2.5 2.5 2.5 2.5 0 2.5 2.5 )
0 4 4 4 4 4 0 0 4 4 4 4 4 0 0 4 4 4 4 4 0 0 )

10 7.5 10 5 10 5 10 5 7.5 5 10 )