diff options
Diffstat (limited to '3588/CH2')
-rw-r--r-- | 3588/CH2/EX2.1/EX2_1.sav | bin | 0 -> 25720 bytes | |||
-rw-r--r-- | 3588/CH2/EX2.1/EX2_1.sce | 35 | ||||
-rw-r--r-- | 3588/CH2/EX2.6/EX2_6.sav | bin | 0 -> 26264 bytes | |||
-rw-r--r-- | 3588/CH2/EX2.6/EX2_6.sce | 43 |
4 files changed, 78 insertions, 0 deletions
diff --git a/3588/CH2/EX2.1/EX2_1.sav b/3588/CH2/EX2.1/EX2_1.sav Binary files differnew file mode 100644 index 000000000..3f7947524 --- /dev/null +++ b/3588/CH2/EX2.1/EX2_1.sav diff --git a/3588/CH2/EX2.1/EX2_1.sce b/3588/CH2/EX2.1/EX2_1.sce new file mode 100644 index 000000000..a9358fcb1 --- /dev/null +++ b/3588/CH2/EX2.1/EX2_1.sce @@ -0,0 +1,35 @@ +//Clearing Console +clc +clear + +//Node 1 Displacement +U1=0 +//Stiffness of Springs +K1=50 +K2=75 +//Nodal Forces +F2=75 +F3=75 +//varible decleration +K=[] +F=[] +U=[] + +//Constructing Stiffness and Force matrices +K(1,1)=K1+K2 +K(1,2)=-K2 +K(2,1)=-K2 +K(2,2)=K2 +F(1,1)=F2 +F(2,1)=F3 + +//Solving for Nodal Displacements U2 and U3 +U=linsolve(K,-F) //K*U=F (equlibrium equation) + +//Solving for Nodal force F1 +F1=-50*U(1,1) + +//Printing Results +printf('\nResults\n') +printf('\nNodal displacements \nU1=%fin \nU2=%fin \nU3=%fin\n',U1,U(1,1),U(2,1)) +printf('\nNodal Force F1=%flb\n',F1) diff --git a/3588/CH2/EX2.6/EX2_6.sav b/3588/CH2/EX2.6/EX2_6.sav Binary files differnew file mode 100644 index 000000000..d4010820b --- /dev/null +++ b/3588/CH2/EX2.6/EX2_6.sav 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) |