Análise Dinâmica Do Mecanismo de Uma Prensa Mecânica
Análise Dinâmica Do Mecanismo de Uma Prensa Mecânica
Análise Dinâmica Do Mecanismo de Uma Prensa Mecânica
Coordenação de Mecânica
Curso Superior de Eng. Mecânica
Dinâmica das Máquinas
Prof. Marco Aurélio B. C. Badan
Aluno: Olair Rosa de Souza Júnior
PRENSA MECÂNICA
1. Parâmetros de Projeto
N° de matrícula: 20132010970138
1200−138
𝐿2 = 3
⟹ 𝐿2 = 354𝑚𝑚
𝐿3 = 2 × 354 ⟹ 𝐿3 = 708𝑚𝑚
𝜔2 = 1 + 3 + 8 + 18 ⟹ 𝜔2 = +30𝑟𝑝𝑚
𝐹𝑝 = 100 × 30 ⟹ 𝐹𝑝 = 3000𝑁
ρ = 7870 kg/m3 (Aço 1020)
µ=0,15
2. Análise de Posição
⃗⃗⃗⃗2 + 𝑅
𝑅 ⃗⃗⃗⃗1 − 𝑅
⃗⃗⃗⃗3 = 0
𝐿 . 𝑐𝑜𝑠𝜃2 + 𝐿1 . 𝑐𝑜𝑠𝜃1 − 𝐿3 . 𝑐𝑜𝑠𝜃3 = 0
𝐿2 . 𝑒 𝑖𝜃2 + 𝐿1 . 𝑒 𝑖𝜃1 − 𝐿3 . 𝑒 𝑖𝜃3 = 0 → { 2
𝐿2 . 𝑠𝑖𝑛𝜃2 + 𝐿1 . 𝑠𝑖𝑛𝜃1 − 𝐿3 . 𝑠𝑖𝑛𝜃3 = 0
3. Análise de Velocidade
𝑑 𝑑 𝑑
(𝐿2 . 𝑒 𝑖𝜃2 ) + (𝐿1 . 𝑒 𝑖𝜃1 ) − (𝐿3 . 𝑒 𝑖𝜃3 ) = 0
𝑑𝑡 𝑑𝑡 𝑑𝑡
4. Análise de Aceleração
𝑑 𝑑 𝑑
(𝐿2 . 𝑖. 𝜔2 𝑒 𝑖𝜃2 ) + (𝑉1 . 𝑒 𝑖𝜃1 ) − (𝐿3 . 𝑖. 𝜔3 . 𝑒 𝑖𝜃3 ) = 0
𝑑𝑡 𝑑𝑡 𝑑𝑡
5. Análise de Forças
5.1. Bloco da Prensa (Elo 4)
Para análise de forças, inicia-se com análise da aceleração do centro de massa do
elo. Assumindo o Bloco como um cubo (3 dimensões de valor igual), sua aceleração
vetorial terá o mesmo valor do Elo 1 (Elo Terra). Assim,
⃗⃗⃗⃗⃗⃗
𝑅𝑔4 = 𝐿1 . 𝑒 𝑖𝛩1
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑅 .
= 𝐿1 . 𝑒 𝑖𝛩1 → ⃗⃗⃗
𝑔4
𝑑𝑡
𝑉4 = 𝑉1 . 𝑒 𝑖𝜃1
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑉 ..
= 𝐿1 . 𝑒 𝑖𝛩1 → ⃗⃗⃗⃗
𝑔4
𝑑𝑡
𝐴4 = 𝐴1 . 𝑒 𝑖𝜃1
Portanto,
𝐴 . 𝑐𝑜𝑠𝜃1 𝑎𝑔4𝑥 = 𝐴1 . 𝑐𝑜𝑠𝜃1 = 0
𝐴𝑔4 = { 1
⃗⃗⃗⃗⃗⃗ →
𝐴1 . 𝑠𝑒𝑛𝜃1 𝑎𝑔4𝑦 = 𝐴1 . 𝑠𝑒𝑛𝜃1 = −𝐴1
∑ 𝐹𝑥 = 𝑚4 . 𝑎𝑔4𝑥
⃗⃗⃗⃗⃗⃗
𝑅𝑔3 = 𝐿2 . 𝑒 𝑖𝜃2 + 𝐿𝑔3 . 𝑒 𝑖𝜃3
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑅 𝑔3
𝑑𝑡
= 𝑖. 𝐿2 . 𝜔2 . 𝑒 𝑖𝜃2 + 𝑖. 𝐿𝑔3 . 𝜔3 . 𝑒 𝑖𝜃3
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑉 𝑔3
𝑑𝑡
= −𝐿2 . 𝜔22 . 𝑒 𝑖𝛩2 − 𝐿𝑔3 . 𝜔32 . 𝑒 𝑖𝛩3 + 𝑖. 𝛼3 . 𝑒 𝑖𝜃2
Portanto,
Análise de Momentos
∑ 𝑀 = 𝐼3 . 𝛼3
−𝐹34𝑥 . 𝑅34𝑦 + 𝐹34𝑦 . 𝑅34𝑥 + 𝐹23𝑥 . 𝑅23𝑦 − 𝐹23𝑦 . 𝑅23𝑥 = 𝐼3 . 𝛼3
⃗⃗⃗⃗⃗⃗
𝑅𝑔2 = 𝐿2 . 𝑒 𝑖𝜃2
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑅 𝑔2
𝑑𝑡
= 𝑖. 𝐿2 . 𝜔2 . 𝑒 𝑖𝜃2
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑉 𝑔2
𝑑𝑡
= −𝐿2 . 𝜔22 . 𝑒 𝑖(𝜃2 +180)
Análise de Momentos
∑ 𝑀 = 𝐼2 . 𝛼2 𝛼2 = 0
−𝐹23𝑥 . 𝑅32𝑦 + 𝐹23𝑦 . 𝑅32𝑥 + 𝐹12𝑥 . 𝑅12𝑦 − 𝐹12𝑦 . 𝑅12𝑥 + 𝑇2 = 𝐼2 . 𝛼2 = 0
5.4. Matrizes
Considerando a multiplicação matricial
𝐴 × 𝐵−1 = 𝐶 ,
temos a matriz A de coeficiente, a matriz B de incógnitas e a matriz C de valores
da igualdade das equações de força e momento. Assim,
1 −1 0 0 0 0 0 0
0 0 1 −1 0 0 0 0
𝑅12𝑥 −𝑅32𝑦 −𝑅12𝑦 𝑅32𝑦 1 0 0 0
0 1 0 0 0 −1 0 0
A= 0 0 0 1 0 0 −1 0
0 𝑅23𝑦 0 −𝑅23𝑥 0 −𝑅43𝑦 𝑅43𝑥 0
0 0 0 0 0 1 0 −1
{ 0 0 0 0 0 0 1 𝜇
B = A x C.
6. Resultados
6.748D+00
Insira os dados do material da Biela (Elo 3)
Digite a massa específica [kg/m^3] do elo 3: 7870
Digite a largura do elo 3 em mm: 90
Digite a altura do elo 3 em mm: 90
Massa do elo 3 [kg]:
4.513D+01
Momento de inércia do elo 3 [m^4]:
1.916D+00
Insira os dados do material da Manivela (Elo 2)
Digite a massa específica [kg/m^3] do elo 2: 7870
Digite a largura do elo 2 em mm: 115
Digite a altura do elo 2 em mm: 115
Massa do elo 2 [kg]:
3.684D+01
Momento de inércia do elo 3 [m^4]:
4.254D-01
Torque Máximo em 2 [N.m]:
1.438D+07
Força máxima de 1 em 2 [N]:
2.711D+09
Força máxima de 3 em 2 [N]:
1.553D+07 .
𝐹 2,71.109
𝜎𝑁𝑂𝑅𝑀𝐴𝐿 = 𝐴→ 𝜎𝑁𝑂𝑅𝑀𝐴𝐿 = 0,04425
= 61,42𝑀𝑃𝑎
Assim, o elo 2 suporta a flexão, mas não suporta a tensão normal, sendo
necessárias alterações no projeto quanto ao dimensionamento ou escolha do material.
7. Gráficos
clear;clc;
//Pedindo dados do usuário
mat=input("Digite os três ultimos digitos da sua matricula: ");
printf("Assumindo um numero de matricula 20132010970138, \n");
printf("os três ultimos digitos serão 1, 3 e 8.\n");
printf("Informe os três ultimos digitos da sua matrícula na sequência\n\n");
Mat1=input("Digite o antepenultimo numero: ");
Mat2=input("Digite o penultimo numero: ");
Mat3=input("Digite o ultimo numero: ");
//Análise de Posição
Gama3=acos((L2/L3)*cos(Teta2)); //ângulo de construção Teta3
Teta3=-Gama3
L1=(L2*sin(Teta2)-L3*sin(Teta3));
//Análise de Velocidades
W3=((L2*W2*sin(Teta2))/(L3*sin(Teta3)));
V1=((L2*W2*cos(Teta2))-(L3*W3*cos(Teta3)))
V2=sqrt(((L2*W2*sin(Teta2))^2)+((L2*W2*cos(Teta2))^2));
//Análise de Acelerações
Alfa2=0;
Alfa3=(((L2*(W2^2)*cos(Teta2))-(L3*(W3^2)*cos(Teta3)))/(L3*sin(Teta3)))
A1=((L3*(W3^2)*sin(Teta3))-(L2*(W2^2)*sin(Teta3))-(L3*Alfa3*cos(Teta3)));
A2=sqrt(((L2*(W2^2)*cos(Teta2))^2)+((L2*(W2^2)*sin(Teta2))^2));
//Análise de Forças
////Matrizes de Forças
g=9.81;
mi=0.15;
P2y=(Massa2*g);
P3y=(Massa3*g);
P4y=(Massa4*g);
T2Max=0;
F12xMax=0;
F12yMax=0;
F32xMax=0;
F32yMax=0;
for i=1:360
MatFator = [1 -1 0 0 0 0 0 0 ;
0 0 1 -1 0 0 0 0 ;
R12y(1,i) -R32y(1,i) -R12x(1,i) R32x(1,i) 1 0 0 0 ;
0 1 0 0 0 -1 0 0 ;
0 0 0 1 0 0 -1 0 ;
0 R23y(1,i) 0 -R23x(1,i) 0 -R34y(1,i) R34x(1,i) 0 ;
0 0 0 0 0 1 0 -1 ;
0 0 0 0 0 0 1 mi];
MatIgual=[(Massa2*Ag2x(1,i));
(Massa2*Ag2y(1,i)+P2y);
J2*Alfa2;
(Massa3*Ag3x(1,i));
(Massa3*Ag3y(1,i)+P3y);
J3*Alfa3;
Massa4*Ag4x;
(Massa4*Ag4y(1,i)+P4y-Fp)];
MatIncog=MatFator*MatIgual;
Cont(1,i)=[MatIncog(5,1)];
if MatIncog(5,1)>T2Max then
T2Max=MatIncog(5,1);
end
if MatIncog(1,1)>F12xMax then
F12xMax=MatIncog(1,1);
end
if MatIncog(3,1)>F12yMax then
F12yMax=MatIncog(3,1);
end
if MatIncog(2,1)>F32xMax then
F32xMax=MatIncog(2,1);
end
if MatIncog(4,1)>F32yMax then
F32yMax=MatIncog(4,1);
end
end
F32Max=sqrt((F32xMax^2)+(F32yMax^2));
F12Max=sqrt((F12xMax^2)+(F12yMax^2));
printf("\nTorque Máximo em 2 [N.m]: ");
disp(T2Max)
printf("\nForça máxima de 1 em 2 [N]: ");
disp(F12Max);
printf("\nForça máxima de 3 em 2 [N]: ");
disp(F32Max);
scf(1)
plot(Teta2,L1)
xlabel("Teta 2 [rad]","fontsize", 6, "color", "blue")
ylabel("Tamanho elo terra [mm]","fontsize", 6, "color", "red")
scf(2)
plot(Teta2,V1)
xlabel("Teta 2 [rad]","fontsize", 6, "color", "blue")
ylabel("Velocidade elo terra [mm/s]","fontsize", 6, "color", "blue")
scf(3)
plot(Teta2,A1)
xlabel("Teta 2 [rad]","fontsize", 6, "color", "blue")
ylabel("Aceleração elo terra [mm/s^2]","fontsize", 6, "color", "magenta")
scf(4)
plot(Teta2,Cont)
xlabel("Teta 2 [rad]","fontsize", 6, "color", "blue")
ylabel("Torque Elo 2 [mm.s]","fontsize", 6, "color", "red")
scf(5)
plot(Teta2,V2)
xlabel("Teta 2 [rad]","fontsize", 6, "color", "blue")
ylabel("Velocidade elo 2 [mm/s]","fontsize", 6, "color", "blue")
scf(6)
plot(Teta2,A2)
xlabel("Teta 2 [rad]","fontsize", 6, "color", "blue")
ylabel("Aceleração elo 2 [mm/s^2]","fontsize", 6, "color", "magenta")
9. Referencial Bibliográfico