summaryrefslogtreecommitdiff
path: root/3864/CH5/EX5.16/Ex5_16.sce
blob: 244888a3d9c4e965cb3fb400f38d136c1f01f949 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
clear
//
//

//Initilization of Variables

L_AC=2 //m //Length of BD,CB,AC
L_BD=2
L_CB=2
F_C=40 //KN //Force at C
F_D=10 //KN Force at D
L=6 //m spna of beam

//EI is constant in this problem

//Calculations

//Let V_A & V_B be the reactions at A & B Respectively
//V_A+V_B=50

//Taking Moment at Pt A
V_B=(F_D*L+F_C*L_AC)*(L_AC+L_CB)**-1
V_A=50-V_B

//Now Taking Moment at distance x from A,M_x
//M_x=15*x-40*(x-2)+35*(x-4)
//EI*(d**2*y/dx**2)=15*x-40*(x-2)+35*(x-4)

//Now Integrating above equation we get
//EI*(dy/dx)=C1+7.5*x**2-20*(x-2)**2+17.5(x-4)**2

//Again Integrating above equation we get
//EI*y=C2+C1*x+2.5*x**2-20*3**-1*(x-2)**3+17.5*(x-4)**3*3**-1

//At
x=0
y=0
//we get
C2=0

//At
x=4 
y=0
//we get
C1=(2.5*4**3-20*3**-1*2**3)*4**-1

//Now Deflection at C
x=2
C1=-26.667
C2=0
y_C=C2+C1*x+2.5*x**3

//Now Deflection at D
C1=-21.667
C2=0
y_D=-26.667*6+2.5*6**3-20*3**-1*4**3+17.5*2**3*3**-1

//Result
printf("\n Deflections Under Loads are:y_D %0.4f  ",y_D)
printf("\n                            :y_C %0.2f  ",y_C)