Análise Dinâmica Do Mecanismo de Uma Prensa Mecânica

Fazer download em docx, pdf ou txt
Fazer download em docx, pdf ou txt
Você está na página 1de 13

Instituto Federal de Goiás - Câmpus Goiânia

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
𝑑𝑡 𝑑𝑡 𝑑𝑡

𝐿2 . 𝑖. 𝜔2 𝑒 𝑖𝜃2 + 𝑉1 . 𝑒 𝑖𝜃1 − 𝐿3 . 𝑖. 𝜔3 . 𝑒 𝑖𝜃3 = 0

−𝐿2 . 𝜔2 . 𝑠𝑖𝑛𝜃2 + 𝑉1 . 𝑐𝑜𝑠𝜃1 + 𝐿3 . 𝜔3 . 𝑠𝑖𝑛𝜃3 = 0


{
𝐿2 . 𝜔2 . 𝑐𝑜𝑠𝜃2 + 𝑉1 . 𝑠𝑖𝑛𝜃1 − 𝐿3 . 𝜔3 . 𝑐𝑜𝑠𝜃3 = 0

4. Análise de Aceleração

𝑑 𝑑 𝑑
(𝐿2 . 𝑖. 𝜔2 𝑒 𝑖𝜃2 ) + (𝑉1 . 𝑒 𝑖𝜃1 ) − (𝐿3 . 𝑖. 𝜔3 . 𝑒 𝑖𝜃3 ) = 0
𝑑𝑡 𝑑𝑡 𝑑𝑡

𝐿2 . 𝑖 2 . 𝜔22 . 𝑒 𝑖𝜃2 + 𝐴1 . 𝑒 𝑖𝜃1 − 𝐿3 . 𝑖. 𝛼3 . 𝑒 𝑖𝜃3 − 𝐿3 . 𝑖 2 . 𝜔32 . 𝑒 𝑖𝜃3 = 0

−𝐿2 . 𝜔22 . 𝑐𝑜𝑠𝜃2 + 𝐴1 . 𝑐𝑜𝑠𝜃1 + 𝐿3 . 𝛼3 . 𝑠𝑖𝑛𝜃3 + 𝐿3 . 𝜔32 . 𝑐𝑜𝑠𝜃3 = 0


{
−𝐿2 . 𝜔22 . 𝑠𝑖𝑛𝜃2 + 𝐴1 . 𝑠𝑖𝑛𝜃1 − 𝐿3 . 𝛼3 . 𝑐𝑜𝑠𝜃3 + 𝐿3 . 𝜔32 . 𝑠𝑖𝑛𝜃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

A análise de Forças no eixo x → +x

∑ 𝐹𝑥 = 𝑚4 . 𝑎𝑔4𝑥

𝐹34𝑥 − 𝐹14𝑥 = 𝑚4 . 𝑎𝑔4𝑥 = 0

Análise de Forças no eixo y ↑ +y


∑ 𝐹𝑦 = 𝑚4 . 𝑎𝑔4𝑦

𝐹34𝑦 ± 𝐹14𝑥 . 𝜇𝑎𝑡 + 𝐹𝑝 − 𝑃4 = 𝑚4 . 𝑎𝑔4𝑦

Como as forças são aplicadas no centro de massa, não gera momento.

5.2. Biela (Bloco 3)

Para a análise da aceleração do centro de massa do elo, foi feito

⃗⃗⃗⃗⃗⃗
𝑅𝑔3 = 𝐿2 . 𝑒 𝑖𝜃2 + 𝐿𝑔3 . 𝑒 𝑖𝜃3
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑅 𝑔3
𝑑𝑡
= 𝑖. 𝐿2 . 𝜔2 . 𝑒 𝑖𝜃2 + 𝑖. 𝐿𝑔3 . 𝜔3 . 𝑒 𝑖𝜃3
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑉 𝑔3
𝑑𝑡
= −𝐿2 . 𝜔22 . 𝑒 𝑖𝛩2 − 𝐿𝑔3 . 𝜔32 . 𝑒 𝑖𝛩3 + 𝑖. 𝛼3 . 𝑒 𝑖𝜃2

Portanto,

−𝐿2 . 𝜔22 . 𝑐𝑜𝑠𝜃2 − 𝐿𝑔3 . 𝜔32 . 𝑐𝑜𝑠𝜃3 − 𝐿𝑔3 . 𝛼3 . 𝑠𝑒𝑛𝜃3 = 𝑎𝑔3𝑥


⃗⃗⃗⃗⃗⃗
𝐴𝑔3 = {
−𝐿2 . 𝜔22 . 𝑠𝑒𝑛𝜃2 − 𝐿𝑔3 . 𝜔32 . 𝑠𝑒𝑛𝜃3 + 𝐿𝑔3 . 𝛼3 . 𝑐𝑜𝑠𝜃3 = 𝑎𝑔3𝑦

Vetores de posição das forças


𝐿𝑔3 . 𝑐𝑜𝑠(𝜃3 + 180) = 𝑅23𝑥
⃗⃗⃗⃗⃗⃗
𝑅23 = {
𝐿𝑔3 . 𝑠𝑒𝑛(𝜃3 + 180) = 𝑅23𝑦
𝐿𝑔3 . 𝑐𝑜𝑠𝜃3 = 𝑅34𝑥
⃗⃗⃗⃗⃗⃗
𝑅34 = {
𝐿𝑔3 . 𝑠𝑒𝑛𝜃3 = 𝑅34𝑦

A análise de Forças no eixo x → +x

∑ 𝐹𝑥 = 𝑚3 . 𝑎𝑔3𝑥 𝐹43𝑥 = −𝐹34𝑥

𝐹34𝑥 + 𝐹23𝑥 = 𝑚3 . 𝑎𝑔3𝑥 −𝐹43𝑥 + 𝐹23𝑥 = 𝑚3 . 𝑎𝑔3𝑥

Análise de Forças no eixo y ↑ +y

∑ 𝐹𝑦 = 𝑚3 . 𝑎𝑔3𝑦 𝐹43𝑦 = −𝐹34𝑦

𝐹43𝑦 + 𝐹23𝑦 − 𝑃3 = 𝑚. 𝑎𝑔4𝑦 −𝐹34𝑦 + 𝐹23𝑦 − 𝑃3 = 𝑚. 𝑎𝑔4𝑦

Análise de Momentos
∑ 𝑀 = 𝐼3 . 𝛼3
−𝐹34𝑥 . 𝑅34𝑦 + 𝐹34𝑦 . 𝑅34𝑥 + 𝐹23𝑥 . 𝑅23𝑦 − 𝐹23𝑦 . 𝑅23𝑥 = 𝐼3 . 𝛼3

5.3. Manivela (Bloco 2)


Para a análise da aceleração do centro de massa do elo, foi feito

⃗⃗⃗⃗⃗⃗
𝑅𝑔2 = 𝐿2 . 𝑒 𝑖𝜃2
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑅 𝑔2
𝑑𝑡
= 𝑖. 𝐿2 . 𝜔2 . 𝑒 𝑖𝜃2
⃗⃗⃗⃗⃗⃗⃗
𝑑𝑉 𝑔2
𝑑𝑡
= −𝐿2 . 𝜔22 . 𝑒 𝑖(𝜃2 +180)

Foi adicionado 180 ao θ2 pois o vetor se encontra ao contrário do adotado na


análise de posição. Portanto,

−𝐿2 . 𝜔22 . 𝑐𝑜𝑠(𝜃2 + 180) = 𝑎𝑔2𝑥


⃗⃗⃗⃗⃗⃗
𝐴𝑔2 = {
−𝐿2 . 𝜔22 . 𝑠𝑒𝑛(𝜃2 + 180) = 𝑎𝑔2𝑦

Vetores de posição das forças


𝐿𝑔2 . 𝑐𝑜𝑠(𝜃2 + 180) = 𝑅32𝑥
⃗⃗⃗⃗⃗⃗
𝑅32 = {
𝐿𝑔2 . 𝑠𝑒𝑛(𝜃2 + 180) = 𝑅32𝑦
𝐿𝑔2 . 𝑐𝑜𝑠𝜃2 = 𝑅12𝑥
⃗⃗⃗⃗⃗⃗
𝑅12 = {
𝐿𝑔2 . 𝑠𝑒𝑛𝜃2 = 𝑅12𝑦

A análise de Forças no eixo x → +x

∑ 𝐹𝑥 = 𝑚2 . 𝑎𝑔2𝑥 𝐹32𝑥 = −𝐹23𝑥

𝐹12𝑥 + 𝐹32𝑥 = 𝑚2 . 𝑎𝑔2𝑥 −𝐹23𝑥 + 𝐹12𝑥 = 𝑚2 . 𝑎𝑔2𝑥

Análise de Forças no eixo y ↑ +y

∑ 𝐹𝑦 = 𝑚2 . 𝑎𝑔2𝑦 𝐹32𝑦 = −𝐹23𝑦

𝐹12𝑦 + 𝐹32𝑦 − 𝑃2 = 𝑚2 . 𝑎𝑔2𝑦 −𝐹23𝑦 + 𝐹12𝑦 − 𝑃3 = 𝑚2 . 𝑎𝑔2𝑦

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 = {𝐹12𝑥 𝐹23𝑥 𝐹12𝑦 𝐹23𝑦 𝑇2 𝐹34𝑥 𝐹34𝑦 𝐹14𝑥


𝑚2 . 𝑎𝑔2𝑥
𝑚2 . 𝑎𝑔2𝑦 + 𝑃2𝑦
𝐼2 . 𝛼2
𝑚3 . 𝑎𝑔3𝑥
C = 𝑚 .𝑎 + 𝑃
3 𝑔3𝑦 3𝑦
𝐼3 . 𝛼3
𝑚4 . 𝑎4𝑥
{ 𝑚 4 . 𝑎4𝑦

Dos dados, obtemos valores vetoriais para as matrizes A e C, para encontrar a


matriz B, faremos

B = A x C.

6. Resultados

A partir da multiplicação das matrizes, escolhem-se os materiais dos elos para


fazer o cálculo. Para todos os elos, foi escolhido o aço 1020. Com o número de
matrícula 20132010970138, foram obtidos os seguintes resultados

Digite os três ultimos digitos da sua matricula: 138


Assumindo um numero de matricula 20132010970138,
os três ultimos digitos serão 1, 3 e 8.
Informe os três ultimos digitos da sua matrícula na sequência
Digite o antepenultimo numero: 1
Digite o penultimo numero: 3
Digite o ultimo numero: 8
L2 será: 354 [mm]
L3 será: 708 [mm]
W2 será: 30 [rad/s]
Fp será: 3000 [N]
Insira os dados do material do Bloco (Elo 4)
Digite a massa específica [kg/m^3] do elo 4: 7870
Digite o comprimento do elo 4 em mm: 95
Digite a largura do elo 4 em mm: 95
Digite a altura do elo 4 em mm: 95
Massa do elo 4 [kg]:

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 .

A partir dos resultados, podemos conferir as forças máximas e determinar se o elo


2 resiste às tensões normal e de flexão.
A carga máxima recebida no vetor normal é no eixo y de valor 2,71.10 9 [N]. A
tensão normal será dada por

𝐹 2,71.109
𝜎𝑁𝑂𝑅𝑀𝐴𝐿 = 𝐴→ 𝜎𝑁𝑂𝑅𝑀𝐴𝐿 = 0,04425
= 61,42𝑀𝑃𝑎

Com o Torque 2 máximo de valor 1,438.107 , a flexão será dada por


𝐿2
.𝑇 177.1,438.107
2 2
𝜎𝐹𝐿𝐸𝑋𝐴𝑂 = 𝐼2
→ 𝜎𝐹𝐿𝐸𝑋𝐴𝑂 = 0,4254.1012
= 598,3307𝑀𝑃𝑎
Assumindo o limite de escoamento de 350 MPa e um Coeficiente de Segurança
valendo 4,
𝜎𝑒𝑠𝑐𝑜𝑎𝑚𝑒𝑛𝑡𝑜 350
𝜎(𝑎𝑑𝑚𝑖𝑠𝑠í𝑣𝑒𝑙) = = = 87,5𝑀𝑃𝑎
𝐶. 𝑆. 4

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

Gráfico 1: Variação do tamanho do elo terra (elo 1) em função do ângulo do elo 2.

Gráfico 2: Variação da velocidade do elo terra (elo 1) em função do ângulo do elo 2.


Gráfico 3: Variação da aceleração do elo terra (elo 1) em função do ângulo do elo 2.

Gráfico 4: Variação do torque do elo 2 em função do ângulo do elo 2.


Gráfico 5: Variação da velocidade do elo 2 em função do ângulo do elo 2.

Gráfico 6: Variação da aceleração do elo 2 em função do ângulo do elo 2.


8. Código do Programa

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: ");

//Calculo dos tamanhos de elos, velocidade e força iniciais


L2= (1200-mat)/3;
L3=2*L2;
printf("L2 será: %d [mm]\n", L2);
printf("L3 será: %d [mm]\n", L3);
W2=(Mat1+Mat2+Mat3+18);
printf("W2 será: %d [rad/s]\n", W2);
Fp=100*W2;
printf("Fp será: %d [N]\n\n", Fp);

//Definindo angulos iniciais (em radianos)


Teta1=3*%pi/2;
Teta2=0:(1*%pi/180):(359*%pi/180);

//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

////Definição das acelerações dos centros de massa


Ag4y=A1;
Lga=L3/2;
Ag3x=((-L2*(W2^2)*cos(Teta2))-(Lga*Alfa3*sin(Teta3))-(Lga*(W3^2)*cos(Teta3)));
Ag3y=((-L2*(W2^2)*sin(Teta2))+(Lga*Alfa3*cos(Teta3))-(Lga*(W3^2)*sin(Teta3)));
Lg2=L2/2;
Ag2x=(-Lg2*(W2^2)*cos(Teta2+%pi));
Ag2y=(-Lg2*(W2^2)*sin(Teta2+%pi));
Ag4x=0;
Ag4y=-A1;

////Definição dos vetores posição das forças


////Elo 3 (Biela)
R23x=(Lga*cos(Teta3+%pi));
R23y=(Lga*sin(Teta3+%pi));
R34x=(Lga*cos(Teta3));
R34y=(Lga*sin(Teta3));
////Elo 2 (Manivela)
R32x=(Lg2*cos(Teta2+%pi));
R32y=(Lg2*sin(Teta2+%pi));
R12x=(Lg2*cos(Teta2));
R12y=(Lg2*sin(Teta2));

////Definição da Massa e Momentos de Inércia


////Elo 4 (Bloco)
printf("Insira os dados do material do Bloco (Elo 4)\n");
Ro4=input("Digite a massa específica [kg/m^3] do elo 4: ");
L4=input("Digite o comprimento do elo 4 em mm: ");
Larg4=input("Digite a largura do elo 4 em mm: ");
Altu4=input("Digite a altura do elo 4 em mm: ");
Massa4=(Ro4*(L4/1000)*(Larg4/1000)*(Altu4/1000));
printf("Massa do elo 4 [kg]:");
disp(Massa4);
////Elo 3 (Biela)
printf("\nInsira os dados do material da Biela (Elo 3)\n");
Ro3=input("Digite a massa específica [kg/m^3] do elo 3: ");
Larg3=input("Digite a largura do elo 3 em mm: ");
Altu3=input("Digite a altura do elo 3 em mm: ");
Massa3=(Ro3*(L3/1000)*(Larg3/1000)*(Altu3/1000));
J3=((((Larg3/1000)^2)+((L3/1000)^2))*Massa3/12);
printf("Massa do elo 3 [kg]:");
disp(Massa3);
printf("Momento de inércia do elo 3 [m^4]:");
disp(J3);
////Elo 2 (Manivela)
printf("\nInsira os dados do material da Manivela (Elo 2)\n");
Ro2=input("Digite a massa específica [kg/m^3] do elo 2: ");
Larg2=input("Digite a largura do elo 2 em mm: ");
Altu2=input("Digite a altura do elo 2 em mm: ");
Massa2=(Ro2*(L2/1000)*(Larg2/1000)*(Altu2/1000));
J2=((((Larg2/1000)^2)+((L2/1000)^2))*Massa2/12);
printf("Massa do elo 2 [kg]:")
disp(Massa2);
printf("Momento de inércia do elo 3 [m^4]:");
disp(J2);

////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);

//Plotagem dos Gráficos

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

CT BORRACHA. Propriedades das borrachas vulcanizadas. Disponível em:


<https://www.ctborracha.com/borracha-sintese-historica/propriedades-das-borrachas-
vulcanizadas/propriedades-tribologicas/>. Acesso em: 20 dez. 2018.

LUZ, Gelson. Aço SAE 10201; propriedades mecânicas. Disponível em:


<https://www.materiais.gelsonluz.com/2017/10/aco-sae-1020-propriedades-
mecanicas.html >. Acesso em: 20 dez. 2018.

Você também pode gostar