summaryrefslogtreecommitdiff
path: root/1898/CH10/EX10.1/Ex10_1.sce
blob: 2b5ba7209da058603257ef8017590e7f70450fa3 (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
clear all; clc;

disp("Scilab Code Ex 10.1 : ")

//Given:
ep_x = 500; //Normal Strain
ep_y = -300; //Normal Strain
gamma_xy = 200; //Shear Strain
theta = 30*(%pi/180);
theta = theta*-1;

ep_x_new = ((ep_x+ep_y)/2) + ((ep_x - ep_y)/2)*cos(2*theta) + (gamma_xy/2)*sin(2*theta);

gamma_xy_new = -((ep_x - ep_y)/2)*sin(2*theta) + (gamma_xy/2)*cos(2*theta);
gamma_xy_new = 2*gamma_xy_new;

phi = 60*(%pi/180);
ep_y_new = (ep_x+ep_y)/2 + ((ep_x - ep_y)/2)*cos(2*phi) + (gamma_xy/2)*sin(2*phi);

//Display:


printf('\n\nThe equivalent strain acting on the element in the x plain oriented at 30 degrees clockwise    = %1.1f *10^-6',ep_x_new);
printf('\nThe equivalent strain acting on the element in the y plain oriented at 30 degrees clockwise    = %1.1f *10^-6',ep_y_new);
printf('\nThe equivalent shear strain acting on the element                                              = %1.0f *10^-6',gamma_xy_new);

//----------------------------------------------------------------------END--------------------------------------------------------------------------------