clc

clear

all

N=0.5;

n=2600;

omega=n*2*pi/60;

z1=12;

z2=36;

z3=14;

z2=36;

z3=14;

z4=28;

i=z2/z1;

beta=5*pi/180;

alfa=20*pi/180;

m=1.25*1e-3;

d_z1=z1*m

d_z2=z2*m;

d_z3=z3*m;

d_z4=z4*m;

a=0.03;

b=0.04;

c=0.03;

d_z1 =

0.0150

M_2=N/omega*i*1e3

M_2 =

5.5092

P_o2=M_2/d_z2*2

% Sila wzdluzna

P_w2=P_o2*tan(beta)

% Sila Promieniowa

P_r2=P_o2*tan(alfa)/cos(beta)

P_o2 =

244.8538

P_w2 =

21.4219

P_r2 =

89.4599

P_o3=M_2/d_z3*2

% Sila Promieniowa

P_r3=P_o3*tan(alfa)/cos(beta)

P_o3 =

629.6240

P_r3 =

230.0397

R_By=(a*P_r2 - P_r3*(a+b) + P_w2*d_z2/2)/ (a+b+c)

% Sila w podporze Ay, suma momentow wzg punktu B

R_Ay=( - P_r3*c + P_r2*(c+b) - P_w2*d_z2/2)/(a+b+c)

% Sprawdzenie

spr=round( R_By + R_Ay - P_r2+P_r3)

if

(spr == 0)

disp(

'Dobrze policzone warunki'

)

else

disp(

'Zle policzone warunki'

)

end

%siła w podporze A wzdluz osi x

R_Ax=P_w2

R_By =

-129.3699

R_Ay =

-11.2099

-11.2099

spr =

0

Dobrze policzone warunki

R_Ax =

21.4219

R_Bz=(P_o2*a+P_o3*(a+b))/(a+b+c)

R_Az=(P_o2*(b+c)+P_o3*(c))/(a+b+c)

% Sprawdzenie

spr=round(-R_Bz-R_Az+P_o2+P_o3)

if

(spr == 0)

disp(

'Dobrze policzone warunki'

)

else

disp(

'Zle policzone warunki'

)

end

R_Bz =

514.1929

R_Az =

360.2848

spr =

0

Dobrze policzone warunki

x1=0:0.01:a

x1=0:0.01:a

Mg_x1_xy=R_Ay*x1

x1 =

0 0.0100 0.0200 0.0300

Mg_x1_xy =

0 -0.1121 -0.2242 -0.3363

x2=a:0.01:(a+b)

Mg_x2_xy= R_Ay*x2 + P_w2*d_z2/2 - P_r2*(x2-a)

x2 =

0.0300 0.0400 0.0500 0.0600 0.0700

Mg_x2_xy =

0.1457 -0.8610 -1.8677 -2.8744 -3.8811

x3=-c:0.01:0

Mg_x3_xy=abs(R_By)*x3

% Przedzial sprawdzajacy

x4=a+b:0.01:(a+b+c)

Mg_x4_xy=R_Ay*x4 +P_w2*d_z2/2- P_r2*(x4-a)+ P_r3*(x4-(a+b))

x3 =

-0.0300 -0.0200 -0.0100 0

Mg_x3_xy =

-3.8811 -2.5874 -1.2937 0

x4 =

0.0700 0.0800 0.0900 0.1000

0.0700 0.0800 0.0900 0.1000

Mg_x4_xy =

-3.8811 -2.5874 -1.2937 0

x1=0:0.01:a

Mg_x1_xz=-abs(R_Az)*x1

x1 =

0 0.0100 0.0200 0.0300

Mg_x1_xz =

0 -3.6028 -7.2057 -10.8085

x2=a:0.01:(a+b)

Mg_x2_xz=-abs(R_Az)*x2+P_o2*(x2-a)

x2 =

0.0300 0.0400 0.0500 0.0600 0.0700

Mg_x2_xz =

-10.8085 -11.9629 -13.1172 -14.2715 -15.4258

x3=-c:0.01:0

Mg_x3_xz=abs(R_Bz)*x3

% Przedzial sprawdzajacy

