Member A: Start Node 2 End Node 3

I (in^4) 875 E (ksi) 29000 Length (in) 180

EI 25375000

v2 θ2 v3 θ3
K v2 52.2119341564 4699.074 -52.211934156 4699.074
θ2 4699.07407407 563888.9 -4699.0740741 281944.4
v3 -52.2119341564 -4699.074 52.2119341564 -4699.074
θ3 4699.07407407 281944.4 -4699.0740741 563888.9

Member A: Start Node 2 End Node 3

I (in^4) 1750 E (ksi) 29000 Length (in) 240
EI 50750000

v2 θ2 v3 θ3
K v2 44.0538194444 5286.458 -44.053819444 5286.458
θ2 5286.45833333 845833.3 -5286.4583333 422916.7
v3 -44.0538194444 -5286.458 44.0538194444 -5286.458
θ3 5286.45833333 422916.7 -5286.4583333 845833.3

Member 1: Start Node 1 End Node 2

I (in^4) 1 E (ksi) 1 Length (in) 16
EI 1

v1 θ1 v2 θ2
K v1 0.0029296875 0.023438 -0.0029296875 0.023438
θ1 0.0234375 0.25 -0.0234375 0.125
v2 -0.0029296875 -0.023438 0.0029296875 -0.023438
θ2 0.0234375 0.125 -0.0234375 0.25

Member 2: Start Node 2 End Node 3

I (in^4) 1 E (ksi) 1 Length (in) 12
EI 1

v2 θ2 v3 θ3
K v2 0.00694444444 0.003472 -0.0069444444 0.003472
θ2 0.00347222222 0.333333 -0.0034722222 0.166667
v3 -0.00694444444 -0.003472 0.00694444444 -0.003472
θ3 0.00347222222 0.166667 -0.0034722222 0.333333
Member 2: Start Node 3 End Node 4
I (in^4) 1 E (ksi) 1 Length (in) 8
EI 1

v3 θ3 v4 θ4
K v3 0.0234375 0.011719 -0.0234375 0.011719
θ3 0.01171875 0.5 -0.01171875 0.25
v4 -0.0234375 -0.011719 0.0234375 -0.011719
θ4 0.01171875 0.25 -0.01171875 0.5

v1 θ1 v2 θ2 v3 θ3 v4
v1 0.00293 0.0234375 -0.00293 0.0234375 0 0 0
θ1 0.023438 0.25 -0.023438 0.125 0 0 0
v2 -0.00293 -0.0234375 0.009874 -0.0199652778 -0.006944 0.003472 0
θ2 0.023438 0.125 -0.019965 0.58333333333 -0.003472 0.166667 0
v3 0 0 -0.006944 -0.0034722222 0.030382 0.008247 -0.023438
θ3 0 0 0.003472 0.16666666667 0.008247 0.833333 -0.011719
v4 0 0 0 0 -0.023438 -0.011719 0.023438
θ4 0 0 0 0 0.011719 0.25 -0.011719

Simplifying Global Matrix

θ2 θ3
θ2 0.583333 0.16666666667
θ3 0.166667 0.83333333333

Member 1: Start Node 1 End Node 2

I (in^4) 510 E (ksi) 29000 Length (in) 240
EI 14790000

v1 θ1 v2 θ2
K v1 12.8385416667 1540.625 -12.838541667 1540.625
θ1 1540.625 246500 -1540.625 123250
v2 -12.8385416667 -1540.625 12.8385416667 -1540.625
θ2 1540.625 123250 -1540.625 246500

Member 2: Start Node 2 End Node 3

I (in^4) 510 E (ksi) 29000 Length (in) 180
EI 14790000
 v2 θ2 v3 θ3
K v2 30.4320987654 2738.889 -30.432098765 2738.889
θ2 2738.88888889 328666.7 -2738.8888889 164333.3
v3 -30.4320987654 -2738.889 30.4320987654 -2738.889
θ3 2738.88888889 164333.3 -2738.8888889 328666.7

Member 3: Start Node 3 End Node 4

I (in^4) 510 E (ksi) 29000 Length (in) 240
EI 14790000

v3 θ3 v4 θ4
K v3 12.8385416667 1540.625 -12.838541667 1540.625
θ3 1540.625 246500 -1540.625 123250
v4 -12.8385416667 -1540.625 12.8385416667 -1540.625
θ4 1540.625 123250 -1540.625 246500

