Scribd
Upload
28 views11 pages

Analisa Struktur - Matriks

This document contains the steps to analyze a simple truss structure with 2 members and 3 nodes. It defines the node coordinates, degrees of freedom, member properties, and assembles the global stiffness matrix. It then solves for member end displacements and forces by applying a load vector and inverting the stiffness submatrix. Reactions are also calculated. Displacements are transformed to member local axes to determine axial strains and forces in each member.

Uploaded by

Irma Anggraini
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
28 views11 pages

Analisa Struktur - Matriks

This document contains the steps to analyze a simple truss structure with 2 members and 3 nodes. It defines the node coordinates, degrees of freedom, member properties, and assembles the global stiffness matrix. It then solves for member end displacements and forces by applying a load vector and inverting the stiffness submatrix. Reactions are also calculated. Displacements are transformed to member local axes to determine axial strains and forces in each member.

Uploaded by

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

SOAL 1

Global DOF Number for Each Member


Local Global DOF Number for Member
DOF
1 2
Nodal Coordinates Number
NODE X(m) Y(m) 1 1 1
1 0 0 2 2 2
2 -4 -2 3 3 5
3 0 -2 4 4 6

Member Data
Input Data Compute Data
Member
S Node E Node EA (Kn) ΔX ΔY L C S EA/L
1 1 2 1 -4 -2 4,4721 -0,8944 -0,4472 0,223609
2 3 2 1 -4 0 4 1 0 0,25

>> koord= [1 0 0; 2 -4 -2; 3 0 -2]

koord =

1 0 0

2 -4 -2

3 0 -2

>> Dof= [1 2; 1 1; 2 2; 3 5; 4 6]
Dof =

1 2

1 1

2 2

3 5

4 6

>> L1= sqrt((koord(2,3)-koord(1,3))^2+ (koord(2,2)-koord(1,2))^2)

L1 =

4.4721

>> L2= sqrt ((koord(3,3)-koord(2,3))^2 + (koord(3,2)-koord(2,2))^2)

L2 =

>> c1=(koord(2,2)- koord(1,2))/L1

s1=(koord(2,3)-koord(1,3))/L1

c1 =

-0.8944

s1 =

-0.4472

>> c2=(koord(3,2)-koord(2,2))/L2

s2=(koord(3,3)-koord(2,3))/L2
c2 =

s2 =

>> kc1= c1^2

kcs1= c1*s1

ks1= s1^2

kc1 =

0.8000

kcs1 =

0.4000

ks1 =

0.2000

>> kc2=c2^2

kcs2=c2*s2

ks2=s2^2

kc2 =

kcs2 =

0
ks2 =

>> EA= 1

EA =

>> K1=(EA/L1)*[kc1 kcs1 -kc1 -kcs1; kcs1 ks1 -kcs1 -ks1; -kc1 -kcs1 kc1 kcs1; -kcs1 -ks1 kcs1 ks1]

>> K2=(EA/L2)*[kc2 kcs2 -kc2 -kcs2; kcs2 ks2 -kcs2 -ks2; -kc2 -kcs2 kc2 kcs2; -kcs2 -ks2 kcs2 ks2]

>> K20= [0 0]

K20 =

0 0

>> K40=[0 0 0 0]

K40 =

0 0 0 0
>> K60=[0 0 0 0 0 0]

K60 =

0 0 0 0 0 0

>> KG_1=[K1 K40' K40'; K60; K60]

>> KG_2= [K2(1:2,1:2) K20' K20' K2(1:2,3:4); K60; K60; K2(3:4,1:2) K20' K20' K2(3:4,3:4)]

>> KGS= KG_1+KG_2


>> Pf= [0 ; -24]

Pf =

-24

>> KGS_F= KGS (1:2,1:2)

KGS_F =

0.4289 0.0894

0.0894 0.0447

>> d= inv(KGS_F)*Pf

d=

192.0000

-920.6563

>> Kff_new = KGS (1:2,1:2)

Kff_new =

0.4289 0.0894

0.0894 0.0447

>> Ksf_new = KGS (3:6,1:2)


>> df= inv(Kff_new)*Pf

df =

192.0000

-920.6563

>> Ps= Ksf_new * df

Ps =

48.0000

24.0000

-48.0000

>> KG_SF= KGS(3:6,:)

>> d_all=[d;0;0;0;0]

d_all =

192.0000

-920.6563

0
>> Reaksi = KG_SF*d_all

Reaksi =

Reaksi =

48.0000

24.0000

-48.0000

Cara pemendekan dan perpanjangan

>> dG= [d' 0 0 0 0]'

dG =

192.0000

-920.6563

>> T1 = [-c1 -s1 c1 s1]

T1 =

0.8944 0.4472 -0.8944 -0.4472

>> T2 = [-c2 -s2 c2 s2]

T2 =

-1 0 1 0
>> dl_1= T1*dG([1 2 3 4])

dl_1 =

-240.0000

>> F1= EA*L1*dl_1

F1 =

-1.0733e+03

>> dl_2= T2*dG([1 2 5 6])

dl_2 =

-192

>> F2= EA*L2*dl_2

F2 =

-768

>> Hasil_delta_aksial= [dl_1 F1; dl_2 F2]

Hasil_delta_aksial =

1.0e+03 *

-0.2400 -1.0733

-0.1920 -0.7680

>> PfPs= [Pf;Ps]


PfPs =

-24.0000

48.0000

24.0000

-48.0000

>> dG_x = dG([1 3 5])

dG_x =

192

>> dG_y= dG([2 4 6])

dG_y =

-920.6563

>> PfPs_x= PfPs ([1 3 5])

PfPs_x =

48.0000

-48.0000

>> PfPs_y= PfPs([2 4 6])


PfPs_y =

-24.0000

24.0000

>> Hasil_Node= [dG_x dG_y PfPs_x PfPs_y]

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy