Stiffness Method of Structural Analysis (Space Truss)

Givin Data Load Vector

E 1E+08 Gpa P1 200 cos²ɵx
A 0.0036 m2 P2 -150 cosɵxcosɵy

Coordinates 1 -10 0 0 P3 -250 K EA cosɵxcosɵz

of Joints 2 0 0 4 L
3 10 0 0
4 0 0 -4
5 0 12 0

Structural Stiffnes Matrics

Global stiffness matrics of member 1 member 2
L 15.62049935 EA 23046.63839 L
Cosɵx 0.64018 L Cosɵx
Cosɵy 0.76822 Cosɵy
Cosɵz 0 Cosɵz
4 5 6 1 2 3
9445.3436 11334.41232 0 -9445.34 -11334.4123 04
K1 11334.412 13601.29479 0 -11334.4 -13601.295 05 K2
0 0 0 0 0 06
-9445.344 -11334.4123 0 9445.344 11334.41232 01
-11334.41 -13601.295 0 11334.41 13601.29479 02
0 0 0 0 0 03

member 3
L 15.62049935 EA 23046.63839
Cosɵx -0.6401844 L
Cosɵy 0.76822128
Cosɵz 0
10 11 12 1 2 3
9445.3436 -11334.4123 0 -9445.34 11334.41232 0 10
K3 -11334.41 13601.29479 0 11334.41 -13601.2948 0 11
0 0 0 0 0 0 12 K4
-9445.344 11334.41232 0 9445.344 -11334.4123 01
11334.412 -13601.295 0 -11334.4 13601.29479 02
0 0 0 0 0 03

Structural Stiffness Matrix

1 2 3
18890.6872 0 01
S 0 78431.49 02
0 0 5692.1 3
Determination Member End
Displacement
 P=Sd 200 18890.6872 0 0
-150 = 0 78431.49 0 *
-250 0 0 5692.1

d=(S-1)P d1 5.29361E-05 0 0
d2 = 0 1.27E-05 0 *
d3 0 0 0.000176

Member Displecment In Global Coordinat System

04 07
05 08
V1 = 06 V2 = 09 V3 =
0.0105872 1 0.010587227 1
-0.001912 2 -0.0019125 2
-0.043921 3 -0.04392052 3

Member Desplecement in Local Coordinate System

member 1 Member 2
T1 = 0.6401844 0.76822128 0 0 0 0 T2
0 0 0 0.640184 0.76822128 0

u1 = T1*V1 = 0 u2
0.005309
member 3 Member 4
T3 = -0.640184 0.76822128 0 0 0 0 T4
0 0 0 -0.64018 0.76822128 0

u3 = T3*V3 = 0 u4
-0.00825

Stiffness matrix in Local Coordinate system

member 1
member 2 member 3
k1 = 23046.638 -23046.6384 k2= 28460.49894 -28460.5 k3=
-23046.64 23046.63839 -28460.4989 28460.5

member end force in Local Coordinate System

Q= k*u
member 1 member 2 member 3
Q1 = -122.3444 Q2= -343.647 Q3= 190.0656
122.34439 343.6473 -190.066

Member End Force In Global Coordinat System

F=T(t)*Q
member 1
F1 0.6401844 0 -78.323 4
F2 0.76822128 0 -122.344387 = -93.9876 5
F1=TtU1 F3 = 0 0 * 122.3443875 06
F4 0 0.640184 78.32297 1
F5 0 0.768221 93.98756 2
F6 0 0 03

member 3
uss)

cosɵxcosɵy cosɵxcosɵz
cos²ɵy cosɵycosɵz
cosɵycosɵz cos²ɵz
cos²ɵx
cos²ɵy
cos²ɵz

member 2
12.64911 EA 28460.5
0 L
0.948683
-0.31623
7 8 9 1 2 3
0 0 0 0 0 07
0 25614.449 -8538.15 0 -25614.4 8538.15 8
0 -8538.1497 2846.05 0 8538.15 -2846.05 9
0 0 0 0 0 01
0 -25614.449 8538.15 0 25614.45 -8538.15 2
0 8538.1497 -2846.05 0 -8538.15 2846.05 3

member 4
L 12.649111 EA 28460.5
Cosɵx 0 L
Cosɵy 0.9486833
Cosɵz 0.3162278
13 14 15 1 2 3
0 0 0 0 0 0 13
0 25614.449 8538.15 0 -25614.4 -8538.15 14
0 8538.1497 2846.05 0 -8538.15 -2846.05 15
0 0 0 0 0 01
0 -25614.449 -8538.15 0 25614.45 8538.15 2
0 -8538.1497 -2846.05 0 8538.15 2846.05 3

Cofactor Of S Cofactor Transpose

