summaryrefslogtreecommitdiff
path: root/3864/CH2/EX2.45/Ex2_45.sce
blob: bcb78b03e815cd19be05210bfb5fc5fa7768d396 (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
61
62
63
64
65
66
clear
//
//

//Initilization of Variables

//Bar-A
d1=30 //mm //Diameter of bar1
L=600 //mm //length of bar1

//Bar-B
d2=30 //mm //Diameter of bar2
d3=20 //mm //Diameter of bar2
L2=600 //mm //length of bar2

//Calculations

//Area of bar-A
A1=%pi*4**-1*d1**2

//Area of bar-B
A2=%pi*4**-1*d2**2
A3=%pi*4**-1*d3**2

//let SE be the Strain Energy
//Strain Energy stored in Bar-A
//SE=p**2*(2*E)**-1*V
//After substituting values and simolifying further we get
//SE=P**2*E**-1*0.4244

//Strain Energy stored in Bar-B
//SE2=p1**2*V1*(2*E)**-1+p2**2*V2*(2*E)**-1
//After substituting values and simolifying further we get
//SE2=0.6897*P**2*E**-1

//Let X be the ratio of SE in Bar-B and SE in Bar-A
X=0.6897*0.4244**-1

//Part-2

//When Max stress is produced is same:Let p be the max stress produced

//Stress in bar A is p throughout 
//In bar B:stress in 20mm dia.portion=p2=p

//Stress in 30 mm dia.portion
//p1=P*A2*A3**-1
//After substituting values and simolifying further we get
//p1=4*9**-1*p

//Strain Energy in bar A
//SE_1=p**2*(2*E)**-1*A1*L1
//After substituting values and simolifying further we get
//SE_1=67500*p**2*%pi*E**-1

//Strain Energy in bar B
//SE_2=p1**2*V1*(2*E)**-1+p2**2*V2*(2*E)**-1
//After substituting values and simolifying further we get
//SE_2=21666.67*%pi*p**2*E**-1

//Let Y be the Ratio of SE in bar B and SE in bar A
Y=21666.67*67500**-1

//result
printf("\n Gradually applied Load is %0.2f  ",X)
printf("\n Gradually applied Load is %0.2f  ",Y)