v1 θ1 v2 θ2 v3 θ3 v4
v1 12.83854 1540.625 -12.83854 1540.625 0 0 0
θ1 1540.625 246500 -1540.625 123250 0 0 0
v2 -12.83854 -1540.625 43.27064 1198.26388889 -30.4321 2738.889 0
θ2 1540.625 123250 1198.264 575166.666667 -2738.889 164333.3 0
v3 0 0 -30.4321 -2738.8888889 43.27064 -1198.264 -12.83854
θ3 0 0 2738.889 164333.333333 -1198.264 575166.7 -1540.625
v4 0 0 0 0 -12.83854 -1540.625 12.83854
θ4 0 0 0 0 1540.625 123250 -1540.625

Simplifying Global Matrix

θ2 θ3
θ2 575166.7 164333.333333
θ3 164333.3 575166.666667

Identifying Load vector

0 v1 15 k -15
P 0 θ1 Pfem 900 k-in P -900
0 v2 15 k -15
0 θ2 -900 k-in 900
0 v3 15 k -15
0 θ3 600 k-in -600
0 v4 15 k -15
 0 θ4 -600 k-in 600

P 900 θ2
-600 θ3

Solving for displacements:

d=Inv(K)*P using mmult(minverse(K)*P)

0.002028 θ2 0 v1
d -0.001623 θ3 d 0 θ1
0 v2
0.00202839757 θ2
0 v3
-0.0016227181 θ3
0 v4
0 θ4

Solving for Reactions:

R=Kd-P

3.125 18.125
250 1150
Kd -2.013889 R= 12.9861111111
900 0
-3.611111 11.3888888889
-600 0
2.5 17.5
-200 -800
θ4
0
0
0
0
0.011719
0.25
-0.011719
0.5
θ4
0
0
0
0
1540.625
123250
-1540.625
246500
Member 1: Start Node 1 End Node 2
I (in^4) 1 E (ksi) 1 Length (in) 3
EI 1

v1 θ1 v2 θ2
K v1 0.444444444444 0.666667 -0.4444444444 0.666667
θ1 0.666666666667 1.333333 -0.6666666667 0.666667
v2 -0.44444444444 -0.666667 0.44444444444 -0.666667
θ2 0.666666666667 0.666667 -0.6666666667 1.333333

Member 2: Start Node 2 End Node 3

I (in^4) 1 E (ksi) 2 Length (in) 4
EI 2

v2 θ2 v3 θ3
K v2 0.375 0.75 -0.375 0.75
θ2 0.75 2 -0.75 1
v3 -0.375 -0.75 0.375 -0.75
θ3 0.75 1 -0.75 2

v1 θ1 v2 θ2 v3 θ3
v1 0.444444 0.666666666667 -0.444444 0.66666666667 0 0
θ1 0.666667 1.333333333333 -0.666667 0.66666666667 0 0
v2 -0.444444 -0.66666666667 0.819444 0.08333333333 -0.375 0.75
θ2 0.666667 0.666666666667 0.083333 3.33333333333 -0.75 1
v3 0 0 -0.375 -0.75 0.375 -0.75
θ3 0 0 0.75 1 -0.75 2

Simplifying Matrix
 θ2 θ3
θ2 3.333333 1
θ3 1 2

Identifying load vector

Pfem
Member 1 Member 2
v1 59.25926 v2 100
θ1 35.55556 θ2 66.66667
v2 20.74074 v3 100
θ2 -17.77778 θ3 -66.66667

0 v1 59.2592592593 v1 -59.25926
Pj 0 θ1 35.5555555556 θ1 -35.55556
40 v2 pfem 120.740740741 v2 P= -80.74074
-10 θ2 48.8888888889 θ2 -58.88889
0 v3 100 v3 -100
0 θ3 -66.666666667 θ3 66.66667

P -58.88889 θ2
66.66667 θ3

Solving for displacements:

d=Inv(K)*P using mmult(minverse(K)*P)

0 v1
d= -32.54902 θ2 d= 0 θ1
49.60784 θ3 0 v2
-32.549019608 θ2
0 v3
49.6078431373 θ3

Solving for Reactions

R=Kd-P

kd= -21.69935 v1 37.559912854

-21.69935 θ1 13.8562091503
34.49346 v2 R= 115.234204793
-58.88889 θ2 0
-12.79412 v3 87.2058823529
66.66667 θ3 0
v1
θ1
v2
θ2
v3
θ3