4.46E+08 0 0 4.46E+08 0 0
det S 8.434E+12 0 1.08E+08 0 0 1.08E+08 0
0 0 1.48E+09 0 0 1.48E+09

S inverse
5.29E-05 0 0
 d1 0 1.27E-05 0
d2 0 0 0.000176
d3

200 0.010587 1
-150 = -0.00191 2
-250 -0.04392 3

0 10 0 13
0 11 0 14
0 12 V4 = 0 15
0.010587 1 0.010587 1
-0.00191 2 -0.00191 2
-0.04392 3 -0.04392 3

Member 2
= 0 0.948683 -0.31623 0 0 0
0 0 0 0 0.948683 -0.31623

= T2*V2 = 0
0.012075
Member 4
= 0 0.948683 0.316228 0 0 0
0 0 0 0 0.948683 0.316228

= T4*V4 = 0
-0.0157

member 3 member 4
23046.64 -23046.638 k4= 28460.5 -28460.5
23046.64 23046.638 -28460.5 28460.5

member 4
Q4= 446.92213
-446.92213
member 2
F1 0 0 0
F2 0.948683 0 -343.647 = -326.012
F2=TtU2 = F3 = -0.31623 0 * 343.6473 108.6708
F4 0 0 0
F5 0 0.948683 326.0124
F6 0 -0.31623 -108.671

member 4
7
8
9
1
2
3
 Stiffness Method of Structural Analysis (Grid)
1 Develop Structural Stiffness Matrix Givin Data E
Coordinates of joints 1 0 0 0 I
2 7 0 0 EI=
3 0 0 8
4 7 0 8
Local Stiffness Matrix
member 1
L 10.6301458 m EI 57.77525
Cosɵ 0.65850461 L^3
Sin ɵ 0.75257669

12 0 6L
k1= EI 0 GJL^2/EI 0
L^3 6L 0 4L^2 =
-12 0 (-)6L
0 (-)GJL^2/EI 0
6L 0 2L^2

Transformation matrix
1 0 0 0 0 0 1 0
0 0.658505 0.75257669 0 0 0 0 0.6585046
T= 0 -0.75258 0.65850461 0 0 0 T^t= 0 0.7525767
0 0 0 1 0 0 0 0
0 0 0 0 0.658505 0.752577 0 0
0 0 0 0 -0.75258 0.658505 0 0
Stiffness Matrix In Global Coordinate System
5 6 7 12 3
693.303 -2773.21 2426.56 -693.303 -2773.212
K1= -2773.21 15146.99 -12534.2 2773.2118 7038.7066
2426.56 -12534.2 11789.61 -2426.56 -6878.285
-693.303 2773.212 -2426.56 693.30296 2773.2118
-2773.21 7038.707 -6878.28 2773.2118 15146.988
2426.56 -6878.28 5196.309 -2426.56 -12534.2
member 2
L 7 EI 202.3324
Cosɵ 1 L^3
Sin ɵ 0

k2=
 Transformation matrix
1 0 0 0 0 0 1 0
0 1 0 0 0 0 0 1
T= 0 0 1 0 0 0 T^t= 0 0
0 0 0 1 0 0 0 0
0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0

5 6 7 1 2
2427.988 0 8497.959 -2427.988 0
K2= 0 1248.571 0 0 -1248.571
8497.959 0 39657.14 -8497.959 0
-2427.99 0 -8497.96 2427.9883 0
0 -1248.57 0 0 1248.5714
8497.959 0 19828.57 -8497.959 0
member 3
L 8 EI 135.5469
Cosɵ 0 L^3
Sin ɵ 1

k3=

Transformation matrix
1 0 0 0 0 0 1 0
0 0 1 0 0 0 0 0
T= 0 -1 0 0 0 0 T^t= 0 1
0 0 0 1 0 0 0 0
0 0 0 0 0 1 0 0
0 0 0 0 -1 0 0 0

1 2 8 3 4
1626.563 -6506.25 0 -1626.563 -6506.25
K3= -6506.25 34700 0 6506.25 17350
0 0 1092.5 0 0
-1626.56 6506.25 0 1626.5625 6506.25
-6506.25 17350 0 6506.25 34700
0 0 -1092.5 0 0
member 4
L 7 EI 202.3324
Cosɵ -1 L^3
Sin ɵ 0
 k4=

Transformation matrix
1 0 0 0 0 0 1 0
0 -1 0 0 0 0 0 -1
T= 0 0 -1 0 0 0 T^t= 0 0
0 0 0 1 0 0 0 0
0 0 0 0 -1 0 0 0
0 0 0 0 0 -1 0 0

