summaryrefslogtreecommitdiff
path: root/3289/CH8/EX8.1/Ex8_1.sce
blob: c5999edcfa07b0f200c8e0abe9b5c43508b09be3 (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
clc

di=0.3 //m
de=0.4 //m
v=0.3
sigmathetamax=250*10^6 //Pa
p0=0
pi=0

//solution a:
a=0.15
b=0.2
r=a
//sigmathetamax=pi*((b^2+a^2)/(b^2-a^2))
pi=sigmathetamax*((b^2-a^2)/(b^2+a^2))
disp(pi,"in Pa is= ")

//solution b:
r=a
//sigmathetamax=-2*p0*(b^2/(b^2-a^2))
p0=-(-sigmathetamax)*((b^2-a^2)/(2*b^2))
disp(p0,"in Pa is= ")

//solution c:
u=((a^3*pi)/(b^2-a^2))*(0.7+1.3*(b^2/a^2))
disp(u,"in per E meter is= ")
sigmaz=(pi*a^2-p0*b^2)/(b^2-a^2)
disp(sigmaz,"for longitudinal stress is")