diff options
Diffstat (limited to '3588/CH2/EX2.6/EX2_6.sce')
-rw-r--r-- | 3588/CH2/EX2.6/EX2_6.sce | 43 |
1 files changed, 43 insertions, 0 deletions
diff --git a/3588/CH2/EX2.6/EX2_6.sce b/3588/CH2/EX2.6/EX2_6.sce new file mode 100644 index 000000000..416b7c1d5 --- /dev/null +++ b/3588/CH2/EX2.6/EX2_6.sce @@ -0,0 +1,43 @@ +//Clearing Console +clc +clear + +//Node 1 Displacement +U1=0 + +//Stiffness of Springs +k1=4 +k2=6 +k3=3 + +//Nodal Forces +F2=-30 +F3=0 +F4=50 + +//varible decleration +K=zeros(3,3) + +//Constructing Stiffness and Force matrices +K(1,1)=k1 + 2*k2 +K(1,2)=-2*k2 +K(2,1)=-2*k2 +K(2,2)=2*k2 +k3 +K(2,3)=-k3 +K(3,2)=-k3 +K(3,3)=k3 + +F(1,1)=F2 +F(2,1)=F3 +F(3,1)=F4 + +//Solving for Nodal Displacements U2, U3 and U4 +U=linsolve(K,-F) //K*U=F (equlibrium equation) + +//Solving for Nodal force F1 +F1=-4*U(1,1) + +//Printing Results +printf('\nResults\n') +printf('\nNodal displacements \nU1=%fmm \nU2=%fmm \nU3=%fmm \nU4=%fmm\n',U1,U(1,1),U(2,1),U(3,1)) +printf('\nNodal Force F1=%fN\n',F1) |