Presentasi Elemen Segitiga Ok

Unduh sebagai ppt, pdf, atau txt
Unduh sebagai ppt, pdf, atau txt
Anda di halaman 1dari 20

ELEMEN SEGITIGA 3 NODE

2
u
3

3
u
2

v
2

v
3

u
1

v
1

1
x
y
- Elemen 2 Dimensi
- 3 Node
- Tiap Node memiliki 2
displacement dalam bentuk
Translasi (u, v)
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
1. MEMILIH FUNGSI PENDEKATAN
FUNGSI INTERPOLASI LINEAR : u = a
1
+ a
2
x + a
3
y . (1)
v = a
1
+ a
2
x + a
3
y . (2)
UTK DISPLACEMENT ARAH u
u
1
= a
1
+ a
2
x
1
+ a
3
y
1

u
2
= a
1
+ a
2
x
2
+ a
3
y
2
u
3
= a
1
+ a
2
x
3
+ a
3
y
3

(
(
(

3
2
1
3 3
2 2
1 1
3
2
1
1
1
1
a
a
a
y x
y x
y x
u
u
u
DALAM BENTUK MATRIK
SETELAH DILAKUKAN INVERS

(
(
(

3
2
1
3 2 1
3 2 1
3 2 1
3
2
1
det
1
u
u
u
a
a
a

| | |
o o o
o
1
= x
2
y
3
x
3
y
2

|
1
= y
2
y
3

1
= x
3
x
2

o
2
= x
3
y
1
x
1
y
3
; o
3
= x
1
y
2
x
2
y
1

|
2
= y
3
y
1
; |
3
= y
1
y
2

2
= x
1
x
3
;
3
= x
2
x
1

MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
Konstanta a1, a2 dan a3 dimasukkan
ke persamaan (1), SEHINGGA DIDAPATKAN :
SHAPE FUNCTION : N
1
= 1/det ( o
1
+ x |
1
+ y
1
)
N
2
= 1/det (o
2
+ x |
2
+ y
2
)
N
3
= 1/det (o
3
+ x |
3
+ y
3
)
=
)
`

v
u
(

3 2 1
3 2 1
0 0 0
0 0 0
N N N
N N N

3
3
2
2
1
1
v
u
v
u
v
u
u = N
1
u
1
+ N
2
u
2
+ N
3
u
3

Untuk Perpindahan arah v juga sama, sehingga :
v = N
1
v
1
+ N
2
v
2
+ N
3
v
3

MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
MATRIK B
| | { } q N
x y
y
x
v
u
x y
y
x
xy
y
x
(
(
(
(
(
(

c
c
c
c
c
c
c
c
=
)
`

(
(
(
(
(
(

c
c
c
c
c
c
c
c
=

0
0
0
0

c
c
(
(
(




=
2 1 1 2 1 3 3 1 3 2 2 3
1 2 3 1 2 3
2 1 1 3 3 2
0 0 0
0 0 0
det
1
y y x x y y x x y y x x
x x x x x x
y y y y y y

3
3
2
2
1
1
v
u
v
u
v
u
det = x
2
y
3
x
3
y
2
+ x
1
(y
2
y
3
) + y
1
(x
3
x
2
)
{ } | | { } q B = c
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
[k] = matrik kekakuan elemen
h = tebal struktur
[B] = matrik hub. regangan & perpindahan.
[C] = matrik hub. tegangan dan regangan.
2. MENURUNKAN PERS. MATRIK K
}
= dA [B] [C] [B]
T
h k
Utk Kasus Plane Stress :
Dari Prinsip Energi Regangan
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
3. MENGHITUNG MATRIK BEBAN TOTAL
R = Q
NF
+ Q
BF
+ Q
T


Q
NF
= beban pada konsentrasi nodal
Q
BF
= beban body force (akibat beban sendiri)
Q
T
= beban traksi (surface traction)
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
4. ASSEMBLY ELEMEN
[K]. {q} = {R}

[K] = matrik kekakuan elemen GLOBAL.
{q} = vektor perpindahan simpul.
{R} = matrik beban total.
k
1.1
q
1
+ k
1.2
q
2
+ ....................... + k
1.n
q
n
= R
1

k
2.1
q
1
+ k
2.2
q
2
+ ....................... + k
2.n
q
n
= R
2

.
k
n.1
q
1
+ k
n.2
q
2
+ ....................... + k
n.n
q
n
= R
n

MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
5. MENDAPATKAN OUTPUT UTAMA
BERUPA PERPINDAHAN
{q} = [K]
-1
{R}
Syarat batas dimasukkan pada perpindahan
simpul (q) dimana syarat batas memberikan
informasi bagaimana struktur ditopang
dalam ruang, dengan memasukkan nilai
perpindahan yang telah ditetapkan sesuai
kondisi pada struktur.
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
6. OUTPUT SEKUNDER BERUPA
REGANGAN DAN TEGANGAN

e
e
(
(
(
(

t xy
y
x
2
xy
y
x
2
v 1
0 0
0 1 v
0 v 1
) v 1 (
E

PLANE STRESS PLANE STRAIN

e
e
(
(
(
(

+
=

t xy
y
x
xy
y
x
2
v 2 1
0 0
0 v 1 v
0 v v 1
) v 2 1 )( v 1 (
E

{ o }= [C] {c}
{c} = [ B] {q}
{o} = [C] [B] {q}
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
2
(6, 5)
3
(2, 1)
(6, 2)
1
x
y
T = W (lb/in
2
)
CONTOH KASUS (1)
Diketahui : v = 0.3 dan E = 30 x 10
6
psi
tebal h
x y
2 1
6 2
6 5
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
[K]. {q} = {R}
[K]. {q} = Q
NF
+ Q
T
+ Q
BF

{ } | |
} }
)
`

=
S
y
x T
y
x
y
x
A
dS
T
T
N h
F
F
F
F
q dA h
0
0
[B] [C] [B]
2
2
1
1
T
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
| |
(
(
(
(

=
2
1
0 0
0 1
0 1
) 1 (
2
v
v
v
v
E
C
Diketahui : v = 0.3
E = 30 x 10
6
psi

| |
(
(
(
(

=
2
3 . 0 1
0 0
0 1 3 . 0
0 3 . 0 1
) 3 . 0 1 (
) 10 ( 30
2
6
C
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
det = x
2
y
3
x
3
y
2
+ x
1
(y
2
y
3
) + y
1
(x
3
x
2
) = 12
| |
(
(
(




=
2 1 1 2 1 3 3 1 3 2 2 3
1 2 3 1 2 3
2 1 1 3 3 2
0 0 0
0 0 0
det
1
y y x x y y x x y y x x
x x x x x x
y y y y y y
B
Dengan memasukkan koordinat x dan y
masing-masing node maka :
| |
(
(
(
(
(
(


=
12
1
3
1
4
1
3
1
4
1
0
3
1
0
3
1
0 0 0
0
12
1
0
4
1
0
4
1
B
| |
(
(
(




=
2 1 2 6 2 5 6 2 5 2 6 6
2 6 0 6 2 0 6 6 0
0 2 1 0 2 5 0 5 2
12
1
B
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
A h dA h k
A
[B] [C] [B] [B] [C] [B] ] [
T T
}
= =
Area = 0.5 det = 6
(
(
(
(
(
(
(




=
4 . 16
6 . 2 6 . 6
4 . 17 8 . 6 6 . 21
2 . 6 6 . 9 4 . 10 6 . 21
05 . 1 2 . 4 2 . 4 2 . 4 15 . 3
6 . 3 3 6 . 3 12 0 9
) 10 ( 374 . 1 ] [
6
simetri
h k
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
| |
} }
)
`

(
(
(
(
(
(
(
(

=
)
`

=
S S
y
x T
T
dy
W
N
N
N
N
N
N
h dS
T
T
N h Q
0
0
0
0
0
0
0
3
3
2
2
1
1
Q
T
= Beban Traksi
(Surface Traction)

=
}
}
}
0
0
0
3
2
1
dy N
dy N
dy N
Wh
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
=

=
}
}
}
0
0
0
3
2
1
dy N
dy N
dy N
Wh Q
T

0
1
0
1
0
0
2
3
Wh
2
3
)) 4 )( ( ) 4 )( 6 ( 4 (
12
1
) (
det
1
5
2
2 2
5
2
2 2
= + + = + + =
} } }
dy y dy y x dy N
S
| o
2
3
)) 4 )( ( ) 1 )( 6 ( 2 (
12
1
) (
det
1
5
2
3 3
5
2
3 3
= + + = + + =
} } }
dy y dy y x dy N
S
| o
0 )) 0 )( ( ) 3 )( 6 ( 18 (
12
1
) (
det
1
5
2
1 1
5
2
1 1
= + + = + + =
} } }
dy y dy y x dy N
S
| o
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.

0
1
0
1
0
0
2
3
0
0
2
2
1
1
3
3
2
2
1
1
Wh
F
F
F
F
v
u
v
u
v
u
y
x
y
x
(
(
(
(
(
(
(




4 . 16
6 . 2 6 . 6
4 . 17 8 . 6 6 . 21
2 . 6 6 . 9 4 . 10 6 . 21
05 . 1 2 . 4 2 . 4 2 . 4 15 . 3
6 . 3 3 6 . 3 12 0 9
) 10 ( 374 . 1
6
simetri
h
[K]. {q} = {R}

Dg memasukkan kondisi batas yaitu di node 1 dan 2 yang dikenai
tumpuan sendi shg u
1
= v
1
= u
2
= v
2
= 0, maka :

(
(
(
(
(
(
(




4 . 16
6 . 2 6 . 6
4 . 17 8 . 6 6 . 21
2 . 6 6 . 9 4 . 10 6 . 21
05 . 1 2 . 4 2 . 4 2 . 4 15 . 3
6 . 3 3 6 . 3 12 0 9
) 10 ( 374 . 1
6
simetri
h

+
=

0
2
3
2
3
0
0
0
0
2
2
1
1
3
3
Wh
F
Wh F
F
F
v
u
y
x
y
x
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
PARTISIONING UKURAN MATRIK

u
3
= 1.76 (10
-7
)W psi dan v
3
= 2.8 (10
-8
)W psi

=
)
`


0
2
3
4 . 16 6 . 2
6 . 2 6 . 6
) 10 ( 374 . 1
3
3 6
Wh
v
u
h
UTK MENDAPATKAN HARGA u
3
DAN v
3
DILAKUKAN TEKNIK
PARTISIONING SHG UKURAN MATRIKNYA HANYA 2X2
(
(
(
(
(
(
(




4 . 16
6 . 2 6 . 6
4 . 17 8 . 6 6 . 21
2 . 6 6 . 9 4 . 10 6 . 21
05 . 1 2 . 4 2 . 4 2 . 4 15 . 3
6 . 3 3 6 . 3 12 0 9
) 10 ( 374 . 1
6
simetri
h

+
=

0
2
3
2
3
0
0
0
0
2
2
1
1
3
3
Wh
F
Wh F
F
F
v
u
y
x
y
x
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.
2
4
1
x
y
CONTOH KASUS (2) - PR
Diketahui : v = 0.5 dan E = 30 x 10
6
psi
tebal h
3
1000 lb
1000 lb
1
2
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.

=
1000
0
0
1000
2
2
1
1
y
x
y
x
NF
F
F
F
F
Q

=
4
4
3
3
0
0
0
0
v
u
v
u
q

=
8 8x
K
MEH by Dr. Eng., Moch. Agus Choiron, ST., MT.

Anda mungkin juga menyukai