9 10 11 3 4
2427.988 0 -8497.96 -2427.988 0
K4= 0 0 0 0 0
-8497.96 0 39657.14 8497.9592 0
-2427.99 0 8497.959 2427.9883 0
0 0 0 0 0
-8497.96 0 19828.57 8497.9592 0
Structural Stiffness matrix
1 2 3 4
4054.551 -6506.25 -1626.56 -6506.25 1
S= -6506.25 35948.5714 6506.25 17350 2 det S 5.375E+16
-1626.56 6506.25 19201.54 -6027.95 3
-6506.25 17350 -6027.95 46489.61 4

End Displecement Member displecement in Global coordinate

P-Pf=Sd d=S*(P-Pf) V1 V2
0 5 0 5
-0.0560186 1 0 6 0 6
d= -0.0036041 2 0 7 0 7
-0.0151813 3 -0.01518 3 -0.056019 1
-0.0028537 4 -0.00285 4 -0.003604 2
0 12 0 8

Member displecement in Local coordinate

u1 u2
0 5 0 5
u=TV 0 6 0 6
0 7 0 7
-0.01518 3 -0.056019 1
-0.00188 4 -0.003604 2
0.002148 12 0 8
member end force (Local)
 Q1 Q2 Q3 Q4
7.914005 5 170.3875 5 75.5916 1 71.204909 9
1.545056 6 4.499945 6 0 2 0 10
Q=ku+Qf 28.04234 7 541.6689 7 -91.1244 8 194.63468 11
-7.91401 3 -120.388 1 124.4084 3 -21.23491 3
-1.54506 4 -4.49995 2 0 4 0 4
56.08468 12 476.0439 8 -104.143 12 129.00968 12

Reaction force(R)

R5 5
R6 6
R7 7
R= R8 8 =
R9 9
R10 10
R11 11
R12 12
 Member Load
200 GPA G 76 GPA 3
347000000 mm4 J 115000000 mm4 2
69400 4

5 6 7 12 3 4 Fixed End force In Lo

693.302955 0 3684.9557522 -693.302955 0 3684.956 5
0 822.19003897 0 0 -822.19 06
3684.95575 0 26114.411306 -3684.955752 0 13057.21 7 Qf1=
-693.30296 0 -3684.955752 693.302955 0 -3684.96 12
0 -822.190039 0 0 822.19 03
3684.95575 0 13057.205653 -3684.955752 0 26114.41 4

0 0 0 0
-0.7525767 0 0 0 693.303 0 3684.9557522
0.65850461 0 0 0 -2773.21 541.4159292 -19653.09735
0 1 0 0 T^tk1= 2426.56 618.7610619 17196.460177
0 0 0.6585046079 -0.752576695 -693.303 0 -3684.955752
0 0 0.7525766947 0.658504608 -2773.21 -541.4159292 -9826.548673
Global Coordinate System 2426.56 -618.7610619 8598.2300885
4 Fixed End force In Global Coordinat
2426.56034 5 0 5
-6878.2846 6 Ff1= 0 6
5196.30898 7 0 7
-2426.5603 12 0 3
-12534.198 3 0 4
11789.6134 4 0 12

5 6 7 1 2 8 Fixed End force In Lo

2427.98834 0 8497.9591837 -2427.988338 0 8497.959 5
0 1248.5714286 0 0 -1248.57 06
8497.95918 0 39657.142857 -8497.959184 0 19828.57 7 Qf2=
-2427.9883 0 -8497.959184 2427.988338 0 -8497.96 1
0 -1248.571429 0 0 1248.571 02
8497.95918 0 19828.571429 -8497.959184 0 39657.14 8
 Stiffness Matrix In Global Coordinate S
0 0 0 0
0 0 0 0 2427.988 0 8497.9591837
1 0 0 0 0 1248.571429 0
0 1 0 0 T^tk2= 8497.959 0 39657.142857
0 0 1 0 -2427.99 0 -8497.959184
0 0 0 1 0 -1248.571429 0
8497.959 0 19828.571429
8 Fixed End force In Global Coordinat
8497.95918 5 34.375 5
06 Ff2= 0 6
19828.5714 7 65.625 7
-8497.9592 1 15.625 1
02 0 2
39657.1429 8 0 8

1 2 8 3 4 12 Fixed End force In Lo

1626.5625 0 6506.25 -1626.5625 0 6506.25 1
0 1092.5 0 0 -1092.5 02
6506.25 0 34700 -6506.25 0 17350 8 Qf3=
-1626.5625 0 -6506.25 1626.5625 0 -6506.25 3
0 -1092.5 0 0 1092.5 04
6506.25 0 17350 -6506.25 0 34700 12

Stiffness Matrix In Global Coordinate S

0 0 0 0
-1 0 0 0 1626.563 0 6506.25
0 0 0 0 -6506.25 0 -34700
0 1 0 0 T^tk3= 0 1092.5 0
0 0 0 -1 -1626.56 0 -6506.25
0 0 1 0 -6506.25 0 -17350
0 -1092.5 0
12 Fixed End force In Global Coordinat
01 100 1
02 Ff3= 0 2
-1092.5 8 0 8
03 100 3
04 0 4
1092.5 12 0 12

