blob: 99d992598fde06fda0d1983c9a2fa79e2e4864cc (
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
|
clear
//
//
//Initilization of Variables
//strains
e_A=500 //microns
e_B=250 //microns
e_C=-150 //microns
E=2*10**5 //N/mm**2 //Modulus of Elasticity
mu=0.3 //Poissons ratio
theta=45 //Degrees
//Calculations
e_x=500
e_A=500
e_45=250
e_B=250
e_y=-150
e_C=-150
//e_45=(e_x+e_y)*2**-1+(e_x-e_y)*2**-1*cos(2*theta)+rho_x_y*2**-1*sin(2*theta)
//After sub values and further simplifying we get
rho_x_y=(e_45-(e_x+e_y)*2**-1-(e_x-e_y)*2**-1*cos(2*theta*%pi*180**-1))*(sin(2*theta*%pi*180**-1))**-1*2
//Principal strains are given by
e1=(e_x+e_y)*2**-1+(((e_x-e_y)*2**-1)**2+(rho_x_y*2**-1)**2)**0.5 //microns
e2=(e_x+e_y)*2**-1-(((e_x-e_y)*2**-1)**2+(rho_x_y*2**-1)**2)**0.5 //microns
//Principal Stresses
sigma1=E*(e1+mu*e2)*(1-mu**2)**-1*10**-6 //N/mm**2
sigma2=E*(e2+mu*e1)*(1-mu**2)**-1*10**-6 //N/mm**2
//Result
printf("\n Principal Strains are:e1 %0.2f N/mm**2",e1)
printf("\n :e2 %0.2f N/mm**2",e2)
printf("\n Principal Stresses are:sigma1 %0.2f N/mm**2",sigma1)
printf("\n :sigma2 %0.2f N/mm**2",sigma2)
|