x4=a+b:0.01:(a+b+c)

Mg_x4_xz=-R_Az*x4 +P_o2*(x4-(a)) + P_o3*(x4-(a+b))

x3 =

x3 =

-0.0300 -0.0200 -0.0100 0

Mg_x3_xz =

-15.4258 -10.2839 -5.1419 0

x4 =

0.0700 0.0800 0.0900 0.1000

Mg_x4_xz =

-15.4258 -10.2839 -5.1419 0

l=[x1,x2,x4]

Mg_XY=[Mg_x1_xy,Mg_x2_xy,Mg_x4_xy]

figure(1)

plot(l,Mg_XY,

'linewidth'

,2);

xlabel(

'Dlugosc walu'

)

ylabel(

'Momenty gnace'

)

l =

Columns 1 through 5

0 0.0100 0.0200 0.0300 0.0300

Columns 6 through 10

0.0400 0.0500 0.0600 0.0700 0.0700

Columns 11 through 13

0.0800 0.0900 0.1000

Mg_XY =

Columns 1 through 5

0 -0.1121 -0.2242 -0.3363 0.1457

Columns 6 through 10

Columns 6 through 10

-0.8610 -1.8677 -2.8744 -3.8811 -3.8811

Columns 11 through 13

-2.5874 -1.2937 0

l=[x1,x2,x4]

Mg_XZ=[Mg_x1_xz,Mg_x2_xz,Mg_x4_xz]

figure(2)

plot(l,Mg_XZ,

'linewidth'

,2);

xlabel(

'Dlugosc walu'

)

ylabel(

'Momenty gnace'

)

l =

Columns 1 through 5

0 0.0100 0.0200 0.0300 0.0300

Columns 6 through 10

Columns 6 through 10

0.0400 0.0500 0.0600 0.0700 0.0700

Columns 11 through 13

0.0800 0.0900 0.1000

Mg_XZ =

Columns 1 through 5

0 -3.6028 -7.2057 -10.8085 -10.8085

Columns 6 through 10

-11.9629 -13.1172 -14.2715 -15.4258 -15.4258

Columns 11 through 13

-10.2839 -5.1419 0

Mg_w=sqrt(Mg_XY.^2+Mg_XZ.^2)

figure(3)

plot(l,Mg_XZ,

'linewidth'

,2);

xlabel(

'Dlugosc walu'

)

ylabel(

'Moment gnacy zastepczy'

)

Mg_w =

Columns 1 through 5

0 3.6046 7.2092 10.8138 10.8095

Columns 6 through 10

11.9938 13.2495 14.5581 15.9065 15.9065

Columns 11 through 13

10.6044 5.3022 0

M_s=M_2

ll=l;

for

(aa=1:length(l))

if

(l(aa)<0.03 || l(aa)>0.07)

ll(aa)=0;

else

ll(aa)=1;

end

M_ss(aa)=M_s*ll(aa);

end

M_s=M_ss;

M_s =

5.5092

k_sj=85

k_so=45

k_sj=85

k_so=45

k_go=80

k_sj =

85

k_so =

45

k_go =

80

%strona 177 osinski

M_s_prim=sqrt(3)/2*k_so/k_sj*M_s

M_s_prim =

Columns 1 through 5

0 0 0 2.5259 2.5259

Columns 6 through 10

2.5259 2.5259 2.5259 2.5259 2.5259

Columns 11 through 13

0 0 0

M_z=sqrt(Mg_w.^2+M_s_prim.^2)

figure(4)

plot(l,M_z,

'linewidth'

,2);

xlabel(

'Dlugosc walu'

)

ylabel(

'Moment gnacy zredukowany'

)

M_z =

M_z =

Columns 1 through 5

0 3.6046 7.2092 11.1049 11.1007

Columns 6 through 10

12.2569 13.4881 14.7756 16.1058 16.1058

Columns 11 through 13

10.6044 5.3022 0

d_teor=(M_z*32/pi/(k_go*1e6)).^(1/3);

figure(6)

plot(l,d_teor/2*1e3,

'linewidth'

,2);

xlabel(

'Dlugosc walu'

)

ylabel(

'Teoretyczny promien walu (mm)'

)