Scribd
Upload
26 views8 pages

Contoh A. Frame 2D: STEP 1: Discretize and Global Degrees of Freedom (DOF)

1) The document describes the process of modeling a 2D frame structure for structural analysis using finite element method. 2) It involves discretizing the frame into elements, defining the element stiffness matrices in local and global coordinate systems, assembling the global stiffness matrix, and solving the system of equations to obtain displacements and internal forces in the elements. 3) The steps include defining degrees of freedom, establishing element properties, transforming between local and global coordinate systems, applying loads, solving for displacements, and determining element forces.

Uploaded by

Rendi Teguh Prayitno
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
26 views8 pages

Contoh A. Frame 2D: STEP 1: Discretize and Global Degrees of Freedom (DOF)

1) The document describes the process of modeling a 2D frame structure for structural analysis using finite element method. 2) It involves discretizing the frame into elements, defining the element stiffness matrices in local and global coordinate systems, assembling the global stiffness matrix, and solving the system of equations to obtain displacements and internal forces in the elements. 3) The steps include defining degrees of freedom, establishing element properties, transforming between local and global coordinate systems, applying loads, solving for displacements, and determining element forces.

Uploaded by

Rendi Teguh Prayitno
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 8

Contoh a.

Frame 2D

2 t/m'

10

0
5

3m
2

6m

STEP 1 : Discretize and Global Degrees of Freedom (DOF)


D3y
Y

Defined DOF
after boundary conditions

1
D3X

#1

3
R3Z
#2

X
2

global axis

STEP 2 : Define Element Stiffness Matrix in LOCAL AXIS


Element #1
E=
A=
I=
L=
Angel =

[k]1=

2.00E+06
0.18
0.0054
6
0

t/m2
m2
m4
m
degree

60000
0
0
-60000
0
0

0
600
1800
0
-600
1800

Cos =
Sin =
Cos2 =
Sin2 =

0
1800
7200
0
-1800
3600

1
0
1
0

-60000
0
0
60000
0
0

local axis

0
-600
-1800
0
600
-1800

0
1800
3600
0
-1800
7200

Element #2
E=
A=
I=

2
2.00E+06 t/m
m2
0.09
4
0.000675 m

Cos =
Sin =
Cos2 =

Page 1

0.0000
1.0000
0.0000

Local axis

L=
Angel =

[k]2=

3
90
60000
0
0
-60000
0
0

Sin2 =

m
degree
0
600
900
0
-600
900

0
900
1800
0
-900
900

1.0000

-60000
0
0
60000
0
0

0
-600
-900
0
600
-900

0
900
900
0
-900
1800

STEP 3 : Define Element Stiffness Matrix in GLOBAL AXIS [K] = [T]T.[k].[T]


Element #1

[T]1=

[K]1=

1
0
0
0
0
0

0
1
0
0
0
0

0
0
1
0
0
0

0
0
0
1
0
0

0
0
0
0
1
0

0
0
0
0
0
1

D 1x

D 1y

R 1z

D 3x

D 3y

R 3z

60000

-60000

D 1x

600

1800

-600

1800

D 1y

1800

7200

-1800

3600

R 1z

-60000

60000

D 3x

-600

-1800

600

-1800

D 3y

1800

3600

-1800

7200

R 3z

0
-1
0
0
0
0

1
0
0
0
0
0

0
0
1
0
0
0

0
0
0
0
-1
0

0
0
0
1
0
0

0
0
0
0
0
1

D 2x

D 2y

R 2z

D 3x

D 3y

R 3z

600

-900

-600

-900

D 2x

60000

-60000

D 2y

-900

1800

900

900

R 2z

-600

900

600

900

D 3x

-60000

60000

D 3y

-900

900

900

1800

R 3z

Element #2

[T]2=

[K]2=

STEP 4 : Assemble GLOBAL STRUCTURE STIFFNESS MATRIX

[K] =

D 3x

D 3y

R 3z

60600

900

D 3x

60600

-1800

D 3y

Page 2

900

-1800

R 3z

9000

STEP 5.a : Compute Equivalent Span Load (each element) in LOCAL axis and GLOBAL axis
0
-6
{fo}1=
lokal

-6

{Fo}1=

global

lokal

D 1x

-6

D 1y

-6

R 1z

D 3x

-6

-6

D 3y

R 3z

D 2x

D 2y

R 2z

0
{fo}2=

{Fo}2=

global

D 3x

D 3y

R 3z

STEP 5.b : Global structure load vector in GLOBAL axis

0
-6
6

{Fo} =

0
-10
5

{F} =

{Fo} + {F} =

0
-16
11

STEP 6 : Solve the Global Equation to find GLOBAL DISPLACEMENT


{D} = [K]-1 * ( {Fo} + {F} )
D 3x
D 3y

1.65263E-05 -4.9E-08 -1.663E-06


=

R 3z

-4.93815E-08

1.7E-05

3.325E-06

-16

-1.66251E-06

3.3E-06

0.0001119

11

Formula : {d} = [T] . {D}


Element #1 :

global

D 1x

D 1y

R 1z

{d} #1 =

-0.00002

D 3x

lokal

-0.00002

-0.00023

D 3y

-0.00023

0.00118

R 3z

0.00118

D 2x

D 2y

0
0

Element #2 :

Page 3

-0.00023
0.00118

STEP 7 : Compute DISPLACEMENT VECTOR (each element) in LOCAL axis

{D} #1 =

-0.00002

{D} #2 =
global

R 2z

{d} #2 =

lokal

-0.00023

-0.00002

D 3x

-0.00023

D 3y

0.00002

0.00118

R 3z

0.00118

STEP 9 : Compute ELEMENT FORCES (in LOCAL axis)


Formula : {f} = [k].{d} - {fo}
Element #1 :
f1 =
f2 =
f3 =
f4 =
f5 =
f6 =

1.0499
8.2581
10.6537
-1.0499
3.7419
2.8950

t
t
t.m
t
t
t.m

10.65365

1.04985

8.25812

Element #2 :
f1 =
f2 =
f3 =
f4 =
f5 =
f6 =

2.89505

1.04985

3.74188

13.74188
13.7419
1.0499
1.0446
-13.7419
-1.0499
2.1050

t
t
t.m
t
t
t.m

2.10495

1.04985

1.04985

1.04460

13.74188

Page 4

Page 5

TT =

TT.k =

TT =

1
0
0
0
0
0

0
1
0
0
0
0

0
0
0
1
0
0

0
0
0
0
1
0

0
0
0
0
0
1

60000

-60000

600

1800

-600

1800

1800

7200

-1800

3600

-60000

60000

-600

-1800

600

-1800

1800

3600

-1800

7200

0
1
0
0
0
0

-1
0
0
0
0
0

0
0
1
0
0
0

0
TT.k =

0
0
1
0
0
0

-600

0
0
0
0
1
0

0
0
0
-1
0
0

0
0
0
0
0
1

-900

600

-900

60000

-60000

900

1800

-900

900

600

900

-600

900

-60000

60000

900

900

-900

1800

Page 6

m
m
radian

Page 7

[k].{d}

{fo}

1.0499
2.2581
4.6537
-1.0499
-2.2581
8.8950

0
-6
-6
0
-6
6

13.7419
1.0499
1.0446
-13.7419
-1.0499
2.1050

0
0
0
0
0
0

Page 8

You might also like