4 5 decimal

Originator: L.Weekes&[DATE] &[TIME] - &[FILENAME]

Stiffness Method for Pin-Jointed Frames

nn 4 :

i 1 nn range :

nm 5 :

j 1 nm range :

x 330041 mat m * :

y 303041 mat m * :

CON 11321231442452353 mat :

Area 50050075075050051 mat mm 2^* :

E 20020020020020051 mat kN * mm 2 -^* :

for

j 1 nm range L j el x CON j 3 el el x CON j 2 el el - (2^y CON j 3 el el y CON j 2 el el - (2^+ sqrt : for

L 3 m * 3 m * 4.2426 m * 3 m * 4.2426 m * 51 mat

j 1 nm range K j el E j el Area j el * L j el / 11 - 1 - 122 mat * : for

j 1 nm range θ j el x CON j 3 el el x CON j 2 el el - 0 > y CON j 3 el el y CON j 2 el el - L j el / asin 180 deg * y CON j 3 el el y CON j 2 el el - L j el / asin - if : for

θ deg 18027022518013551 mat

0

OK HERE

j 1 nm range z 12 nn * range t j el 1 z el z 2 CON j 2 el * 1 - ≡ θ j el cos z 2 CON j 2 el * ≡ θ j el sin 0 if if : for for

j 1 nm range z 12 nn * range t j el 2 z el z 2 CON j 3 el * 1 - ≡ θ j el cos z 2 CON j 3 el * ≡ θ j el sin 0 if if : for for

t 00 0.7071 - 0.7071 00000000 0.7071 - 0.7071 0028 mat

Only gives the last transformation matix when j=5 - is there a problem with the nesting?

doing the z loop ok but not the j loop

Alternative ways of programming this?

j 1 nm range temp j el 23 : for

temp 232323232351 mat