summaryrefslogtreecommitdiff
path: root/1205/CH7/EX7.2/S_7_2.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /1205/CH7/EX7.2/S_7_2.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '1205/CH7/EX7.2/S_7_2.sce')
-rw-r--r--1205/CH7/EX7.2/S_7_2.sce53
1 files changed, 53 insertions, 0 deletions
diff --git a/1205/CH7/EX7.2/S_7_2.sce b/1205/CH7/EX7.2/S_7_2.sce
new file mode 100644
index 000000000..795251925
--- /dev/null
+++ b/1205/CH7/EX7.2/S_7_2.sce
@@ -0,0 +1,53 @@
+clc;
+//Drawing of shear and bending moment diagram
+printf("Given problem is for drawing diagram, this diagram is drawn by step by step manner.\n ");
+F_A=-20;//kN, force applied at A
+F_C=-40;//kN, force applied at C
+AB=2.5;//m, perpendicular distance between A and B
+BC=3;//m, perpendicular distance between C and B
+CD=2;//m, perpendicular distance between C and D
+//By free body of entire beam
+//By sum(m_D)=0
+R_B=-(CD*F_C+(AB+BC+CD)*F_A)/(BC+CD);//kN, Reaction atB
+//By sum(m_A)=0
+R_D=-(BC*F_C-(AB)*F_A)/(BC+CD);//kN, Reaction atB
+//For section 1
+//Applying sum(Fy)=0
+V1=F_A;//kN
+//Applying sum(M1)=0
+M1=V1*0;//kN.m
+
+//For section 2
+//Applying sum(Fy)=0
+V2=F_A;//kN
+//Applying sum(M1)=0
+M2=F_A*AB;//kN.m
+
+//For section 3
+//Applying sum(Fy)=0
+V3=R_B+F_A;//kN
+//Applying sum(M1)=0
+M3=F_A*AB;//kN.m
+
+//For section 4
+//Applying sum(Fy)=0
+V4=R_B+F_A;//kN
+//Applying sum(M1)=0
+M4=F_A*(AB+BC)+R_B*BC //kN.m
+
+//For section 5
+//Applying sum(Fy)=0
+V5=R_B+F_A+F_C;//kN
+//Applying sum(M1)=0
+M5=F_A*(AB+BC)+R_B*BC//kN.m
+
+//For section 6
+//Applying sum(Fy)=0
+V6=R_B+F_A+F_C;//kN
+//Applying sum(M1)=0
+M6=V6*0//kN.m
+X=[0,2.5,2.5,5.5,5.5,7.5]
+V=[V1,V2,V3,V4,V5,V6];//Shear matrix
+M=[M1,M2,M3,M4,M5,M6];//Bending moment matrix
+plot(X,V);//Shear diagram
+plot(X,M,'r');//Bending moment diagram