blob: c5938e1639910a5bf8a3a9a3a924238fb4af7be7 (
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
|
clear all; clc;
disp("Scilab Code Ex 6.18 : ")
//Given:
M = 12; //kNm
l_bc = 0.2; //m
l_be = 0.4; //m
//Internal Moment Components:
My = (-4/5)*M;
Mz = (3/5)*M;
Iy = (1/12)*(l_be*l_bc^3);
Iz = (1/12)*(l_bc*l_be^3);
//Bending Stress:
sigma_B = (-Mz*1000*(l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;
sigma_B = sigma_B/10^6;
sigma_C = (-Mz*1000*(l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;
sigma_C = sigma_C/10^6;
sigma_D = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;
sigma_D = sigma_D/10^6;
sigma_E = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;
sigma_E = sigma_E/10^6;
//Orientation of Nuetral Axis:
z = (0.45)/(sigma_E + sigma_B);
//theta = -atan(4/3);
tanA = (Iz/Iy)*(-4/3);
alpha = atan(tanA);
alpha = alpha*(180/%pi);
//Display:
printf("\n\nThe normal stress at B = %1.2f MPa',sigma_B);
printf("\nThe normal stress at C = %1.2f MPa',sigma_C);
printf("\nThe normal stress at D = %1.2f MPa',sigma_D);
printf("\nThe normal stress at E = %1.2f MPa',sigma_E);
printf("\nThe orientation of the nuetral axis = %1.1f degrees',alpha);
//------------------------------------------------------------------------END---------------------------------------------------------------------------------------
|