9 10 11 3 4 12 Fixed End force In Lo

 2427.98834 0 8497.9591837 -2427.988338 0 8497.959 9
0 0 0 0 0 0 10
8497.95918 0 39657.142857 -8497.959184 0 19828.57 11 Qf4=
-2427.9883 0 -8497.959184 2427.988338 0 -8497.96 3
0 0 0 0 0 04
8497.95918 0 19828.571429 -8497.959184 0 39657.14 12

Stiffness Matrix In Global Coordinate S

0 0 0 0
0 0 0 0 2427.988 0 8497.9591837
-1 0 0 0 0 0 0
0 1 0 0 T^tk4= -8497.96 0 -39657.14286
0 0 -1 0 -2427.99 0 -8497.959184
0 0 0 -1 0 0 0
-8497.96 0 -19828.57143
12 Fixed End force In Global Coordinat Fixed End Forces
-8497.9592 9 34.345 9 115.625
0 10 Ff4= 0 10 Pf4= 0
19828.5714 11 -65.625 11 115.625
8497.95918 3 15.625 3 0
04 0 4
39657.1429 12 0 12

Adjugate matrix Inverse Of S

2.4314E+13 2.486724E+12 1.724643E+12 2.32901E+12 0.000452396 4.626885E-05
2.4867E+12 2.408649E+12 -8.11469E+11 -6.5611E+11 4.626885E-05 4.481615E-05
1.7246E+12 -8.11469E+11 5.331945E+12 1.00253E+12 3.20893E-05 -1.509848E-05
2.329E+12 -6.5611E+11 1.002534E+12 1.85687E+12 4.333442E-05 -1.220781E-05

V3 V4
-0.0560186 1 0 9
-0.0036041 2 0 10
0 8 0 11
-0.0151813 3 -0.015181255 3
-0.0028537 4 -0.002853733 4
0 12 0 12

Local coordinate
u3 u4
-0.0560186 1 0 9
0 2 0 10
0.00360408 8 0 11
-0.0151813 3 -0.015181255 3
0 4 0.0028537326 4
0.00285373 12 0 12
member end force (global)
 F1 F2 F3
7.9140050015 5 170.3875 5 75.59159975 1
-20.08658665 6 4.499945 6 91.12438012 2
F=TtQ= 19.628785027 7 541.6689 7 0 8
-7.914005001 3 -120.388 1 124.4084003 3
-43.22545336 4 -4.49995 2 104.1428219 4
35.769249983 12 476.0439 8 0 12

178.301538 5
-15.586641 6
561.29765 7
476.043865 8
71.2049092 9
0 10
-194.63468 11
-93.240432 12
Member Load

Fixed End force In Local Coordinat

Qsb 0 5
Qtb 0 6
Qmb = 0 7
Qse 0 3
Qte 0 4
Qme 0 12

-693.302955 0 3684.956
2773.2118201 -541.4159292 -9826.55
-2426.560343 -618.7610619 8598.23
693.30295504 0 -3684.96
2773.2118201 541.4159292 -19653.1
-2426.560343 618.76106195 17196.46

Fixed End force In Local Coordinat

Qsb 34.375 5
Qtb 0 6
Qmb = 65.625 7
Qse 15.625 1
Qte 0 2
Qme 0 8
Matrix In Global Coordinate System

-2427.988338 0 8497.959
0 -1248.571429 0
-8497.959184 0 19828.57
2427.9883382 0 -8497.96
0 1248.5714286 0
-8497.959184 0 39657.14

Fixed End force In Local Coordinat

Qsb 100 1
Qtb 0 2
Qmb = 0 8
Qse 100 3
Qte 0 4
Qme 0 12

Matrix In Global Coordinate System

-1626.5625 0 6506.25
6506.25 0 -17350
0 -1092.5 0
1626.5625 0 -6506.25
6506.25 0 -34700
0 1092.5 0

Fixed End force In Local Coordinat

 Qsb 34.345 9
Qtb 0 10
Qmb = 65.625 11
Qse 15.625 3
Qte 0 4
Qme 0 12

Matrix In Global Coordinate System

-2427.988338 0 8497.959
0 0 0
8497.9591837 0 -19828.6
2427.9883382 0 -8497.96
0 0 0
8497.9591837 0 -39657.1
d End Forces Joint Load
1 0 1
2 0 2
3 0 3
4 0 4

verse Of S
3.20893E-05 4.333442E-05
-1.509848E-05 -1.220781E-05
0.000099208 1.865349E-05
1.865349E-05 3.454949E-05
 F4
71.204909223 9
0 10
-194.6346823 11
-21.23490922 3
0 4
-129.0096823 12