diff options
author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
---|---|---|
committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /1205/CH7/EX7.5 | |
parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
download | Scilab-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.5')
-rw-r--r-- | 1205/CH7/EX7.5/S_7_5.sce | 49 |
1 files changed, 49 insertions, 0 deletions
diff --git a/1205/CH7/EX7.5/S_7_5.sce b/1205/CH7/EX7.5/S_7_5.sce new file mode 100644 index 000000000..c6370eb0d --- /dev/null +++ b/1205/CH7/EX7.5/S_7_5.sce @@ -0,0 +1,49 @@ +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 ");
+
+w=20;//kN/m, distributed load applied at D to E
+AB=6;//m, perpendicular distance between A and B
+BC=3;//m, perpendicular distance between C and B
+
+F_B=w*AB;//kN, force exerted at AB/2 from A
+
+//By free body of entire beam
+//By sum(m_C)=0
+RA=(F_B*(AB/2+BC))/(AB+BC);//kN, Reaction at A
+
+//By sum(m_A)=0
+RC=(F_B*(AB/2)/(AB+BC));//kN, Reaction at C
+
+//For section 1
+//Applying sum(Fy)=0
+VA=RA;//N, shear force just to right to A
+
+//For section 2
+//Applying sum(Fy)=0
+VB=VA-F_B;//kN, shear force just left to B
+
+//For section 3
+//Applying sum(Fy)=0
+VC=VB;//kN, shear force from B to C
+
+
+//Bending moment at each end is zero
+// Maximum bending moment is at D where V=0
+VD=0;//kN
+
+x=-(VD-VA)/w;//m, location of maximum bending moment
+printf("Maximum bending moment is at D x= %.0f m from A\n",x);
+MA=0;//kN.m
+MD=MA+1/2*VA*x;//kN.m, maximum bending moment is at D
+MB=MD+1/2*VB*(AB-x);//N.m
+MC=MB+VB*BC;//N.m
+
+printf("Maximum bending moment is at MD= %.0fkN. m from A\n",MD);
+X=[0,x,AB,AB+BC];//m,
+V=[VA,VD,VB,VC];//kN,Shear matrix,
+
+plot(X,V);//Shear diagram
+X=[0,x,AB,AB+BC];//m
+M=[MA,MD,MB,MC];//kN.m,Bending moment matrix
+plot(X,M,'r');//Bending moment diagram
|