Plataforma Automatizada de Misiles
Plataforma Automatizada de Misiles
REALIZADO POR:
SANGOLQUÍ, 2012-09-19
CERTIFICACIÓN DE LA ELABORACIÓN DEL
PROYECTO
_______________________ ________________________
DIRECTORA CODIRECTOR
Sangolquí, 2012-09-06
ii
LEGALIZACIÓN DEL PROYECTO
MILITAR SEMEJANTE.”,
ELABORADO POR:
_______________________ _______________________
_______________________
Sangolquí, 2012-09-19
iii
CARRERA DE INGENIERÍA MECATRÓNICA
AUTORIZACIÓN
_______________________ _______________________
Sangolquí, 2012-09-19
iv
DEDICATORIA
v
DEDICATORIA
vi
AGRADECIMIENTOS
vii
AGRADECIMIENTOS
Agradezco a mi Dios por todo su amor, por darme la sabiduría y por llenarme de
su gracia para poder culminar una etapa más de mi vida, le agradezco por que
ha sido bueno y en su voluntad esta puesta mi esperanza.
Agradezco a mis papis Freddy y Sandra porque he visto en ellos todo el amor y
el esfuerzo diario para darnos lo mejor a mi y a mis hermanos y por darme toda
la confianza.
Agradezco a todos los ingenieros del CICTE por todo su apoyo y amistad.
viii
ÍNDICE DE CONTENIDOS
ÍNDICE DE TABLAS...........................................................................................xx
ix
1.11.1. RECONOCIMIENTO DE LA NECESIDAD. ........................................ 6
1.11.2. DEFINICIÓN DEL PROBLEMA. ......................................................... 6
1.11.3. SÍNTESIS. .......................................................................................... 7
1.11.4. ANÁLISIS Y OPTIMIZACIÓN. ............................................................ 7
1.11.5. EVALUACIÓN..................................................................................... 7
1.11.6. PRESENTACIÓN. .............................................................................. 7
CAPITULO 2: DISEÑO DE LA PLATAFORMA ................................................. 8
2.1 CONSIDERACIONES GENERALES DE DISEÑO. ................................... 8
2.1.1 GENERALIDADES DEL MISIL IGLA. .................................................... 8
2.1.2 ESPECIFICACIONES TÉCNICAS DEL MISIL IGLA. ............................ 8
2.1.3 PARTES DEL MISIL IGLA. .................................................................... 9
2.1.4 PROCEDIMIENTO DE DISPARO DEL MISIL IGLA. ........................... 12
2.2 DESCRIPCIÓN DE PLATAFORMAS DE LANZAMIENTO DE MISILES
IGLA SEMEJANTES. ........................................................................................ 13
2.2.1 SISTEMA ANTIAÉREO GHIBKA ......................................................... 13
2.2.2 SISTEMA ANTIAÉREO ZIPKIN KMS .................................................. 14
2.2.3 SISTEMA DE LANZAMIENTO DE MISILES STRELA SA-7A 9K32-2 . 14
2.3 ANÁLISIS DE LOS SISTEMAS DE MOVIMIENTO DE EQUIPOS
MILITARES AFINES. ........................................................................................ 15
2.3.1 ANÁLISIS DE LOS SISTEMAS DE MOVIMIENTO DE LA
PLATAFORMA PARA LA AMETRALLADORA MÚLTIPLE CALIBRE 0.50. ..... 15
2.3.2 ANÁLISIS DE LOS SISTEMAS DE MOVIMIENTO DEL SISTEMA DE
PUNTERÍA HMOS UTILIZADO EN HELICÓPTEROS...................................... 20
2.3.3 SELECCIÓN DE LOS SISTEMAS MÁS ADECUADOS PARA LA
PLATAFORMA. ................................................................................................. 23
2.4 ESTUDIO DE MOVIMIENTO DE LA PLATAFORMA. ............................. 35
2.4.1 CINEMÁTICA DE LA PLATAFORMA. ................................................. 35
2.4.2 SIMULACIÓN DE MOVIMIENTO EN MATLAB. .................................. 44
2.4.3 DINÁMICA DEL SISTEMA DE ACOPLAMIENTO DE LA
PLATAFORMA. ...................................................................................... 48
2.5 SELECCIÓN DE ACTUADORES. ........................................................... 64
2.6 SELECCIÓN DE LAS CAJAS REDUCTORAS. ...................................... 70
x
2.7 DISEÑO DEL SISTEMA DE ORIENTACIÓN DE LA PLATAFORMA...... 73
2.7.1. ESQUEMA. .......................................................................................... 73
2.7.2. ANÁLISIS DE CARGAS ESTÁTICAS. ................................................. 74
2.7.3. DISEÑO DE ENGRANAJE CORONA-SATÉLITES. ............................ 93
2.7.4. DISEÑO DEL AFUSTE INFERIOR. ................................................... 108
2.7.5. DISEÑO DE LA BASE. ...................................................................... 115
2.7.6. DISEÑO DEL AFUSTE SUPERIOR. ................................................. 119
2.7.7. DISEÑO DE EJES. ............................................................................ 127
2.7.8. SELECCIÓN DE RODAMIENTOS. .................................................... 134
2.7.9. SELECCIÓN DE ACOPLES. ............................................................. 137
2.7.10. ANÁLISIS DE LOS DIÁMETROS REQUERIDOS PARA LOS
PERNOS……...………………………………………………………….139
2.8 DISEÑO DEL SISTEMA DE ELEVACIÓN DE LA PLATAFORMA …….147
2.8.1 ESQUEMA. ........................................................................................ 147
2.8.2 ANÁLISIS DE CARGAS ESTÁTICAS. ............................................... 148
2.8.3 DISEÑO DEL BRAZO DE LA PLATAFORMA. .................................. 169
2.8.4 DISEÑO DE EJES Y BUJES. ............................................................ 170
2.8.5 DISEÑO DE ENGRANAJES HELICOIDALES. .................................. 179
2.8.6 DISEÑO DE LA BANCADA................................................................ 182
2.8.7 SELECCIÓN DE RODAMIENTOS. .................................................... 185
2.8.8 SELECCIÓN DE ACOPLES. ............................................................. 187
2.9 SIMULACIÓN DE ESFUERZOS EN LA PLATAFORMA ....................... 188
2.10 PLANOS CONSTRUCTIVOS DE LA PLATAFORMA. .......................... 191
CAPITULO 3: AUTOMATIZACIÓN DE LA PLATAFORMA .......................... 198
3.1 DESCRIPCIÓN DEL SISTEMA DE CONTROL DE LA PLATAFORMA.198
3.2 DESCRIPCIÓN DE LOS REQUERIMIENTOS DEL SISTEMA. ............ 198
3.3 DIAGRAMA DE BLOQUES DE LA PLATAFORMA. ............................. 199
3.4 ANÁLISIS DEL TIEMPO DE RESPUESTA MÍNIMO............................. 201
3.5 SELECCIÓN Y DESCRIPCIÓN DE LOS SERVO DRIVES. ................. 202
3.6 REQUERIMIENTOS DEL SISTEMA DE COMUNICACIONES. ............ 206
3.7 SELECCIÓN DE LOS COMPONENTES DEL SISTEMA
DE CONTROL ....................................................................................... 210
xi
3.7.1 SELECCIÓN DE SENSORES. .......................................................... 210
3.7.2 SELECCIÓN DEL CONTROLADOR. ................................................ 215
3.7.3 SELECCIÓN DE LOS PERIFÉRICOS. .............................................. 228
3.7.4 SELECCIÓN DE LAS FUENTES DE ALIMENTACIÓN. .................... 231
3.8 DIAGRAMAS DE CONEXIONES ELÉCTRICAS................................... 234
3.9 DESARROLLO DEL SOFTWARE DE LA MÁQUINA............................ 234
3.9.1 CONFIGURACIÓN INICIAL DE LOS SERVO DRIVES. .................... 234
3.9.2 ALGORITMO DE PROGRAMACIÓN. ................................................ 236
3.10 DESARROLLO DE LA INTERFAZ HUMANO-MÁQUINA (HMI). .......... 243
3.11 SIMULACIÓN DEL SISTEMA DE CONTROL SIN CARGA. ................. 251
3.11.1 SIMULACIÓN DEL ALGORITMO DE SEGUIMIENTO. ..................... 251
CAPITULO 4: PRUEBAS DE FUNCIONAMIENTO ....................................... 254
4.1. PRUEBAS DE COMUNICACIÓN. ......................................................... 254
4.1.1. PRUEBAS DE COMUNICACIÓN COMPUTADOR - PLC.................. 254
4.1.2. PRUEBAS DE COMUNICACIÓN COMPUTADOR - SERVO
DRIVE................................................................................................... 258
4.1.3. PRUEBAS DE COMUNICACIÓN COMPUTADOR – PANTALLA
TÁCTIL. ................................................................................................ 260
4.1.4. PRUEBAS DE COMUNICACIÓN COMPUTADOR - ENCODER. ...... 262
4.1.5. PRUEBAS DE COMUNICACIÓN SERVO DRIVE - PLC. .................. 264
4.1.6. PRUEBAS DE COMUNICACIÓN PANTALLA TÁCTIL - PLC. ........... 265
4.1.7. PRUEBAS DE COMUNICACIÓN ENCODER - PLC. ........................ 266
4.2. PRUEBAS DE FUNCIONAMIENTO SIN CARGA. ................................ 266
4.3. DESCRIPCIÓN DE LA BASE DE FUEGO ZU-23mm. .......................... 268
4.4. PRUEBAS DE MOVIMIENTO DE LA PLATAFORMA........................... 268
4.4.1. DISEÑO DEL ACOPLAMIENTO PROVISIONAL............................... 269
4.4.2. PRUEBA DE MOVIMIENTO DESDE JOG......................................... 271
4.4.3. PRUEBA DE MOVIMIENTO DESDE EL JOYSTICK. ........................ 272
4.4.4. PRUEBA DE MOVIMIENTO EN MODO AUTOMÁTICO. .................. 273
4.4.5. PRUEBA DE MOVIMIENTO EN MODO POSICIÓN. ......................... 273
4.4.6. PRUEBAS MOVIMIENTO PARA EVALUACIÓN DE
PARÁMETROS DE OPERACIÓN. ....................................................... 273
xii
CAPITULO 5: CONCLUSIONES Y RECOMENDACIONES .......................... 293
5.1 CONCLUSIONES.................................................................................. 293
5.2 RECOMENDACIONES. ........................................................................ 295
xiii
ÍNDICE DE FIGURAS
xiv
xv
xvi
xvii
xviii
xix
ÍNDICE DE TABLAS
xx
ÍNDICE DE APÉNDICES
xxi
RESUMEN
El diseño parte de las características técnicas del misil Igla con el que cuenta la
fuerza terrestre, para lo cual se describen sus partes y procedimientos de
disparo. Tales procedimientos son materia importante para establecer la
robustez y el tiempo de respuesta del sistema.
xxii
CAPÍTULO 1
INTRODUCCIÓN
1
• Modernizar y repotenciar equipos militares abarcando estudios en las áreas
antes mencionadas.
Debido a las dificultades planteadas por los convenios internacionales así como
por el alto costo, la adquisición de los sistemas automáticos de lanzamiento de
misiles antiaéreos no ha sido viable para nuestro país.
2
1.6. JUSTIFICACIÓN DEL PROYECTO.
3
1.8. MARCO INSTITUCIONAL.
1.9.1. GENERAL.
1.9.2. ESPECÍFICOS.
4
• Integrar los componentes necesarios para la automatización en una base de
fuego antiaéreo.
• Realizar pruebas de funcionamiento de la base de fuego antiaéreo
automatizada.
1
SHYGLEY Joseph, “Diseño en Ingeniería Mecánica”, P. 4
5
Figura 1.1. Fases del Diseño2.
2
Ibíd., P. 4
6
1.11.3. SÍNTESIS.
1.11.5. EVALUACIÓN.
1.11.6. PRESENTACIÓN.
7
CAPITULO 2
DISEÑO DE LA PLATAFORMA
El misil Igla es de tipo tierra – aire (SAM, surface air misil) de fabricación rusa
desarrollado en 1971, con un sistema de seguimiento automático a través de
cámaras infrarrojas. Está diseñado para proporcionar protección a tropas de
infantería, blindados, instalaciones militares y zonas estratégicas, frente a
aeronaves que vuelan a baja altura en cursos de encuentro o persecución.
Fabricación: Rusa.
Disparo: Manual.
8
Distancia mínima de disparo: 1000 [m].
[Link] Tubo.
9
• Aparatos de puntería, miras alineadas para compensar la caída del misil en
su salida del tubo.
• Mecanismo de encendido, palanca que conecta la fuente de alimentación.
• Enchufe del mecanismo de disparo, terminales de conexión entre el misil y
el disparador.
• Iniciador, enciende motores de expulsión y de marcha.
• Abrazaderas, para el transporte manual del misil.
• Tapas, protegen al misil de contacto con polvo y humedad.
[Link] Misil.
10
[Link] Mecanismo de disparo.
11
La fuente de poder está constituida por las siguientes partes:
1. Sacar las tapas anterior y posterior del tubo del misil en la posición de
cuclillas.
2. Colocar el misil al hombro.
3. Colocar la correa al costado derecho para no dificultar la visibilidad.
4. Desplazar hacia afuera los aparatos de puntería.
5. Esperar visualización con la aeronave.
6. Una vez que el blanco se encuentra visible, activa la fuente de
alimentación con el mecanismo del tubo.
7. Se apunta al blanco y se presiona el rabillo disparador en la posición
media.
8. Los indicadores (chicharra y lámpara) serán intermitentes hasta
enganchar al blanco, en ese momento su señal será continua. Este
procedimiento de disparo se conoce como manual.
9. Una vez enganchado el blanco se pasa el rabillo disparador a la posición
final, de disparo.
10. El misil saldrá del tubo y activará el motor de marcha para dar
seguimiento o encuentro a la aeronave a través de la cabeza de
detección térmica.
12
Se puede presionar el rabillo disparador hacia la posición de disparo
inmediatamente después de apuntar a la aeronave, este procedimiento de
emergencia se conoce como disparo automático. La efectividad de este último
es menor a la del disparo manual.
El sistema Ghibka (ver figura 2.6) es un sistema de defensa antiaéreo ruso para
el lanzamiento de misiles Igla desarrollado por el Instituto Radioelectrónico
Altair y la fábrica radiotécnica Ratep; ambas empresas forman parte del Grupo
Almaz-Antey, el consorcio Almaz-Antey es la mayor empresa del sector militar-
industrial de Rusia en la fabricación de armamento de defensa antiaérea,
defensa antimisil y equipos electrónicos de guerra de aplicación diversa.
13
El sistema aprovecha los módulos de lanzamiento tipo Strelets, para seis
misiles antiaéreos guiados Igla, desarrollados por la oficina de diseño de
Kolomna (en la provincia de Moscú). El peso total del sistema es de unas 1,3
toneladas.
El sistema antiaéreo de corto alcance Zipkin KMS (ver figura 2.7), es un sistema
desarrollado en Turquía, capaz de lanzar seis misiles Igla, posee un sistema
óptico-electrónico que permite visualizar a la aeronave y dar seguimiento
continuo de la misma. Está diseñado para trabajar conjuntamente con otros
sistemas de defensa, y puede ser instalado sobre tanques, camiones y buques
de guerra.
14
el hombro, o desde la plataforma automatizada (ver figura 2.8). El SA-7 es la
primera generación de manpads soviéticos; es un misil distraído fácilmente por
el calor solar. Este misil debe ser disparado a la cola del avión para que pueda
ser efectivo, es guiado por un sistema óptico-infrarrojo. Su sistema de
lanzamiento automático es muy semejante a las plataformas de lanzamiento del
Igla, este sistema además posee un radar instalado, a diferencia de los
sistemas lanzadores de misil Igla que poseen el radar externo.
15
El sistema consta de una plataforma dirigida por un operario a través de un
joystick mecánico que permite movimientos de orientación y elevación de las
ametralladoras. Este sistema es capaz de disparar a blancos a una distancia de
1 [km] con gran efectividad y precisión. (Ver figura 2.9).
16
Este movimiento lento lo logra gracias a un sistema de transmisión piñón-
corona (ver figura 2.10), en donde la corona de gran tamaño permite girar en su
interior un piñón pequeño.
17
mecanismos de 4 barras para la conexión y desconexión de las mismas (ver
figura 2.11), las mismas que son controladas por el joystick principal.
19
2.3.2 ANÁLISIS DE LOS SISTEMAS DE MOVIMIENTO DEL SISTEMA
DE PUNTERÍA HMOS UTILIZADO EN HELICÓPTEROS.
20
Figura 2.14. Sistema de transmisión en elevación del HMOS.
21
Figura 2.15. Lado derecho del HMOS.
22
Figura 2.16. Vista frontal del HMOS.
23
orientación debe ser preciso y debe proporcionar una base robusta para los
sistemas que se localicen en la parte superior.
24
3. Utilizar un sistema de engranaje epicicloidal corona-satélites, con porta-
planetas fijo. Con tres piñones en el interior de la corona, uno de estos
transmite el torque del servomotor, hace girar a la corona alrededor de los
planetas y los otros dos trabajan como engranes locos para mantener la
estabilidad del sistema. (Ver figura 2.19).
25
Tabla 2.1. Ventajas y desventajas de los sistemas de engranaje epicicloidal.
26
En base a la tabla comparativa anterior y teniendo como referencia el sistema
de orientación de la plataforma de la ametralladora 0.50, se ha considerado que
la mejor opción es el sistema de engranaje epicicloidal corona-satélites con
corona fija, debido a que permite tener un mayor espacio para manipular los
elementos y dispositivos internos ya que no requiere de eje central, posee la
corona fija y esto permite mayor estabilidad en el sistema, además de que su
acoplamiento con la estructura superior es mucho más sencillo con relación a
las otras opciones posibles, se acopla mediante el afuste superior y gracias a la
presencia del rodamiento la corona no se ve forzada a soportar todo el peso de
la estructura superior.
27
Figura 2.22. Afuste inferior.
28
Figura 2.23. Sistema de rodamiento.
29
Otro factor importante considerado para la selección de esta relación de
transmisión, es que la misma nos permite tener cinco dientes en contacto en
todo momento entre la corona y cada uno de los piñones, esta ventaja permite
tener mayor precisión a la hora de rotar lentamente la plataforma en busca de
un objetivo.
Los tres piñones se encuentran sujetos al afuste superior (ver figura 2.25), uno
de ellos se conecta con el eje del motor, mientras los otros dos serán engranes
locos que ayuden a transmitir el movimiento y a dar estabilidad a todo el
sistema. Para acoplar los piñones con el afuste superior, se utilizará tres ejes,
los mismos que para el caso de los engranes locos irán en un extremo unidos a
los engranes mediante un acople, y en su otro extremo se ubicarán sobre un
par de rodamiento ubicados dentro del afuste superior, esto debido a las
condiciones del ensamblaje de la plataforma. Para el caso del engrane que
transmite el torque, su eje lo acoplará al mismo con el eje del motor mediante
un acople sin chavetas (ver figura 2.26).
30
Figura 2.26. Acople sin chavetas y rodamientos de los engranes locos.
31
Figura 2.28. Tapa superior del afuste superior.
∙ 2.1
32
En base al estudio realizado en sistemas afines se plantean dos alternativas
para el movimiento de elevación:
33
Los engranajes más adecuados debido a la precisión y velocidades de
funcionamiento son los helicoidales, ya que el movimiento es suave debido al
ángulo de hélice.
34
2.4 ESTUDIO DE MOVIMIENTO DE LA PLATAFORMA.
35
!
tan 2.4]
$
% %
tan " = = [2.5]
!" + $"
!
tan, - .
$
%
" tan, / 0
!" # $"
!1 ∙ $ − ! ∙ $1
1 [2.6]
!" + $"
% 1 ∙ !" + $" ! ∙ !1 + $ ∙ $1 ∙ %
"1 = − [2.7]
!^2 + $^2 + %^2 !" + $" + %" ∙ !" + $"
!4 ∙ $ − ! ∙ $4 ∙ ! " + $ " − 2 ∙ !1 ∙ $ − ! ∙ $1 ∙ ! ∙ !1 + $ ∙ $1
4 = [2.8]
! " + $ " "
!" + $"
"4 = %4 ∙ 5 6
!" + $" + %"
! ∙ !1 + $ ∙ $1 2 ∙ ! " + $ " ∙ ! ∙ !1 + $ ∙ $1 + % ∙ %1
+%1 ∙ 5 −
!" + $" + %" ∙ !" + $" ! " + $ " + % " "
! ∙ !1 + $ ∙ $1
− 6
!" + $" + %" ∙ !" + $"
36
! ∙ !4 + $ ∙ $4 + !1 " + $1 " 2 ∙ ! ∙ !1 + $ ∙ $1 ∙ ! ∙ !1 + $ ∙ $1 + % ∙ %1
−% ∙ 7 −
!" + $" + %" ∙ !" + $" ! " + $ " + % " " ∙ ! " + $ "
! ∙ !1 + $ ∙ $1 "
− 89 : [2.9]
!" + $" + %" ∙ !" + $" "
!
= tan, - .
$
%
" = tan, / 0
!" + $"
!1 ∙ $ − ! ∙ $1
1 =
!" + $"
% 1 ∙ !" + $" ! ∙ !1 + $ ∙ $1 ∙ %
"1 = −
!^2 + $^2 + %^2 !" + $" + %" ∙ !" + $"
−2 ∙ !1 ∙ $ − ! ∙ $1 ∙ ! ∙ !1 + $ ∙ $1
4 =
! " + $ " "
! ∙ !1 + $ ∙ $1 2 ∙ ! " + $ " ∙ ! ∙ !1 + $ ∙ $1 + % ∙ %1
"4 = %1 ∙ 5 −
!" + $" + %" ∙ !" + $" ! " + $ " + % " "
! ∙ !1 + $ ∙ $1
− 6
!" + $" + %" ∙ !" + $"
!1 " + $1 " 2 ∙ ! ∙ !1 + $ ∙ $1 ∙ ! ∙ !1 + $ ∙ $1 + % ∙ %1
−% ∙ 7 −
!" + $" + %" ∙ !" + $" ! " + $ " + % " " ∙ ! " + $ "
! ∙ !1 + $ ∙ $1 "
− 89 :
!" + $" + %" ∙ !" + $" "
37
Utilizando Matlab se calculan los resultados para las trayectorias de la
aeronave. El archivo m escrito para llevar a cabo los cálculos y mostrar los
resultados es el siguiente:
function trayectoria(x,y,z,t0,tmax)
clc
m=10000;
%Vector tiempo.
deltha=(tmax-t0)/m;
t =t0:deltha:tmax;
%Trayectoria de la aeronave [m].
x=x(t);
y=y(t);
z=z(t);
if(length(z)==1)
z(1:m+1)=z;
end
if(length(y)==1)
y(1:m+1)=y;
end
if(length(z)==1)
z(1:m+1)=z;
end
%Velocidad de la aeronave [ m / s ].
vx = diff(x)/deltha;
vy = diff(y)/deltha;
vz = diff(z)/deltha;
%Aceleración de la aeronave [ m / s ^ 2 ].
ax = diff(vx)/deltha;
ay = diff(vy)/deltha;
az = diff(vz)/deltha;
%Igualar dimensiones.
38
t=t(1:m-1);
x=x(1:m-1);
y=y(1:m-1);
z=z(1:m-1);
vx=vx(1:m-1);
vy=vy(1:m-1);
vz=vz(1:m-1);
%Rango de la gráfica.
% Create the plot using the parametric functions
figure(1);set(gcf,'Name','Trayectoria de la Aeronave’,
‘NumberTitle’,’off’,’DefaultLineLineWidth',3);
plot3(x,y,z,'LineWidth',2,'Color','b');
hold on
plot3(0,0,0,'go','MarkerFaceColor',[0.2 0.5 0.3],'MarkerSize',15)
hold off
xlabel('X [m]');
ylabel('Y [m]');
zlabel('Z [m]');
grid(gca,'minor')
grid on
%%Cálculo de movimiento de la plataforma para las trayectorias dadas.
theta1 = ( atan ( x ./ y ) ) * 180 / pi;
theta2 = atan ( z ./ sqrt ( x .^ 2 + y .^ 2 ) ) * 180 / pi;
ohmega1 = ( ( vx .* y - x .* vy ) ./ ( x .^ 2 + y .^ 2) ) * 60 / ( 2 * pi );
ohmega2 = ( ( ( vz .* ( x .^ 2 + y .^ 2 ) ) - ( ( x .* vx + y .* vy ) .* z ) ) ./ ( ( x .^ 2
+ y .^ 2 + z .^ 2 ) .* sqrt ( x .^ 2 + y .^ 2 ) ) ) * 60 / ( 2 * pi );
alpha1 = ( ( ( ax .* y - x .* ay ) .* ( x .^ 2 + y .^ 2 ) - 2 * ( vx .* y - x .* vy ) .* ( x
.* vx + y .* vy ) ) ./ ( ( x .^ 2 + y .^ 2 ) .^ 2 ) );
x1 = ( x .^ 2 + y .^ 2 + z .^ 2 );
x2 = ( x .^ 2 + y .^ 2 );
x3 = ( x .* vx + y .* vy + z .* vz );
39
x4 = ( x .* vx + y .* vy );
x5 = ( x .* ax + y .* ay + vx .^ 2 + vy .^ 2 );
alpha2 = az .* ( x1 .^ ( - 1 ) .* x2 .^ ( - 0.5 ) ) + vz .* ( x4 .* x1 .^ ( - 1 ) .* x2 .^ (
- 0.5 ) - 2 * x3 .* x2 .^ ( 0.5 ) .* x1 .^ ( - 2 ) - x4 .* x1 .^ ( - 1 ) .* x2 .^ ( - 0.5 ) ) - z
.* ( x5 .* x1 .^ ( - 1 ) .* x2 .^ ( - 0.5 ) - 2 * x4 .* x3 .* x1 .^ ( - 2 ) .* x2 .^ ( - 0.5 ) -
x4 .^ ( 2 ) .* x1 .^ ( - 1 ).* x2 .^ ( - 3 / 2 ) );
%Gráfica de la solución.
figure(2);subplot(3,2,1);plot(t(1:m-1),theta1(1:m-1));
xlabel('Time (s)');
ylabel('Theta 1 (grado)');
grid(gca,'minor')
grid on
figure(2);subplot(3,2,3);plot(t,ohmega1);
xlabel(' Time (s) ');
ylabel(' Ohmega 1 (rpm) ');
grid(gca,'minor')
grid on
figure(2);subplot(3,2,5);plot(t,alpha1);
xlabel('Time (s) ');
ylabel('Alpha 1 (rad/s^2)');
grid(gca,'minor')
grid on
figure(2);subplot(3,2,2);plot(t,theta2);
xlabel(' Time (s) ');
ylabel(' Theta 2 (grado) ');
grid(gca,'minor')
grid on
figure(2);subplot(3,2,4);plot(t,ohmega2);
xlabel(' Time (s) ');
ylabel(' Ohmega 2 (rpm) ');
grid(gca,'minor')
40
grid on
figure(2);subplot(3,2,6);plot(t,alpha2);
xlabel(' Time (s) ');
ylabel(' Alpha 2 (rad/s^2) ');
grid(gca,'minor')
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Cinemática de la Plataforma');
set(gca,'Visible','off');
set(h,'Visible','on');
[max(ohmega1) max(ohmega2);max(alpha1) max(alpha2)]
41
y=@(t)1000;
z=@(t)1000;
trayectoria(x,y,z,0,10000/360)
a=input('Siguiente Trayectoria')
%Trayectoria 2.
x=@(t)110*sin(3*t);
y=@(t)(100*cos(4*t)+200);
z=@(t)150*sin(2*t)+200;
trayectoria(x,y,z,-0.5,0.5)
Parámetros Valores
Fuente: Propia
43
Figura 2.36. Movimientos (ɵ1, ɵ2) para la segunda trayectoria de la
aeronave.
Parámetros Valores
Fuente: Propia
44
cinemático se considera además un tercer grado de libertad prismático
correspondiente a la distancia entre el foco de las cámaras de seguimiento y la
aeronave.
B
; A 2 0 ℎE
5;" 6 @@B D
0 " 0D
;8 @2 D
?0 0 0C
Donde:
Matlab cuenta con una herramienta de robótica, ventajosa para realizar estudios
de movimiento de robots de varios grados de libertad. Sobre estas coordenadas
podemos aplicar el siguiente programa para obtener el movimiento del robot
siguiendo una trayectoria.
45
%%Definición del robot.
L1=link([pi/2 0 0 700]);
L2=link([pi/2 0 0 0]);
L3=link([0 0 0 0 1]);
r=robot({L1,L2,L3});
%%Cálculo de la cinemática inversa.
tic
q = ikine(r, T);
toc;
qd=diff(q)./delthat;
qd=[qd(1,:); qd];
qdd=diff(qd)./delthat;
qdd=[qdd(1,:); qdd];
%%Gráficas de resultados.
theta1=q(:,1)*180/pi+90;
res=mod(max(abs(theta1)),360);
n=(max(abs(theta1))-res)/360;
theta1=-(abs(theta1)-n*360);
figure(1);subplot(3,2,1);plot(t,theta1);
xlabel('Time (s)');
ylabel('Theta 1 (grado)')
grid on
figure(1);subplot(3,2,3);plot(t,-qd(:,1)*60/(2*pi));
xlabel('Time (s)');
ylabel('Ohmega 1 (rpm)')
grid on
figure(1);subplot(3,2,5);plot(t,-qdd(:,1));
xlabel('Time (s)');
ylabel('Alpha 1 (rad/s^2)')
grid on
theta2=q(:,2)*180/pi;
46
res=mod(max(abs(theta2)),360);
n=(max(abs(theta2))-res)/360;
theta2=90-(abs(theta2)-n*360);
figure(1);subplot(3,2,2);plot(t,theta2);
xlabel('Time (s)');
ylabel('Theta 2 (grado)')
grid on
figure(1);subplot(3,2,4);plot(t,qd(:,2)*60/(2*pi));
xlabel('Time (s)');
ylabel('Ohmega 2 (rpm)')
grid on
figure(1);subplot(3,2,6);plot(t,qdd(:,2));
xlabel('Time (s)');
ylabel('Alpha 2 (rad/s^2)')
grid on
figure(2);plot(r,q);
xlim([-5000000 5000000]);
ylim([-1500000 0]);
zlim([0 1500000]);
47
Figura 2.38. Resultados de la simulación (ɵ1, ɵ2) en Matlab.
48
Figura 2.39. Esquema general del acoplamiento de actuadores3.
GH
F 2.10
GI
KLMN OLMN ∙ GI
J 2.11
K O ∙ GH
OLMN − OI PI ∙ <I
49
PI ∙ <I # OI
OH # PH ∙ <H
F∙J
PI ∙ <I # OI
OH # PH ∙ F ∙ <I
F∙J
PI OI
OH - # F ∙ PH . ∙ <I #
F∙J F∙J
PI OI
OH - # F ∙ PH . ∙ <I + [2.14]
F F
>OH
=0 [2.15]
>F
PI OI
-− + PH . ∙ <I − " = 0
F " F
1 OI
FQRS = T ∙ -PI + . [2.16]
PH <I
PI 1 OI OI
OH = V +T ∙ -PI + . ∙ PH ` ∙ <I +
PH <I
WX ∙ ZPI + ]\ ^ WX ∙ ZPI + ]\ ^
[ [
U Y \ _ Y \
50
Más adelante se verá la utilidad de despejar la relación del torque del motor
sobre la raíz de su inercia, separando los miembros de la ecuación de acuerdo
a la procedencia de sus términos.
OH
2 ∙ PI ∙ <I " + OI ∙ <I
PH
KH = OH ∙ GH = OH ∙ FQRS ∙ GI [2.17]
KH = KI [2.18]
PI OI
KI = a- + F ∙ PH . ∙ <I + b ∙ F ∙ GI
F F
KI = [ PI + F " ∙ PH ∙ <I + OI ] ∙ GI
1 OI
KH = a-PI + ∙ -PI + . ∙ PH . ∙ <I + OI b ∙ GI
PH <I
KH = 2 ∙ PI ∙ <I + OI ∙ GI
51
Figura 2.40. Ley de movimiento del servomotor. Triangular o trapezoidal.
ef ∙
c= = [2.19]
d − d d
e ∙
∝= = " [2.20]
d 2 d ∙ d d
Donde:
d
ef = [2.21]
d 2 d
d"
e = [2.22]
d 2 d ∙ d
ef
c= ∙<∙d
e
e ∙ c
=d
ef ∙ <
52
Introduciéndola en la ecuación de la potencia podemos obtener su derivada
respecto al tiempo.
ef
KI = 2 ∙ ∙ P ∙ < " + O I ∙ <I ∙ d
e I I
ef
K1 I = 2 ∙ ∙ P ∙ < " + O I ∙ <I [2.23]
e I I
ef
Para el caso de movimiento trapezoidal con 9e = 1910, la potencia transitoria
resulta:
2
K1 I = ∙ P ∙ < " + O I ∙ <I
10 I I
ef
Para el caso de movimiento triangular con 9e = 192, la potencia transitoria
resulta:
OH
= 2 ∙ PI ∙ <I " + OI ∙ <I [2.24]
PH
53
OH "
= 4 ∙ PI ∙ <I " + OI ∙ <I = 4 ∙ K1 I
PH
1
hI = ∙ P ∙ GI " [2.25]
2 I
1
hH = ∙ P ∙ F ∙ GI " [2.26]
2 H
1
hH = ∙P ∙ F ∙ GI "
2 H QRS
1 OI
hH = ∙ -PI + . ∙ GI "
2 <I
54
cercana y disponible en el mercado, para evitar la manufactura de una caja de
transmisión personalizada para la aplicación.
F
l= [2.27]
FQRS
"
OH " 1
= PI ∙ <I + OI ∙ <I ∙ - + l.
"
PH l
1
K1 I = ∙ l " + 1 ∙ PI ∙ <I " + OI ∙ <I
2
"
OH " 1
= PI ∙ <I + OI ∙ <I ∙ - + l. = m ∙ K1 I
"
PH l
1
m = 2∙- + 1.
l"
OH " 1
= 2 ∙ - " + 1. ∙ K1 I
PH l
1 " OI
hH = ∙ l ∙ -PI + . ∙ GI "
2 <I
55
Tabla 2.4. Parámetros de servomotores SureServo.
56
thauM=[320 700 1350 2500 3500 4800 9400 14500]/1000;
%Velocidad máxima del servomotor [rad/s].
ohmegaM=2*pi/60*[3000 3000 3000 3000 3000 2000 2000 2000];
57
WM2(i,:)=thauM(i)^2/JM(i);
end
58
WM1k(i,1:1001)=thauM(i).^2./JM(i);
WM2k(i,1:1001)=thauM(i).^2./JM(i);
end
59
hold on
plot(k,log(WM1k),'-b','LineWidth',2);
hold off
xlabel('k');
ylabel('W [W/s]');
title('Servomotores Orientación');
grid on
figure(2);subplot(1,2,2);
plot(k,log(WL2k),'-r','LineWidth',2);
hold on
plot(k,log(WM2k),'-b','LineWidth',2);
hold off
xlabel('k');
ylabel('W [W/s]');
title('Servomotores Elevación');
grid on
%%Cálculo de las relaciones de transmisión mínimas.
k1(1:8)=0;
k2(1:8)=0;
for j=1:8
for i=1:1001
if(abs(log(WL1k(i))-log(WM1(j)))<100*(kmax-kmin)/1000)
k1(j)=k(i);
end
if(abs(log(WL2k(i))-log(WM2(j)))<2500*(kmax-kmin)/1000)
k2(j)=k(i);
end
end
end
k2=k2+0.085;
i1=iOPT1.*k1
60
i2=iOPT2.*k2
thau1=thauL1MIN./i1
thau2=thauL2MIN./i2
ohmega1=(ohmegaL1*0.5*60/pi).*i1
ohmega2=(ohmegaL2*0.5*60/pi).*i2
W1=thau1.*ohmega1*2*pi/60
W2=thau2.*ohmega2*2*pi/60
Parámetros Valores
Fuente: Propia
Fuente: Propia
61
Figura 2.411. Puntos en K-W de los servomotores y la carga.
62
Tabla 2.7. Relaciones de transmisión mínimas de los servomotores.
Fuente: Propia.
Fuente: Propia.
63
Por recomendación de los expertos en desarrollo de proyectos de la institución
auspiciantes se ha tomado un factor de seguridad de 5 para la selección de
actuadores y cajas de transmisión. Los parámetros a emplear en la selección
son los siguientes:
Parámetros Valores
vp [w] 686
vu [w] 1181
Fuente: Propia.
64
En la tabla anterior se analizan las 3 alternativas en base a varios criterios. De
las 3 marcas se prefiere a SureServo gracias a la confiabilidad que ha
manifestado en aplicaciones anteriores realizadas por la institución auspiciante.
Figura 2.43. Diagrama del valor mínimo de grados para el seguimiento del
helicóptero.
65
Se procede a calcular la precisión que debe tener el servomotor, para esto se
realiza un análisis en donde se obtiene el valor mínimo de grados que debe
girar y elevarse la plataforma para poder seguir al helicóptero sin ningún
inconveniente. Como se observa en la figura 2.43.
ℎF
x =
>=!
3.96 []
5000 []
x =
De esta manera observamos que la plataforma debe girar con una precisión de
2’43’’ para tener al helicóptero en la mira.
66
tren de engranes a implementar es de 3 igualmente, por lo tanto podemos
apreciar la precisión a la salida de la caja de transmisión.
| = x ∙ 3 [2.29]
Donde x son los grados que gira la corona y | son los grados que gira el
piñón, en el caso de orientación.
Se tiene:
| = 0.0453° ∙ 3
| = 0.1359° = 8′9′′
Y se obtiene:
360°
}~ =
0.1359°
67
este tipo de comunicación, y los servomotores Sanyodenki no trabajan con
modbus. Por tal motivo y además por no tener alguna recomendación de
profesionales acerca de los mismos se ha descartado la opción de utilizar los
servomotores Sanyodenki.
Ambas marcas son servomotores robustos de uso industrial, pero pueden ser
considerados para trabajar bajo estándar militar, trabajan con comunicación tipo
modbus, y su alimentación varía entre 200 [VAC] y 400 [VAC].
68
Tabla 2.10. Análisis comparativo de servomotores.
Fuente: Propia.
69
Figura 2.44. Servomotor de marca SureServo4.
4
Ibíd., P. e16-33, Vol. 13
70
Tabla 2.11. Torques con los servomotores seleccionados y sus relaciones
de transmisión mínimas.
Fuente: Propia.
Servomotor r [t ∙ q]
SVM-210 10 48.00
SVM-220 25 235.00
Fuente: Propia.
71
Y por el tipo de engranajes se clasifican en:
Las únicas cajas de transmisión que cumplen con este requisito último de ser
autotrabante son las de tornillo sin fin.
72
Fuente: Catálogo “Cajas Reductoras Motovario”, P. 21-22
73
1) Base.
2) Empaque 1.
3) Afuste Inferior.
4) Corona.
5) Engranes rectos.
6) Ejes de los engranes rectos.
7) Rodamientos de los engranes rectos.
8) Rodamiento de bolas grande.
9) Afuste superior parte 1
10) Empaque 2
11) Afuste superior parte 2.
2.7.2. ANÁLISIS DE CARGAS ESTÁTICAS.
74
Figura 2.46. Diagrama del cuerpo libre 1 del sistema de orientación de la
plataforma.
$ = 0 [2.30
# K=2$ # K $ 2.31
En donde:
K $ K ∗ cos 2.33
75
Además se asume que RA es igual a RB ya que la distribución de la masa de la
tapa del afuste superior es simétrica y el centro de gravedad se encuentra en la
mitad de la tapa. De tal manera se tiene:
2 K=2$ # K $
K=2$ # K $
2
Nota: Todos los cálculos que se realizan para el esquema anterior se han
obtenido asumiendo que solo existe un elemento de cada uno de los expuestos
en el sistema ya que de esta manera logramos el caso extremo en donde solo
se tenga por ejemplo un tornillo en lugar de los 16 tornillos y de esta manera en
lugar de dividir la fuerza para los 16 tornillos se asume que la fuerza recae solo
sobre uno.
! 0 2.34
K ! 2.35
K ! K ∗ sin [2.36
76
Análisis de los engranes rectos:
! 0
S S 2.37
Tanto como para los acoples B-109 y B-WK el análisis de fuerzas es el mismo.
(La selección de estos acoples se explica en la sección 2.7.9 del presente
capítulo).
! 0
=e S
77
En donde Nac es la fuerza normal que se genera en el acople debido a la
fuerza NFT.
Figura 2.50. Diagrama del cuerpo libre del eje de los engranes rectos.
! 0
=e 2.39
S
0
Figura 2.51. Diagrama del cuerpo libre en los rodamientos de los engranes
rectos.
78
! 0
" 2.41
" S
Figura 2.52. Diagrama del cuerpo libre de la tapa del afuste superior.
! 0
K=2! # K ! # S 8
R3 es la reacción generada por la presencia del tornillo que sujeta la tapa del
afuste superior con el afuste superior.
79
Análisis en los tornillos que ensamblan la tapa del afuste superior con el
afuste superior:
Figura 2.53. Diagrama del cuerpo libre del tornillo de unión de la tapa del
afuste superior con el afuste superior.
! 0
8
Remplazando R3 se tiene:
K=2! # K ! # S 2.43
Los pesos de los tornillos, rodamientos, ejes, acoples y engranes rectos son
despreciables, debido a que su valor es mínimo en relación al peso de el
soporte superior y su influencia en los cálculos es mínima.
80
El diagrama del cuerpo libre 2 del sistema de orientación muestra la distribución
de fuerzas entre el afuste superior y el rodamiento grande.
$ 0
# # K=1$ #
En donde:
K=1$
#
2
Remplazando RA se tiene:
K=2$ # K $ # K=1$
2.44
2
RC y RD son las reacciones que se generan a lo largo del afuste superior debido
al contacto con el rodamiento grande.
! 0
2 + K=1! 0
81
# K=1!
Remplazando R4 se tiene:
0
$ 0
# # K $ #
En donde:
K $ K ∗ cos
82
Como RC=RD situación explicada anteriormente, y de manera similar se asume
que RE=RF ya que la distribución de la masa en el rodamiento es simétrica y el
centro de gravedad se encuentra en la mitad del mismo. De tal manera se tiene:
2 # K $ 2
K $
# 2.47
2
Remplazando Rc se tiene:
K=2$ # K $ # K=1$ # K $
2
! 0
K ! K ∗ sin
= K ! +
Remplazando R5 se tiene:
83
Figura 2.56. Diagrama del cuerpo libre 3 del sistema de orientación de la
plataforma.
$ 0
# Ke$
En donde:
Ke$ Ke ∗ cos
Entonces se tiene:
84
Ke$ 2
Ke$
2.49
2
! 0
Ke! Ke ∗ sin
$ 0
# # # # K=F$ # X
En donde:
Wafi es el peso del afuste inferior, RE y RF son las reacciones que se generan a
lo largo del afuste inferior debido al contacto con el rodamiento grande, RG y RH
son las reacciones que se generan por el contacto con la corona, y RI y RJ son
las reacciones que se generan en el afuste inferior por el contacto con la base.
85
K=F$
# # 2.51
2
Reemplazando RE y RG se tiene:
! 0
K=F! + = 8
Reemplazando R6 se tiene;
0
Figura 2.57. Diagrama del cuerpo libre del tornillo de unión de la corona y
el afuste inferior.
86
! 0
Figura 2.58. Diagrama del cuerpo libre del tornillo de unión del afuste
inferior con la base.
! 0
" 8
87
El diagrama del cuerpo libre 4 del sistema de orientación muestra la distribución
de fuerzas en la base y en los tornillos que ensamblan el afuste inferior con la
base.
$ 0
# X # K¤$ ¥ # I
En donde:
K¤$ K¤ ∗ cos
K¤$
¥ # 2.55
2
Remplazando RI se tiene:
88
! 0
K¤! K¤ ∗ sin
0
Una vez obtenidas las ecuaciones respectivas para cada elemento que
conforman la plataforma, se procede a ingresar las ecuaciones en el software
Matlab para obtener los valores respectivos de todas las fuerzas, momentos y
reacciones.
89
Wafi= 297.56; %Peso del afuste inferior en (N)
Wb= 268.8; %Peso de la base en (N)
T=124; %Torque generado en la corona en (Nm)
R=0.2583;
d3=0.039; %Distancia entre el centro de gravedad del eje de los engranes y los
puntos de aplicacion de las fuerzas Fen y Feje en (m)
d4=0.03; %Distancia entre el centro de gravedad del eje de los engranes y el
punto de aplicacion de la fuerza Nac en (m)
d5=0.0287; %Distancia entre el centro de gravedad del afsute superior y el
punto de aplicacion de R6 y R5.
d6=0.0307; %Distancia entre el centro de gravedad del afuste inferior y el punto
de aplicacion de R10 y R9.
d7=0.0377; %Distancia entre el centro de gravedad del afuste inferior y el punto
de aplicacion de R8 y R7.
d8=0.0207; %Distancia entre el centro de gravedad del afuste inferior y el punto
de aplicacion de R14.
d9=0.0153; %Distancia entre el centro de gravedad del afuste inferior y el punto
de aplicacion de R13.
d10=0.098; %Distancia entre el centro de gravedad de la base y el punto de
aplicacion de R12.
% Desarrollo de las operaciones para hallar los valores de las
reacciones,fuerzas y momentos.
FT=T/R;
Wsx= Ws*sin (a);
Wsy= Ws*cos (a);
Waf1x= Waf1*sin (a);
Waf1y= Waf1*cos (a);
Waf2x= Waf2*sin (a);
Waf2y= Waf2*cos (a);
Wex = We*sin (a);
Wey = We*cos (a);
90
Wcx = Wc*sin (a);
Wcy = Wc*cos (a);
Wafix =Wafi*sin (a);
Wafiy =Wafi*cos (a);
Wbx= Wb*sin (a);
Wby= Wb*cos (a);
% Ecuaciones del diagrama de cuerpo libre 1
RA=(Waf2y+Wsy)/2
RB=RA
R1=Wsx
NFT=FT
Nac=NFT
Fen=Nac
M1= Fen*d3-Nac*d4
R2=FT
R3=Waf2x+R1+R2
R4=R3
% Ecuaciones del diagrama de cuerpo libre 2
RC =RA+(Waf1y/2)
RD=RC
R5=Waf1x+R4
M2=(-1)*(R5*d5)
RE =RC+ (Wey/2)
RF=RE
R6=Wex+R5
% Ecuaciones del diagrama de cuerpo libre 3
RG=Wcy/2
RH=RG
R11=FT+Wcx
RI=RE+RG+(Wafiy/2)
RJ=RI
91
R13=R6+Wafix
M3=((R6*d7)+(R13*d9))*(-1)
Ftor=R11
R12=R13
% Ecuaciones del diagrama de cuerpo libre 4
RK=RI+ (Wby/2)
RL=RK
FR=Wbx+R12
M4=(R12*d10)*(-1)
Es importante considerar que Matlab opera con el valor del ángulo en radianes
por tal motivo el ingreso de los ángulos se lo ha hecho en radianes.
El primer caso se lo utiliza para hallar los valores máximos de las reacciones en
la coordenada Y. El segundo caso se lo utiliza para hallar los valores máximos
de las reacciones y fuerzas en la coordenada X. Para la selección del valor de
los momentos se tomará el valor más alto que se obtenga después de realizar
los dos casos expuestos.
Reacciones en Y:
92
X 2004.5 I 2138.9
¥ 2138.9
Reacciones en X:
1 2940
Fuerzas:
480
Momentos:
1 4.32
2 −109.66
3 −224.31
4 −435.5
93
El módulo a utilizarse es de 2, y la relación de transmisión a utilizarse es de
3a1, estos conceptos fueron explicados en el presente capítulo en la sección
[Link].
El material del cuál van hacer los engranes y la corona es el AISI 4340, este
material fue escogido ya que sus propiedades permiten a los engranes y corona
tener una buena combinación entre dureza y tenacidad y también debido a que
este material tiene un gran rendimiento en piezas sometidas a ciclos de trabajo.
94
• La relación de transmisión. (de 3 a 1)
• El módulo de los engranes. (m=2)
• El ángulo de presión (20 grados)
• Si el sistema de engranaje es interno o externo. (interno)
• Número de dientes.
• Distancia entre centros.
• Ancho de cara del diente.
• Diámetro exterior.
• Diámetro de base.
• Diámetro primitivo.
• Diámetro de raíz.
• Addendum
• Deddendum
Parámetros Valores
Dp 132 [mm]
Z2 66
m 2
F 20 [mm]
[Link]ón 20°
Fuente: Propia.
95
Y la corona las siguientes características:
Fuente: Propia.
96
Es importante destacar que el programa Autodesk Inventor trabaja bajo varias
normas establecidas para el diseño de engranes, las mismas que se observa en
la figura 2.62.
Figura 2.62. Normas bajo las cuales trabaja Autodesk Inventor para el
diseño de engranajes.
Para conocer si el diseño de los engranes establecidos para este proyecto van
a fallar se realizaron cálculos de esfuerzo de contacto y de esfuerzo último de
tensión, y posteriormente estos valores obtenidos se comprueban con los
valores que entrega el software de diseño.
c£ l= ∙ l ∙ l ∙ l¤ ∙ lF
¨N ∙ 2.58
∙∙P l©
Donde:
• ¨N es el esfuerzo de tensión.
• c£ es la fuerza tangencial que se genera en el engrane debido al torque.
• l= es el factor de aplicación de la carga para resistencia a la flexión.
• l es el factor de tamaño para resistencia a la flexión.
• l es el factor de distribución de carga para resistencia a la flexión.
97
• l¤ es el factor de espesor de aro.
• lF es el factor de engrane intermedio o loco.
• es el ancho del diente.
• es el módulo del engrane.
• P es el factor geométrico para resistencia a la flexión.
• l© es el factor dinámico para resistencia a la flexión.
La tabla 2.16 muestra los valores del Factor de forma de Lewis, de la AGMA.
En base a la tabla siguiente se tiene que el factor de forma J para los piñones
cuyo número de dientes es de 66, J=0.427 (Valor obtenido por medio de
interpolación). Y para la corona cuyo número de dientes es de 200, J=0.466.
(Valor obtenido por medio de interpolación).
98
Factor dinámico Kv:
[« + ¬f 2 3]"
[2.59]
200
ªN =
« = 50 + 56 1 − [2.60]
12 − ¬f "/8
= [2.61]
4
99
[59.8 # 6 2 3]"
ªN =
200
ªN = 19.72 [/ ]
Cualquier desalineación axial o desviación axial en la forma del diente hará que
la carga transmitida Wt, quede no uniformemente distribuida sobre el ancho de
la cara del diente del engrane. Una manera aproximada y conservadora de
tomar en consideración una distribución de carga inferior a la uniforme es
aplicar el factor Km, a fin de incrementar los esfuerzos en anchos de cara más
grande.
5
NORTON Robert L., “Diseño de Máquinas”, P. 739
100
Tabla 2.17. Factores de distribución de carga Km.
Ancho de la cara
Km
en [mm]
6 (150) 1.7
9 (250) 1.8
101
Tabla 2.18. Factores de aplicación Ka.
Un engrane intermedio o loco está sujeto a la vez a más ciclos de esfuerzo por
unidad de tiempo, y a cargas alternantes más elevadas que sus primos de no
102
giro libre. Debido a que este sistema de engranajes no poseen engranes locos,
el factor Ki será 1.
O
cN = [2.62]
1
96 []
0.2 []
cN =
cN = 480 []
¨N = 81.754 [°=]
¨N = 74.912 [°=]
103
m= ∙ m ∙ m ∙ m ∙ c£
¨x = mR ∙ T [2.63]
m© ∙ ∙ ±} ∙ ²
Donde:
• ¨x es el esfuerzo de contacto.
• mR es el coeficiente elástico del material.
• m= es el factor de aplicación de la carga para resistencia al desgaste.
• m es el factor de tamaño para resistencia al desgaste.
• m es el factor de distribución de carga para resistencia al desgaste.
• m es el factor de terminado superficial.
• m© es el factor dinámico para resistencia al desgaste.
• ±} es el diámetro de paso del engrane.
• ² es el factor geométrico para resistencia al desgaste.
Este factor toma en consideración los radios de curvatura de los dientes del
engrane y el ángulo de presión. Para obtener el valor de I se tienen las
siguientes gráficas tanto para el piñón como la corona respectivamente:
6
MOTT Robert L., “Diseño de Elementos de Máquinas”, P. 455
104
Figura 2.65. Gráfica para hallar el factor geométrico superficial en la
corona7.
²} = 0.121
²e = 0.140
7
Ibíd., P. 455
105
El coeficiente elástico toma en consideración diferencias en el material del
diente y se lo puede obtener a partir del estándar 2001-B88 de la AGMA que
muestra la tabla 2.19.
En donde tanto el material del piñón como el de la corona son de acero, por
tanto se obtiene un coeficiente elástico m| = 191 [°=]³. .
Para los factores Ca, Cm, Cs y Cv se consideran los mismos valores obtenidos
para Ka, Km, Ks y Kv.
¨x = 399.335 [°=]
¨x = 213.267 [°=]
106
un factor de seguridad aproximado de 2.9 para los piñones y un factor de
seguridad de 5 para la corona.
107
un factor de seguridad para el esfuerzo de contacto de 3.181, este último valor
no muy similar al factor obtenido mediante los cálculos que es de 5, este error
se da debido a que el software de diseño asume un mismo valor de factor de
geometría superficial I = 0.119 tanto para los piñones como para la corona, y
en los cálculos se colocan valores distintos, consideración que es la más
efectiva.
Para el diseño del afuste inferior es importante conocer las medidas tanto de la
corona, como del rodamiento grande, ya que son estos elementos los que se
acoplaran interna y externamente en el afuste inferior. Con estas medidas, se
procede a diseñar la forma del afuste inferior, el cuál internamente debe poseer
un diámetro de igual tamaño al diámetro de la corona, el cuál es de 445.50
[mm], con sus respectivas tolerancias, las cuáles en este caso serán de h6J7.
En base a la longitud de los ejes de los piñones, y a la ubicación que debe tener
la corona para que los piñones engranen correctamente con la misma, es
requerido diseñar una grada en el afuste inferior que levante la corona hasta la
posición ideal. Esta ubicación de los engranes y corona, dependen de muchos
factores: del espesor del rodamiento grande, y en base a este espesor,
depende del tamaño del afuste superior, y de su ubicación con respecto al
afuste inferior, además depende de la posición de los pernos de acoplamiento.
Esta ubicación de los piñones y de la corona es primordial para poder tener un
sistema que funcione correctamente, y además para que sea sencillo de
ensamblar. De esta manera se obtiene que el tamaño de la grada requerida
para elevar la corona es de 46 [mm]. La grada tendrá 8 agujeros roscados M6
para los pernos de acoplamiento de la corona, y 2 agujeros de 6 y 8 [mm] de
diámetro respectivamente para las guías que servirán para colocar a la corona
siempre en la misma posición. El diámetro interno de la grada debe ser más
grande que el diámetro interno de la corona, para evitar de esta manera que los
dientes de la corona estén en contacto con la grada, por tal motivo se ha
108
considerado un diámetro interno de la grada de 405 [mm] tomando en cuenta
que el diámetro interno de la corona es de 396 [mm].
°
¨= [2.64]
«
«=¤∙ℎ [2.65]
22.55 []
¨=
445.5 [] ∙ 46[]
22.55 []
¨=
20493 []"
¨ = 1.1 [l°=]
109
Donde F= 22.55 [N], es el valor de RG y RH obtenidos en el punto 2.7.2 del
presente capítulo.
Por tal motivo la pared que forma el afuste inferior tiene un espesor de 14.5
[mm], valor que hay q comprobar que sea capaz de soportar el peso del
rodamiento grande para el caso extremo en el cual la plataforma se encuentre
virada a 90 grados. La reacción que actúan sobre esta pared es R6. De esta
manera se realiza un análisis por esfuerzo cortante:
O= [2.66]
«1
Donde r1 es el radio externo del cilindro igual a 230 [mm] y r2 es el radio interno
del cilindro igual a 222.75 [mm]. De esta manera se tiene:
4146.4 []
O = µ∙"8³¶ [··]¶ ,µ∙""".¡¶[··]¶
"
4146.4 []
5156.05 [" ]
O=
110
O = 804.18 [l°=]
Para ubicar al rodamiento grande en una altura necesaria para poder acoplar el
afuste inferior con el afuste superior, se ha diseñado una grada externa de un
diámetro de 475.5 [mm] con un valor menor al referido en las especificaciones
del rodamiento para su pista interna, y con una altura de 20 [mm], este valor de
la altura se ha considerado analizando varios parámetros, entre los cuales se
encuentran la ubicación de la corona y piñones, la ubicación del afuste superior
y su brida respectiva, y los pernos de acoplamiento y desacoplamiento del
afuste superior.
Para comprobar que la grada no vaya a fallar se realiza el análisis por esfuerzo
a compresión. En este caso las fuerzas que actúan sobre la grada externa del
afuste inferior, son las generadas por el peso del rodamiento grande, RE y
RF=1833.2 [N].
¨1 =
«2
«2 = ¤1 ∙ ℎ1
1833.2 []
¨1 =
475.5 [] ∙ 20[]
111
1833.2 []
¨1 =
9510 []"
¨1 = 192.76 [l°=]
Este disco deberá soportar las fuerzas de las reacciones por contacto con la
base, ejercidas principalmente por el propio peso del afuste inferior.
El valor escogido para el espesor del disco es de 10 [mm], esto debido a que
debe soportar un gran peso, y a que es un valor estándar de espesor para
placas de acero.
Con estos valores se procede a comprobar que el disco no falle por esfuerzo
cortante:
X
O1 =
«3
«3 ¤2 ∙ ℎ2
112
Donde b2 es la base igual a 638 [mm] y h2 es el espesor del disco igual a 10
[mm]. De esta manera se tiene:
2004.5 []
O1 =
638 [] ∙ 10[]
2004.5 []
O1 =
6380[]"
O1 = 314.18 [l°=]
8
O=1 =
«=1
DondeR13 es la reacción que se genera en los agujeros del afuste inferior y Aa1
es el área transversal al punto de aplicación de la fuerza sobre los agujeros.
Donde ba1 es la base del área transversal igual a 8.5 [mm] y ha1 es la altura
del área transversal igual a 10 [mm].
8
O=1 =
¤=1 ∙ ℎ=1
4443.9 []
O=1 =
8.5 [] ∙ 10[]
113
O=1 = 52.3 [°=]
8
O 1=
« 1
« 1=¤ 1∙ℎ 1
Donde be1 es la base del área transversal igual a 7.75 [mm] y he1 es la altura
del área transversal igual a 10 [mm].
8
O 1=
¤ 1∙ℎ 1
4443.9 []
O 1=
7.75 [] ∙ 10[]
O 1 = 57.34 [°=]
114
2.7.5. DISEÑO DE LA BASE.
La base tendrá una brida en la parte superior con un diámetro de 638 [mm] y un
espesor de 6 [mm] para colocar sobre la misma al afuste inferior.
Es importante que para ese espesor de la brida esta no vaya a fallar debido al
peso de todo el sistema que se colocará encima.
Las fuerzas generadas representativas de todo este peso son RI y RJ, con un
valor de 2004.5 [N].
El material del que va hacer la base será el mismo que del afuste inferior, el
AISI 1045, con un esfuerzo mínimo de 370 [MPa].
O2 =
«4
«4 = ¤3 ∙ ℎ3
115
2004.5 []
O2 =
638 [] ∙ 6 []
2004.5 []
O2 =
3828[]"
O2 = 523.64 [l°=]
Se observa que la brida superior no llegará a fallar. Esta brida constará con 16
agujeros pasantes con un diámetro de 8.5 [mm], para el acoplamiento con el
afuste inferior. Es importante comprobar que estos agujeros no vayan a fallar
por cortante debido a la presencia de la reacción R12 generada por los tornillos
de acople. Por lo cual se realiza el siguiente análisis:
"
O=2 =
«=2
Donde ba2 es la base del área transversal igual a 8.5 [mm] y ha2 es la altura
del área transversal igual a 6 [mm].
"
O=2 =
¤=2 ∙ ℎ=2
4443.9 []
O=2 =
8.5 [] ∙ 6[]
116
Se procede a comprobar que a la distancia que se encuentran los agujeros, el
espesor entre el agujero y el diámetro externo de la brida superior no falle por la
acción de R12.
"
O 2=
« 2
« 2=¤ 2∙ℎ 2
Donde be2 es la base del área transversal igual a 7.75 [mm] y he2 es la altura
del área transversal igual a 6 [mm].
"
O 2=
¤ 2∙ℎ 2
4443.9 []
O 2=
7.75 [] ∙ 6 []
O 2 = 95.56 [°=]
Además la base tendrá una brida en la parte inferior con 700 [mm] de diámetro,
10 [mm] de espesor y 16 agujeros de 12 [mm] de diámetro, esto para poder
acoplar a toda la plataforma ya sea al suelo, o algún camión de guerra.
Sobre esta brida inferior se generan reacciones para soportar todo el peso
general de la plataforma. Estas reacciones son RK y RL, ambas iguales de
2138.9 [N]. En base a estas reacciones, se comprueba que el espesor de la
brida inferior no vaya a fallar por esfuerzo cortante:
¥
O3
«5
117
Donde A5 es el área transversal al espesor de la brida inferior.
«5 = ¤4 ∙ ℎ4
2138.9 []
O3
700 [] ∙ 10 []
2138.9 []
O3
7000 []"
O3 305.55 [l°=]
¦
O=3 =
«=3
Donde ba3 es la base del área transversal igual a 12 [mm] y ha3 es la altura del
área transversal igual a 10 [mm].
¦
O=3 =
¤=3 ∙ ℎ=3
4468 []
O=3 =
12 [] ∙ 10[]
118
O=3 = 37.23 [°=]
¦
O 3=
« 3
« 3=¤ 3∙ℎ 3
Donde be3 es la base del área transversal igual a 12 [mm] y he3 es la altura del
área transversal igual a 10 [mm]. Como se observa que el valor de Ae3=Aa3, y
la fuerza FR actúa en ambas situaciones, el valor de O 3 sera:
O 3 = 37.23 [°=]
119
mismo escogido para los elementos anteriores, el AISI 1045, con un valor de
esfuerzo mínimo de 370 [MPa].
Para comprobar que la tapa no vaya a fallar se analiza a la misma por esfuerzo
de compresión, teniendo así:
¨2 =
«6
«6 = ¤5 ∙ ℎ5
1589.2 []
¨2 =
590 [] ∙ 12[]
¨2 = 224.46 [l°=]
8
O=4 =
«=4
120
Donde Aa4 es el área transversal al punto de aplicación de la fuerza sobre los
agujeros.
Donde ba4 es la base igual a 8.5 [mm] y ha4 es la altura igual a 12 [mm].
8
O=4 =
¤=4 ∙ ℎ=4
3658.3 []
O=4 =
8.5 [] ∙ 12[]
8
O 4=
« 4
« 4=¤ 4∙ℎ 4
Donde be4 es la base del área transversal igual a 13.75 [mm] y he4 es la altura
del área transversal igual a 12 [mm]. De esta manera se tiene:
8
O 4=
¤ 4∙ℎ 4
3658.3 []
O 4=
13.75 [] ∙ 12[]
121
O 4 = 22.17 [°=]
Además la tapa del afuste superior tendrá 6 agujeros pasantes de 10.5 [mm] de
diámetro ubicados a una distancia de 185 [mm] del centro, que servirán para
acoplarla a la estructura superior mediante pernos.
O=5 =
«=5
Donde ba5 es la base igual a 10.5 [mm] y ha5 es la altura igual a 12 [mm].
O=5 =
¤=5 ∙ ℎ=5
2940 []
O=5 =
10.5 [] ∙ 12[]
La tapa poseerá tres agujeros a una distancia de 134 [mm] del centro, dos de
ellos con un diámetro de 37 [mm] con una profundidad de 7 mm para albergar a
los rodamientos pequeños, y a partir de esta superficie, se generará un agujero
de 34 [mm] pasante para que se pueda colocar los ejes de los piñones locos, y
el otro agujero tendrá un diámetro de 68 [mm] y será pasante, este agujero es
122
más grande debido a que a través de este pasará el acople que une el eje de la
caja de transmisión al eje del piñón que va a transmitir el movimiento.
"
O=6 =
«=6
"
O=6 =
¤=6 ∙ ℎ=6
480 []
O=6 =
37 [] ∙ 7[]
El diseño final de la tapa del afuste superior se observa en la figura 2.28 del
capítulo [Link].
La parte inferior del afuste superior será un cilindro hueco de 590 [mm] de
diámetro exterior y 580 [mm] de diámetro interior, este valor viene dado por el
diámetro externo del rodamiento grande. El afuste superior será del mismo
123
material que las descritas anteriormente, el AISI 1045, con un esfuerzo mínimo
de 370 [MPa].
Para comprobar que esta pared del cilindro de 10 [mm] no va a fallar por la
carga que recae sobre ella debido al peso del rodamiento, se realiza el
siguiente análisis por esfuerzo cortante:
O4 =
«7
B ∙ 8 " − B ∙ "
«7 =
2
Donde r3 es el radio externo del cilindro actual de la parte inferior del afuste
superior, y r4 es el radio interno de la parte inferior del afuste superior, teniendo
un espesor de la pared de 10 [mm]. De esta manera se tiene:
3821 []
O4 = µ∙"¢¶ ,µ∙"¢³¶
"
3821 []
O4 =
4594.6[" ]
O4 = 831.63 [l°=]
El afuste tendrá además una brida interna en la parte superior, para poder
empernar sobre esta a la tapa del afuste superior. La brida tendrá un diámetro
interno de 530 [mm] y 26 [mm] de espesor este valor de espesor es necesario,
ya que es importante que los pernos tengan espacio en donde enroscar y
además puedan servir de ayuda para el desmontaje del rodamiento grande. Los
530 [mm] de radio interno, dan a la brida una longitud de 30 [mm], este valor es
el considerado para que no llegase a fallar la brida en cualquiera de los
124
agujeros de los pernos. La brida tendrá 16 agujeros M8 roscados para poder
empernar la tapa del afuste superior. Es necesario comprobar que la brida no
va a llegar a fallar por lo que se realiza a continuación un análisis de esfuerzo a
compresión.
¨3 =
«8
«8 = ¤6 ∙ ℎ6
1589.2 []
¨3 =
590 [] ∙ 26[]
¨3 = 103.59 [l°=]
Se comprueba que los agujeros no vayan a fallar por esfuerzo cortante por la
presencia de R4:
O=7 =
«=7
O=7 =
¤=7 ∙ ℎ=7
125
3658.3 []
O=7 =
6.65 [] ∙ 26[]
O 5=
« 5
« 5=¤ 5∙ℎ 5
Donde be5 es la base del área transversal igual a 14.67 [mm] y he5 es la altura
del área transversal igual a 26 [mm]. De esta manera se tiene:
O 5=
¤ 5∙ℎ 5
3658.3 []
O 5=
14.67 [] ∙ 26[]
O 5 = 9.6 [°=]
La altura que tendrá la parte inferior del afuste superior es de 102 [mm], este
valor es el ideal para poder tener un sistema de engranajes colocado en la
correcta posición.
El diseño final del afuste superior se observa en la figura 2.27 del capítulo
[Link].
126
2.7.7. DISEÑO DE EJES.
Es conocido que dos de los ejes de los engranes deben estar conformados por
3 secciones, la primera con el menor diámetro que es la que estará en contacto
con el acople, la segunda sección es en donde se colocará el rodamiento, y
esta debe tener el diámetro del mismo tamaño que el diámetro interno del
rodamiento. Y la tercera sección será en donde se sostendrá el rodamiento en
su pista interna.
127
El software devuelve que el valor mínimo que debe tener el diámetro del eje
para no fallar debido a las fuerzas ejercidas que se observan en la sección 2.7.2
del presente capítulo es de 18.5237 [mm].
9
9 8
32 ∙ ¹ h ∙ º 3 dº
" " "
> = ¸ ∙ 5- . + ∙/ 0 6 ¼ [2.68]
B ¹ 4 ¹»
Donde:
Sy: es el esfuerzo último de fluencia del material, en este caso 855 [MPa]
correspondiente al (AISI 4340), material recomendado para árboles y
engranajes.
Donde Sn’ es el límite de resistencia a la fatiga y es igual a 0.5 Sut para aceros
que tengan un valor de Sut menor a los 1400 [MPa].
8
NORTON Robert L., “Diseño de Máquinas”, P. 999
128
Para los cálculos se utiliza el acero AISI 4340 templado y revenido a 1200°F, ya
que es el que menor resistencia tiene con diferencia a los diferentes aceros AISI
4340 sometidos a otros procesos, esto con el fin de tener consideraciones
extremas.
Cm: es el factor del material, de acuerdo al proceso del cual se obtuvo, en este
caso 1 dado que es un acero de alta calidad, forjado y con varios tratamientos
térmicos.
Cada uno de los factores se explica con mayor claridad en la sección 2.7.9 en el
análisis por fatiga de los pernos.
Es necesario hallar el momento máximo al que estará sujeto el eje, por tal
motivo se realiza un análisis por esfuerzo cortante y momento flector como se
describe en el siguiente programa desarrollado en Matlab.
129
MA=(F1*d1)+(F2*d2);
MB=MA+(F1*0);
MC=MB+(F1*(d-d1));
M1(1:93)=0;
V1(1:93)=0;
x=[0.001:d2/96:d2];
for (d=0:0.001:0.094)
if (d <= d1)
M1(i)=(F1*d1)+(F2*d2);
V1(i)=0;
i=i+1;
end
if (d>d1)&&(d<=d2)
M1(i)=MB+(F1*(d-d1));
V1(i)=F1;
i=i+1;
end
if (d >= 0.074)
M1(i)=0;
V1(i)=0;
i=i+1;
end
end
figure (1);subplot (2,2,1);plot (x,V1);
xlabel('X [m]');
ylabel('V [N]');
title('Esfuerzo Cortante en el eje');
grid on;
figure (1);subplot (2,2,2);plot (x,M1);
xlabel('X [m]');
ylabel('M [Nm]');
130
title('Momento flector en el eje');
grid on;
¹ = 386.48 [°=]
l = 1 + ; ∙ l£ − 1 [2.69]
131
Figura 2.69. Diagrama de sensibilidad a la muesca para aceros9.
El radio del filete del hombre del eje es de 1 [mm], Sut es de 965 [MPa], con
estos valores nos dirigimos a la gráfica anterior y se obtiene el valor de q, el
cual dio un valor aproximado de ; y 0.65.
Para este caso se tiene una flecha sujeta a flexión, de esta manera el valor de
kt se lo obtiene en base a la figura 2.70.
132
Figura 2.70. Factor de concentración de esfuerzos geométrico Kt para una
flecha con un filete de hombro a flexión10.
l = 1 + 0.65 ∙ 2 − 1
l = 1.65
% Programa desarrollado para el cálculo del diámetro del eje de los piñones
% en orientación.
FS=5; % Ingreso del factor de seguridad
kf=1.65; % Ingreso del factor de resistencia a la fatiga
Mmax=33.12; % Ingreso del momento máximo en Nm
133
Tmax=32; % Ingreso del torque máximo en Nm
Sn=386.48; % Ingreso de la resistencia a la fatiga en MPa
Sy=855; % Esfuerzo último de fluencia del material. Acero AISI 4340
% Calculo del diámetro:
d=(((32*FS/3.1416)*(((kf*Mmax*1000/Sn)^2)+((3/4)*((Tmax*1000/Sy)^2)))^(1/2))
^(1/3))
134
mayor tamaño. Como se observa en los resultados obtenidos en el diseño del
afuste inferior, la pared del mismo no va a fallar, por tal motivo el primer
parámetro para seleccionar el rodamiento grande está aprobado, su diámetro
interno será de 460 [mm], no se colocará un rodamiento con un diámetro interno
de mayor tamaño, ya que se estaría sobredimensionando el sistema, y esto
implicaría más gastos, y un mayor peso innecesario. Una vez conocido que se
requiere un rodamiento con un diámetro interno de 460 [mm], es importante
seleccionar de que valor se necesita el diámetro externo, en este caso como es
irrelevante este valor, se escogerá el de menor valor, para de esta manera
optimizar el tamaño tanto del afuste inferior y del afuste superior, solo en el
caso extremo de que la superficie de contacto de uno de los afustes con el
rodamiento llegaría a fallar porque es muy pequeña, se consideraría un
rodamiento con un diámetro externo más grande para así tener mayor
superficie de contacto. Como se observa en el diseño de los afustes, ninguno
de ellos llegará a fallar, por tal motivo se ha seleccionado el rodamiento con un
diámetro externo de 580 [mm].
11
< [Link]
135
Figura 2.71. Cálculo del valor de P en el software en línea de la SKF.
136
Como estos rodamientos solo reciben carga radial, el valor de P será el mismo
valor que la carga radial. La carga real que recae en los rodamientos es de 160
[N]. Como en total son 2 rodamientos, cada uno de ellos recibirá una carga de
80 [N] radialmente. Este valor de carga es menor a los 125 [N] que puede
soportar el rodamiento, de tal manera estos trabajarán sin ningún problema.
Para este análisis se utiliza la carga real, ya que al utilizar un factor de
seguridad los rodamientos llegarían a fallar, y necesitaríamos rodamientos más
grandes, en este caso debido a la altura del rodamiento estos rodamiento no
entrarían correctamente en la tapa del afuste superior, y no existiría una
superficie para que se puedan asentar. Es por eso que fue necesario buscar el
rodamiento más óptimo que permita esta situación. En el apéndice A.11 se
muestra las dimensiones del rodamiento 61805-2RS1.
Para acoplar los piñones con sus respectivos ejes se van a utilizar los acoples
del tipo B-LOC 109. El mismo que se observa en la figura 2.72.
12
Catálogo “Keyless Bushings Fenner Drives”, P. 17
137
El principal parámetro para seleccionar el acople es el tamaño del diámetro del
eje, en este caso el eje posee un diámetro de 24 [mm], para lo cual buscamos
en el catálogo de la Keyless Bushings un acople que cumpla con esta
especificación. (Es importante destacar que estos acoples tienen una tolerancia
de +-0.025 [mm]). Una vez ubicado el posible acople, es necesario comprobar
el segundo parámetro, este es que el diámetro externo D del acople quepa sin
ningún problema dentro del engrane, una vez comprobado esto, se comprueba
que la fuerza que va actuar sobre el acople no sea mayor a la fuerza expuesta
en el catálogo. En el apéndice A.14 se muestra el catálogo con los diversos
tipos de acoples B-109.
La fuerza que va actuar sobre el acople es de 480 [N], mientras que la fuerza
máxima permisible sobre el acople es de 31219 [N], como se puede observar el
acople no va a sufrir ningún tipo de daño por la fuerza aplicada.
Para esta aplicación existen los acoples B-LOC de la serie WK que permiten
acoplar dos ejes sin chavetas. El acople se muestra en la figura 2.73.
138
Figura 2.73. Acople B-LOC WK13.
13
Ibíd., P. 48
139
Para los tornillos de grado 4,6 se tiene que el esfuerzo de fluencia Sy es de 240
[MPa] y el esfuerzo último Sut es de 400 [MPa].
Para obtener esta fuerza se asume el caso extremo en el que la plataforma está
a 90°, ya que es para ese caso en donde se generan los mayores valores de
esfuerzos por cortante en los pernos. Para seleccionar la fuerza de mayor valor
se observan los valores obtenidos en la sección 2.7.2 del presente capítulo,
obteniendo que la mayor fuerza generada es de 4443,9 [N], generada en el
perno de acoplamiento del afuste inferior con la base.
¹MN
O= [2.70]
400 [°=]
O=
5
O = 80 [°=]
140
Como
ª ª
O= = µ∙¶ [2.71]
«
4∙ª
±=T
B∙O
Obteniendo así:
4443.9 []
ª= =
6 6
ª = 740.65 []
4 ∙ 740.65 [N]
±=T
B ∙ 80 ¾¶ À
¿
± y 3.43 []
141
Análisis por esfuerzo de tensión:
Para hallar el valor del diámetro, se utiliza el valor del esfuerzo de fluencia y un
factor de seguridad, ya que para este caso se obtendrá el valor mínimo del
diámetro requerido para que el tornillo no falle.
¹»
¨= [2.72]
240 [°=]
¨=
5
¨ = 48 [°=]
Como
£ £
¨= = µ∙¶ [2.73]
«
4 ∙ S
±=T
B∙´
S = 222.97 []
4 ∙ 222.97 [N]
±=T
B ∙ 48 ¾ À
¿
¶
± y 2.43 []
Para el análisis por fatiga se analiza la figura 2.74 donde se muestra la gráfica
de esfuerzo-tiempo para esta aplicación, y en esta gráfica se puede observar el
valor del esfuerzo de tensión alternante ´= y el valor del esfuerzo medio de
tensión ´. Estos valores nos permitirán realizar un análisis por fatiga por el
método de Goodman.
¨º − ¨
¨ = [2.74]
2
14
NORTON Robert L., “Diseño de Máquinas”, P. 356
143
¨º + ¨
¨ = [2.75]
2
Sº
¨º = [2.76]
«
S
¨ = [2.77]
«
Sº − S
¨ =
2∙«
Sº + S
¨ =
2∙«
Sº − S
¨ = "∙µ∙ ¶
Sº + S
¨ = "∙µ∙ ¶
2 ∙ Sº
¨ =
B ∙ ±"
144
2 ∙ Sº
¨ =
B ∙ ±"
¨ ¨ 1
+ = [2.78]
¹ ¹MN
Donde:
Donde Sn’ es el límite de resistencia a la fatiga y es igual a 0.5 Sut para aceros
que tengan un valor de Sut menor a los 1400 [MPa].
Cm: es el factor del material, de acuerdo al proceso utilizado para fabricar los
tornillos. Se tiene que:
m = = ∙ ¹MN Ä [2.80]
145
Como el proceso de fabricación de los tornillos es por medio de laminación, se
tiene que el valor de a es de 57.7 [MPa] y el valor del exponente b es de -0.718.
m = 57.7 ∙ 400,³.¡Å
m = 0.781
Cr: es el factor de confiabilidad del material, para este caso se utiliza un valor de
0.75 ya que se han considerado tornillos de la más baja calidad.
> ,³.88
m = [2.81]
7.62
Para este caso se asumirá un diámetro mayor a los valores obtenidos en los
análisis como diámetros mínimos, por tal motivo para este análisis se considera
un diámetro =4 [mm]. Y se obtiene:
4 ,³.88
m =
7.62
m = 1.075
¹ = 100.75 [°=]
146
222.97 [] 222.97 [] 1
+ =
100.75 [°=] ∙ B ∙ ± " 400 [°=] ∙ B ∙ ± " 5
± y 2.1 []
147
El sistema de elevación de la plataforma consta de las siguientes partes a las
que se hará referencia para el análisis de cargas y posterior diseño de
elementos.
1) Brazo izquierdo.
2) Soporte de misiles igla izquierdo.
3) Engrane helicoidal.
4) Buje izquierdo.
5) Árbol de transmisión.
6) Caja de cámaras.
7) Soporte de misiles igla derecho.
8) Brazo derecho.
9) Buje derecho.
10) Eje derecho.
2.8.2 ANÁLISIS DE CARGAS ESTÁTICAS.
Figura 2.76. Diagrama del cuerpo libre del misil Igla del lado izquierdo de
la plataforma.
148
Aplicando las ecuaciones de equilibrio obtenemos las reacciones.
= 0
» + » + K sin = 0
Æ + Æ + K cos = 0
= 0
2 ∙ » ∙ % − K cos ∙ $ − Ç ∙ $ + $ = 0 [2.82]
149
= 0
» = » + » + K sin
Æ = Æ + Æ + K cos [2.83]
= 0
º = − » ∙ % − Æ ∙ $ − » ∙ % + Æ ∙ $ − » ∙ % + Æ ∙ $ + º
+ º
» = Æ ∙ ! + º ∙ % + Æ ∙ ! + º ∙ % + Æ ∙ ! + º ∙ % + » + »
Æ = » ∙ ! − º ∙ $ + » ∙ ! + º ∙ $ + » ∙ ! + º ∙ $ + Æ + Æ
Figura 2.78. Diagrama del cuerpo libre del conjunto ejes - cámaras.
150
Mediante las ecuaciones de equilibrio estático obtenemos las reacciones en los
rodamientos.
= 0
º =
= 0
O = 2 ∙ º
−N ∙ ! − !R + K sin ∙ ! − !È − » ∙ ! − 2 ∙ !
» = [2.86]
! − !
−Ã ∙ ! − !R + K cos ∙ ! − !È − » ∙ ! − 2 ∙ !
Æ = [2.87]
! − !
151
O
N = [2.88]
= N ∙ tan É [2.89]
N ∙ tan Ê
à = [2.90]
cos É
= 0
º = º
» = » # K sin
Æ = Æ # K cos
152
= 0
» = Æ ∙ ! # Æ ∙ !
º = » ∙ ! # » ∙ !
= 0
º =
» = » # N 2 » 2 K sin
153
Æ = Æ # Ã 2 Æ 2 K cos
= 0
» = Ã ∙ !| 2 Æ ∙ ! 2 Æ ∙ ! # Æ ∙ !
Æ = 2N ∙ !| # » ∙ ! # » ∙ ! 2 » ∙ !
clear
clc
%%Posición del robot en rad.
theta1=0;
theta2=20*pi/180:(50*pi/180)/100:70*pi/180;
%%Dimensiones del IGLA en m.
xa=0;
ya=(236.769-25)/1000;
za=99.080/1000;
xb=xa;
yb=(615.231+25)/1000;
zb=za;
%%Peso del conjunto de 2 IGLA con soportes en N.
%%Se toma en cuenta un 25% de peso adicional por motivos de acoples o
%%elementos adicionales que se puedan acoplar para la automatización.
Wigla=1.25*2*36*9.81;
%%Cálculo de las reacciones sobre el brazo en N.
Rax(1:101)=0;
Ray=Wigla.*sin(theta2)/2;
154
Raz=-(2*Ray*za-Wigla*yb.*cos(theta2))/(ya+yb);
Rbx(1:101)=0;
Rby=Ray;
Rbz=Wigla.*cos(theta2)-Raz;
May(1:101)=0;
Max(1:101)=0;
Maz(1:101)=0;
Mby(1:101)=0;
Mbx(1:101)=0;
Mbz(1:101)=0;
Ra=sqrt(Rax.^2+Ray.^2+Raz.^2);
Rb=sqrt(Rbx.^2+Rby.^2+Rbz.^2);
Ma=sqrt(Max.^2+May.^2+Maz.^2);
Mb=sqrt(Mbx.^2+Mby.^2+Mbz.^2);
%%Gráficas de las reacciones del IGLA.
figure(1);subplot(4,2,1);plot(theta2*180/pi,Rax);
xlabel('Theta 2 [deg]');
ylabel('Rax [N]');
grid on
figure(1);subplot(4,2,3);plot(theta2*180/pi,Ray);
xlabel('Theta 2 [deg]');
ylabel('Ray [N]');
grid on
figure(1);subplot(4,2,5);plot(theta2*180/pi,Raz);
xlabel('Theta 2 [deg]');
ylabel('Raz [N]');
grid on
figure(1);subplot(4,2,2);plot(theta2*180/pi,Rbx);
xlabel('Theta 2 [deg]');
ylabel('Rbx [N]');
grid on
155
figure(1);subplot(4,2,4);plot(theta2*180/pi,Rby);
xlabel('Theta 2 [deg]');
ylabel('Rby [N]');
grid on
figure(1);subplot(4,2,6);plot(theta2*180/pi,Rbz);
xlabel('Theta 2 [deg]');
ylabel('Rbz [N]');
grid on
figure(1);subplot(4,2,7);plot(theta2*180/pi,Ra);
xlabel('Theta 2 [deg]');
ylabel('Ra [N]');
grid on
figure(1);subplot(4,2,8);plot(theta2*180/pi,Rb);
xlabel('Theta 2 [deg]');
ylabel('Rb [N]');
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Reacciones IGLA');
set(gca,'Visible','off');
set(h,'Visible','on');
%%Dimensiones del brazo en m.
xa=54.763/1000;
ya=355.617/1000;
za=(52.2-35)/1000;
xb=xa;
yb=544.383/1000;
zb=za;
xc=47.737/1000;
156
yc=55.617/1000;
zc=17.8/1000;
%%Peso del brazo en N.
%%Se toma en cuenta un 25% que pudiera excederse por motivos de
%%manufactura.
Wbrazo=1.25*10*9.81;
%%Ecuaciones de equilibrio del brazo en N, m.
Rcy=Ray+Rby+Wbrazo.*sin(theta2);
Rcx=Rax+Rbx;
Rcz=Raz+Rbz+Wbrazo.*cos(theta2);
Mcx=-Rax*za-Raz*ya-Rbx*zb+Rbz*yb-Rcy*zc+Rcz*yc+Max+Mbx;
Mcy=Raz*xa+Rax*za+Rbz*xb+Rbx*zb+Rcz*xc+Rcx*zc+May+Mby;
Mcz=Rax*xa-Ray*ya+Rbx*xb+Rbx*yb+Rcy*xc+Rcx*yc+Maz+Mbz;
Rc=sqrt(Rcx.^2+Rcy.^2+Rcz.^2);
Mc=sqrt(Mcx.^2+Mcy.^2+Mcz.^2);
%%Gráficas de las reacciones y momentos del brazo.
figure(2);subplot(4,2,1);plot(theta2*180/pi,Rcx);
xlabel('Theta 2 [deg]');
ylabel('Rcx [N]');
grid on
figure(2);subplot(4,2,3);plot(theta2*180/pi,Rcy);
xlabel('Theta 2 [deg]');
ylabel('Rcy [N]');
grid on
figure(2);subplot(4,2,5);plot(theta2*180/pi,Rcz);
xlabel('Theta 2 [deg]');
ylabel('Rcz [N]');
grid on
figure(2);subplot(4,2,2);plot(theta2*180/pi,Mcx);
xlabel('Theta 2 [deg]');
ylabel('Mcx [N*m]');
157
grid on
figure(2);subplot(4,2,4);plot(theta2*180/pi,Mcy);
xlabel('Theta 2 [deg]');
ylabel('Mcy [N*m]');
grid on
figure(2);subplot(4,2,6);plot(theta2*180/pi,Mcz);
xlabel('Theta 2 [deg]');
ylabel('Mcz [N*m]');
grid on
figure(2);subplot(4,2,7);plot(theta2*180/pi,Rc);
xlabel('Theta 2 [deg]');
ylabel('Rc [N]');
grid on
figure(2);subplot(4,2,8);plot(theta2*180/pi,Mc);
xlabel('Theta 2 [deg]');
ylabel('Mc [N]');
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Reacciones en el Brazo');
set(gca,'Visible','off');
set(h,'Visible','on');
%%Dimensiones del conjunto ejes-cámara en m.
xp=50.000/1000;
xd=102.500/1000;
xw=484.500/1000;
xe=866.500/1000;
xc=969.000/1000;
yw=0;
zw=0;
158
%%Peso del conjunto ejes-cámara en N.
%%Se toma en cuenta un 25% que pudiera excederse en la manufactura o
%%elementos adicionales necesarios para la sujeción de las cámaras.
Wej_cam=1.25*(88.080+2*8.126)*9.81;
%%Cálculo del torque estático del motor.
thau2=2*Mcx+Wej_cam*yw.*cos(theta2)+Wej_cam*zw.*sin(theta2);
%%Gráfica del torque estático de inclinación.
figure(3);plot(theta2*180/pi,thau2);
figure(3);title('Torque Estático de Inclinación');
xlabel('Theta 2 [deg]');
ylabel('Thau 2 [N*m]');
grid on
%%Transmisión Vertical.
Fi=20*pi/180; %%Ángulo de hélice en rad.
Psi=20*pi/180; %%Ángulo de presión en rad.
R1=172.667/2000; %%Radio del engrane 1 en m.
R2=339.380/2000; %%Radio del engrane 2 en m.
R3=238.162/2000; %%Radio del engrane 3 en m.
%%Fuerzas en los dientes.
Ft=thau2./R3;
Fa=Ft.*tan(Psi);
Fr=Ft.*tan(Fi)./cos(Psi);
%%Ecuaciones de equilibrio del conjunto ejes-cámara en N, m.
Rdx=Fa;
Rdy=(-Rcy*(xc-2*xe)+Wej_cam.*sin(theta2)*(xe-xw)-Ft*(xe-xp))./(xe-xd);
Rdz=(-Rcz*(xc-2*xe)+Wej_cam.*cos(theta2)*(xe-xw)-Fr*(xe-xp))./(xe-xd);
Rex(1:101)=0;
Rey=Wej_cam.*sin(theta2)+2*Rcy-Ft-Rdy;
Rez=Wej_cam.*cos(theta2)+2*Rcz-Fr-Rdz;
Rd=sqrt(Rdx.^2+Rdy.^2+Rdz.^2);
Re=sqrt(Rex.^2+Rey.^2+Rez.^2);
159
%%Gráficas de las reacciones y momentos del conjunto ejes-cámara.
figure(4);subplot(4,2,1);plot(theta2*180/pi,Rdx);
xlabel('Theta 2 [deg]');
ylabel('Rdx [N]');
grid on
figure(4);subplot(4,2,3);plot(theta2*180/pi,Rdy);
xlabel('Theta 2 [deg]');
ylabel('Rdy [N]');
grid on
figure(4);subplot(4,2,5);plot(theta2*180/pi,Rdz);
xlabel('Theta 2 [deg]');
ylabel('Rdz [N]');
grid on
figure(4);subplot(4,2,2);plot(theta2*180/pi,Rex);
xlabel('Theta 2 [deg]');
ylabel('Rex [N*m]');
grid on
figure(4);subplot(4,2,4);plot(theta2*180/pi,Rey);
xlabel('Theta 2 [deg]');
ylabel('Rey [N*m]');
grid on
figure(4);subplot(4,2,6);plot(theta2*180/pi,Rez);
xlabel('Theta 2 [deg]');
ylabel('Rez [N*m]');
grid on
figure(4);subplot(4,2,7);plot(theta2*180/pi,Rd);
xlabel('Theta 2 [deg]');
ylabel('Rd [N]');
grid on
figure(4);subplot(4,2,8);plot(theta2*180/pi,Re);
xlabel('Theta 2 [deg]');
160
ylabel('Re [N]');
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Reacciones en el Conjunto Ejes-Cámara');
set(gca,'Visible','off');
set(h,'Visible','on');
%%Dimensiones de la brida en m.
xd=13.887/1000;
xf=12.311/1000;
%%Peso de la brida en N.
%%Se toma en cuenta un 25% adicional debido a posibles cambios de los
%%materiales en la manufactura.
Wbrida=1.25*0.81*9.81;
%%Ecuaciones de equilibrio de la brida izquierda en N, m.
Rfx=Rdx;
Rfy=Rdy+Wbrida.*sin(theta2);
Rfz=Rdz+Wbrida.*cos(theta2);
Mfx=Rfy*xf+Rdy*xd;
Mfy=Rfz*xf+Rdz*xd;
Mfz(1:101)=0;
Rf=sqrt(Rfx.^2+Rfy.^2+Rfz.^2);
Mf=sqrt(Mfx.^2+Mfy.^2+Mfz.^2);
%%Ecuaciones de equilibrio de la brida derecha en N, m.
Rhx=Rex;
Rhy=Rey+Wbrida.*sin(theta2);
Rhz=Rez+Wbrida.*cos(theta2);
Mhx=Rhy*xf+Rey*xd;
Mhy=Rhz*xf+Rez*xd;
Mhz(1:101)=0;
161
Rh=sqrt(Rhx.^2+Rhy.^2+Rhz.^2);
Mh=sqrt(Mhx.^2+Mhy.^2+Mhz.^2);
%%Gráficas de las reacciones y momentos de la brida izquierda.
figure(5);subplot(4,2,1);plot(theta2*180/pi,Rfx);
xlabel('Theta 2 [deg]');
ylabel('Rfx [N]');
grid on
figure(5);subplot(4,2,3);plot(theta2*180/pi,Rfy);
xlabel('Theta 2 [deg]');
ylabel('Rfy [N]');
grid on
figure(5);subplot(4,2,5);plot(theta2*180/pi,Rfz);
xlabel('Theta 2 [deg]');
ylabel('Rfz [N]');
grid on
figure(5);subplot(4,2,2);plot(theta2*180/pi,Mfx);
xlabel('Theta 2 [deg]');
ylabel('Mfx [N*m]');
grid on
figure(5);subplot(4,2,4);plot(theta2*180/pi,Mfy);
xlabel('Theta 2 [deg]');
ylabel('Mfy [N*m]');
grid on
figure(5);subplot(4,2,6);plot(theta2*180/pi,Mfz);
xlabel('Theta 2 [deg]');
ylabel('Mfz [N*m]');
grid on
figure(5);subplot(4,2,7);plot(theta2*180/pi,Rf);
xlabel('Theta 2 [deg]');
ylabel('Rf [N]');
grid on
162
figure(5);subplot(4,2,8);plot(theta2*180/pi,Mf);
xlabel('Theta 2 [deg]');
ylabel('Mf [N]');
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Reacciones en la Brida Izquierda');
set(gca,'Visible','off');
set(h,'Visible','on');
%%Gráficas de las reacciones y momentos de la brida derecha.
figure(6);subplot(4,2,1);plot(theta2*180/pi,Rhx);
xlabel('Theta 2 [deg]');
ylabel('Rhx [N]');
grid on
figure(6);subplot(4,2,3);plot(theta2*180/pi,Rhy);
xlabel('Theta 2 [deg]');
ylabel('Rhy [N]');
grid on
figure(6);subplot(4,2,5);plot(theta2*180/pi,Rhz);
xlabel('Theta 2 [deg]');
ylabel('Rhz [N]');
grid on
figure(6);subplot(4,2,2);plot(theta2*180/pi,Mhx);
xlabel('Theta 2 [deg]');
ylabel('Mhx [N*m]');
grid on
figure(6);subplot(4,2,4);plot(theta2*180/pi,Mhy);
xlabel('Theta 2 [deg]');
ylabel('Mhy [N*m]');
grid on
163
figure(6);subplot(4,2,6);plot(theta2*180/pi,Mhz);
xlabel('Theta 2 [deg]');
ylabel('Mhz [N*m]');
grid on
figure(6);subplot(4,2,7);plot(theta2*180/pi,Rh);
xlabel('Theta 2 [deg]');
ylabel('Rh [N]');
grid on
figure(6);subplot(4,2,8);plot(theta2*180/pi,Mh);
xlabel('Theta 2 [deg]');
ylabel('Mh [N]');
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Reacciones en la Brida Derecha');
set(gca,'Visible','off');
set(h,'Visible','on');
%%Dimensiones del árbol de transmisión en m.
xc=97.014/1000;
xd=4.236/1000;
xg=43.486/1000;
xp=4.236/1000;
%%Peso del árbol de transmisión en N.
%%Se toma en cuenta un 25% adicional debido a posibles cambios de los
%%materiales en la manufactura.
Warbol=1.25*8.126*9.81;
%%Ecuaciones de equilibrio del árbol de transmisión en N, m.
Rgx(1:101)=0;
Rgy=Rdy+Ft-Rcy-Warbol.*sin(theta2);
Rgz=Rdz+Fr-Rcz-Warbol.*cos(theta2);
164
Mgx(1:101)=0;
Mgy=Fr.*xp-Rcz.*xc-Rdz.*xd+Rgz.*xg;
Mgz=-Ft.*xp+Rcy.*xc+Rdy.*xd-Rgy.*xg;
Rg=sqrt(Rgx.^2+Rgy.^2+Rgz.^2);
Mg=sqrt(Mgx.^2+Mgy.^2+Mgz.^2);
%%Gráficas de las reacciones y momentos del árbol de transmisión.
figure(7);subplot(4,2,1);plot(theta2*180/pi,Rgx);
xlabel('Theta 2 [deg]');
ylabel('Rgx [N]');
grid on
figure(7);subplot(4,2,3);plot(theta2*180/pi,Rgy);
xlabel('Theta 2 [deg]');
ylabel('Rgy [N]');
grid on
figure(7);subplot(4,2,5);plot(theta2*180/pi,Rgz);
xlabel('Theta 2 [deg]');
ylabel('Rgz [N]');
grid on
figure(7);subplot(4,2,2);plot(theta2*180/pi,Mgx);
xlabel('Theta 2 [deg]');
ylabel('Mgx [N*m]');
grid on
figure(7);subplot(4,2,4);plot(theta2*180/pi,Mgy);
xlabel('Theta 2 [deg]');
ylabel('Mgy [N*m]');
grid on
figure(7);subplot(4,2,6);plot(theta2*180/pi,Mgz);
xlabel('Theta 2 [deg]');
ylabel('Mgz [N*m]');
grid on
figure(7);subplot(4,2,7);plot(theta2*180/pi,Rg);
165
xlabel('Theta 2 [deg]');
ylabel('Rg [N]');
grid on
figure(7);subplot(4,2,8);plot(theta2*180/pi,Mg);
xlabel('Theta 2 [deg]');
ylabel('Mg [N]');
grid on
%%Título
set(gcf,'NextPlot','add');
axes;
h = title('Reacciones en el Árbol de Transmisión');
set(gca,'Visible','off');
set(h,'Visible','on');
166
Figura 2.83. Resultado de las reacciones y momentos en el brazo
izquierdo de la plataforma.
167
Figura 2.85. Resultados del torque estático necesario para la inclinación
de los misiles y las cámaras.
168
Figura 2.87. Resultados de las reacciones en la brida derecha.
169
Raz = 588.323 [N] Rby = 150.984 [N] Rcx = 0.000 [N]
Rb = 284.670 [N] Rbz = 241.331 [N] Rcy = 343.909 [N]
Rbx = 0.000 [N] Rc = 1005.525 [N] Rcz = 944.884 [N]
ª
¨Ë = [2.91]
«
2 ∙ e
¨H = [2.92]
²
170
Figura 2.90. Definición de la geometría del eje en Autodesk Inventor.
Del diagrama de cuerpo libre definido en [Link] se colocan las cargas externas
a las que se encuentra expuesto el eje a través de las herramientas localizadas
en “loads & supports” en la pestaña de cálculo. Como se muestra en la figura
2.91.
171
Por último, damos clic en “Calcular” para tener resultados sobre las reacciones
y fuerzas internas. En la pestaña de gráficos podemos observar todos los
resultados correspondientes a las fuerzas internas a generarse en el eje y
además el cálculo del diámetro ideal de nuestro eje. (Ver figura 2.92).
Figura 2.92. Resultados acerca del diámetro ideal del eje diseñado en
Autodesk Inventor.
9
9 8
32 ∙ ¹ h ∙ º 3 dº
" " "
> = ¸ ∙ 5- . # ∙/ 0 6 ¼
B ¹ 4 ¹»
Donde:
172
• FS: factor de seguridad, según el estándar MIL no puede ser menor a
5.
• Kf: factor de resistencia a la fatiga, que este caso se tomará de 2.
• Mmax: momento máximo del árbol, se calculó en [Nm].
• Tmax: torque máximo en el árbol, se calculó de [Nm].
• Sy: esfuerzo último de fluencia del material, en este caso 700[MPa]
correspondiente al (AISI 4340), material recomendado por árboles y
engranajes.
• Sn: resistencia a la fatiga, calculada de la siguiente manera:
¹»
¹ = - . ∙ m ∙ mÂN ∙ mà ∙ mÂ
2
Donde:
“En el diseño de flechas deben considerarse tanto los esfuerzos como las
deflexiones. La deflexión suele ser el factor crítico, ya que una deflexión
excesiva puede causar un desgaste rápido de los cojinetes de la flecha.15”
15
NORTON Robert L., “Diseño de Máquinas”, P217.
173
[Link] Programa de cálculo del diámetro mínimo por el método asme
b106, para el árbol de transmisión.
clc
%%Secciones del eje en mm.
x0=0/1000;
x1=10/1000;
x2=60/1000;
x3=112.5/1000;
x4=494/1000;
x5=876.5/1000;
x6=979/1000;
x7=989/1000;
%%Cálculo de fuerza cortante.
dx=x7/100;
x=[0:dx:x7];
vXZ(1:101)=0;
mXZ(1:101)=0;
vXY(1:101)=0;
mXY(1:101)=0;
i=1;
while x(i)<x1
vXZ(i)=0;
mXZ(i)=0;
vXY(i)=0;
mXY(i)=0;
i=i+1;
end
i1=i-1;
174
while x(i)<x2
vXZ(i)=vXZ(i1)-Rcz(101);
mXZ(i)=mXZ(i-1)+vXZ(i)*dx;
vXY(i)=vXY(i1)-Rcy(101);
mXY(i)=mXY(i-1)+vXY(i)*dx;
i=i+1;
end
i2=i-1;
while x(i)<x3
vXZ(i)=vXZ(i2)+Fr(101);
mXZ(i)=mXZ(i-1)+vXZ(i)*dx;
vXY(i)=vXY(i2)+Ft(101);
mXY(i)=mXY(i-1)+vXY(i)*dx;
i=i+1;
end
i3=i-1;
while x(i)<x4
vXZ(i)=vXZ(i3)+Rdz(101);
mXZ(i)=mXZ(i-1)+vXZ(i)*dx;
vXY(i)=vXY(i3)+Rdy(101);
mXY(i)=mXY(i-1)+vXY(i)*dx;
i=i+1;
end
i4=i-1;
while x(i)<x5
vXZ(i)=vXZ(i4)-Wej_cam.*cos(theta2(101));
mXZ(i)=mXZ(i-1)+vXZ(i)*dx;
vXY(i)=vXY(i4)-Wej_cam.*sin(theta2(101));
mXY(i)=mXY(i-1)+vXY(i)*dx;
i=i+1;
end
175
i5=i-1;
while x(i)<x6
vXZ(i)=vXZ(i5)+Rez(101);
mXZ(i)=mXZ(i-1)+vXZ(i)*dx;
vXY(i)=vXY(i5)+Rey(101);
mXY(i)=mXY(i-1)+vXY(i)*dx;
i=i+1;
end
i6=i-1;
while x(i)<x7
vXZ(i)=vXZ(i6)-Rcz(101);
mXZ(i)=mXZ(i-1)+vXZ(i)*dx;
vXY(i)=vXY(i6)-Rcy(101);
mXY(i)=mXY(i-1)+vXY(i)*dx;
i=i+1;
end
i7=i-1;
%%Gráficas de fuerza cortante y momento flector en el árbol de transmisión.
figure(8);subplot(2,2,1);plot(x*1000,vXZ);
xlabel('X [m]');
ylabel('V [N]');
title('Cortante V en XZ [N]');
grid on
figure(8);subplot(2,2,3);plot(x*1000,mXZ);
xlabel('X [m]');
ylabel('M [N*m]');
title('Momento M en XZ [N*m]');
grid on
figure(8);subplot(2,2,2);plot(x*1000,vXY);
xlabel('X [m]');
ylabel('V [N]');
176
title('Cortante V en XY [N]');
grid on
figure(8);subplot(2,2,4);plot(x*1000,mXY);
xlabel('X [m]');
ylabel('M [N*m]');
title('Momento M en XY [N*m]');
grid on
%%Cálculo de las cargas en la zona crítica
[Mxz xz]=max(mXZ);
[Mxy xy]=max(mXY);
if(Mxz>Mxy)
i=xz;
else
i=xy;
end
Mxz=mXZ(i);
Mxy=mXY(i);
Vxz=vXZ(i);
Vxy=vXY(i);
V=sqrt(Vxz^2+Vxy^2);
M=sqrt(Mxz^2+Mxy^2);
%%Diseño del diámetro mínimo método ASME B106
Sy=370*10^6; %Esfuerzo último de fluencia del AISI 1045.
FS=5; %Factor de seguridad estandar MIL.
Kf=2; %Factor de resistencia a la fatiga.
Sn=Sy/2; %Resistencia a la fatiga.
Cm=1; %Factor de forma material laminado en frío.
Cst=1; %Factor de esfuerzo.
Cr=0.9; %Factor de confiabilidad.
Cs=0.89; %Factor de tamaño Tabla...
Sn=Sn*Cm*Cst*Cr*Cs; %Límite de resistencia a la fatiga.
177
dmin=1000*((32*FS/pi)*sqrt((Kf/Sn)^2+(3/4)*(thau2(101)/Sy)^2))^(1/3);
178
Del cálculo de diámetro mínimo se obtuvo un valor de 32.2829mm, por lo tanto
se toman como correcto el diseño de Autodesk Inventor y nuestros cálculos.
Ahora se sabe que el diámetro del árbol debe superar los 32.2829mm.
Como se describe en 2.7.3, inventor cuenta con una poderosa herramienta para
diseño de engranajes. En la figura 2.95, se muestra el sistema de transmisión
empleado para llevar el movimiento de la parte baja, donde se encuentra el
servomotor, hacia la parte alta, donde se encuentran las cámaras.
179
Para los engranajes helicoidales se toman los mismo parámetros que en el
cálculo de esfuerzos de engranajes rectos, más el paso diametral. Además el
factor de forma J es diferente, para lo cual se utilizarán las tablas disponibles en
el libro de Robert L. Norton.
c£ ∙ }> l= ∙ l ∙ l ∙ l¤ ∙ lF
¨N = ∙ [2.93]
∙P l©
Donde:
• ¨N es el esfuerzo de tensión.
• c£ es la fuerza tangencial que se genera en el engrane debido al torque.
• }> es el paso diametral.
• l= es el factor de aplicación de la carga para resistencia a la flexión.
• l es el factor de tamaño para resistencia a la flexión.
• l es el factor de distribución de carga para resistencia a la flexión.
• l¤ es el factor de espesor de aro.
• lF es el factor de engrane intermedio o loco.
• es el ancho del diente.
• es el módulo del engrane.
• P es el factor geométrico para resistencia a la flexión.
• l© es el factor dinámico para resistencia a la flexión.
Los valores para cada engrane se obtuvieron empleando la tabla 2.21 y las
definiciones de 2.7.3. A continuación se presentan los factores y los esfuerzos
obtenidos. Además se calculó los factores de seguridad para el acero AISI 4340
con un límite de fluencia de 705 [MPa].
180
Tabla 2.21. Cálculo de esfuerzo para los engranajes helicoidales.
Engrane Helicoidal
Parámetro
1 2 3
Z 58 80 114
m 3 3 3
Ì 20 20 20
Í 20 20 20
Fuente: Propia.
Se verifica que los factores de seguridad coinciden con aquellos calculados por
Autodesk Inventor.
181
Figura 2.96. Resultados de los factores de diseño de los engranajes
helicoidales.
182
Figura 2.97. Cargas sobre la bancada.
ª ª
¨Ë = = [2.94]
« 2∙ ∙∙>
Donde:
183
El esfuerzo máximo del material a emplear en la bancada, Acero AISI 1045 es
de 370[MPa], obtenido del catálogo referencial de Tata Steel. Se colocarán 6
pernos, pero se realizará el diseño para el caso en el que todos estén flojos,
menos uno. Se toma en cuenta el factor de seguridad mínimo para aplicaciones
militares, FS=5.
ª 2572.600[]
> = ¹ =5∙ = 4.346[]
2∙ ∙ ∙ ¨M 2 ∙ 1 ∙ 4[] ∙ 370[°=]
2 ∙ "
2 ∙ e
¨H = =
² ∙ > ∙ 8
"
Donde:
• I es la inercia de la sección.
• c es la distancia al eje neutro.
• d es el espesor del material.
• l es la longitud de la sección.
Se puede observar un mayor efecto del cortante que de los momentos sobre la
bancada, por lo tanto se selecciona una plancha metálica de acero AISI 1045
de 5[mm] de espesor.
184
El análisis de deflexiones para la geometría propuesta resulta complejo, por lo
cual se ha realizado una simulación en Autodesk Inventor para simular el
comportamiento de la bancada frente a sus cargas, en 2.9.2 se presentan
resultados.
Los rodamientos de bolas a contacto angular son buenos para soportar cargas
combinadas en ambos sentidos de giro, como es el caso al utilizar engranajes
helicoidales. Además tienen una rigidez elevada y buena precisión en
funcionamiento.
16
<[Link]
[Link]&file=1_0_2&maincatalogue=1&lang=en>
185
Una vez seleccionado el tipo de rodamiento, los parámetros de selección más
importantes son: las dimensiones, dados eje y alojamiento, las cargas radial,
axial y la velocidad de operación. En base a la tabla 2.22 se seleccionaron los
rodamientos respectivos.
Velocidad [RPM] 10 15
Fuente: Propia.
Figura 2.99. Cálculo del factor de carga del rodamiento del árbol de
transmisión.
186
Figura 2.100. Cálculo del factor de carga del rodamiento del eje loco.
187
de fácil montaje que disminuyen los tiempos de ensamblaje y mantenimiento de
las máquinas que los utilizan.
El árbol de transmisión debe poseer un agujero pasante para llevar los cables
del disparador de los misiles, por lo tanto su diámetro deberá ser mucho mayor.
El fabricante de FennerDrives recomienda que el agujero no tenga un diámetro
mayor al 35% del diámetro total del eje, debido a las radiales ejercidas sobre el
eje. Se considera un diámetro del agujero de 25[mm] para pasar el cableado del
disparador, con lo que resulta un diámetro mínimo del árbol de 72.429[mm], por
este motivo se seleccionó un rodamiento mayor en el ítem 2.8.7, de 80[mm] de
diámetro interno.
Fuente: Propia.
188
Figura 2.101. Menú de la herramienta Stress Analysis.
Una vez ubicadas todas las cargas, se procede a realizar la simulación dando
click sobre el ícono Simulate del menú de Stress Analysis, como se observa en
la figura 2.103. Se siguen todos los pasos y se generan los resultados.
189
Figura 2.103. Ícono Simulate.
De estos resultados los que tienen mayor importancia para esta aplicación son
el esfuerzo de Von Mises, el de desplazamiento y el de factor de seguridad.
190
Tabla 2.24. Resultados de las simulaciones por análisis de esfuerzos en
los elementos de la plataforma.
Fuente: Propia.
191
Los planos se realizaron en base a las normas DIN para planos de dibujo
técnico.
Los formatos de la serie DIN A (normalizados por las normas DIN 823 y la DIN
476) se pueden subdividir racionalmente así: AO en dos formatos A1; en cuatro
formatos A2; en ocho formatos A3; en dieciséis formatos A4. Esta subdivisión
se identifica como doblez modular, como se muestra en la figura 2.104.
192
Para el desarrollo de los planos se ha utilizado el formato DIN A3, el cuál posee
las siguientes medidas: 420 mm de ancho y 297 mm de alto.
1) Línea gruesa de trazo seguido: Se emplea para todas las aristas visibles.
2) Línea fina de trazo seguido: Este tipo de líneas se emplean en líneas de
cota, líneas auxiliares de cota, diámetros de núcleos de roscas en pernos y
diámetros nominales de roscas en turcas. Además, se emplea en secciones
rebatidas, rayadas, líneas de referencia, cruces de diagonales y líneas de
dobladura. Líneas de trazos: Se emplean para aristas ocultas (no vistas) y
para contornos. Línea de trazo y punto (gruesa): Sirve para indicar el
193
recorrido de las secciones y los planos de sección así como para indicar los
tratamientos superficiales. La línea comienza y acaba con trazo.
3) Línea de trazo y punto (fina): Se emplean para ejes, circunferencias
primitivas, en ruedas dentadas, representación de formas iniciales.
0,18-0,25-0,35-0,5-0,7-1-1,4 y 2 [mm]
194
La relación entre las anchuras de las líneas finas y gruesas en un mismo dibujo,
no deben ser inferior a 2. Deben conservarse la misma anchura de línea para
las diferentes vistas de una pieza, dibujadas con la misma escala.
En la figura 2.107 se dan 6 tipos de líneas, las cuales se indican con un número
sobre ellas que representa su anchura en décimas de milímetros.
Con el fin de alcanzar la armonía del dibujo, se dan cuatro grupos de líneas
como se muestra en la figura 2.108.
195
Para representar las superficies cortadas, se utiliza la norma DIN 201, en donde
el rayado se debe hacer a 45° con las aristas del contorno, con una separación
de 0.7 a 3 [mm] según el tamaño del modelo y con un espesor de 1/3 de las
líneas de las aristas vistas del contorno. Como se muestra en la figura 2.109.
196
Considerando todas estas normas y en base a las bases técnicas de las piezas
descritas en el índice 2.10.1 se desarrollaron los planos en el software Autodesk
Inventor, y se los puede observar en el apéndice C.2.
Para el plegado de los planos se empleó la norma DIN 824, la cual especifica
cómo se deben doblar los planos para cada formato. En este caso el doblado
del formato DIN A3 se muestra en la figura 2.110.
197
CAPITULO 3
AUTOMATIZACIÓN DE LA PLATAFORMA
Para llevar a cabo la tarea de seguimiento el sistema debe contar con una
pantalla donde el operario pueda visualizar a la aeronave y la mira colimada de
la cámara con los misiles (Lugar donde apuntan los misiles). Y para llevar a
cabo las tareas de posicionamiento manual de la plataforma y enganche de la
aeronave, se cuenta con un joystick militar.
El dato del error entregado por el sistema de visión artificial será nuestra
entrada para llevar a cabo los movimientos necesarios en orientación y
elevación que localicen a la aeronave siempre en el centro de la pantalla, que
debe estar colimado con las miras de los misiles que se encontrarán sobre la
misma.
198
En base a lo anterior, podemos plantear como primer requerimiento que los
movimientos de la plataforma no deben permitir que la aeronave abandone el
ángulo sólido de visibilidad de los misiles luego de haber sido enganchada.
Para esto se requerirá de un sistema de visión artificial que entregue muestras
del error en un tiempo menor o igual al tiempo de muestreo mínimo que se
calcula en el ítem siguiente, y de un controlador capaz de comunicarse con los
servo drives para cambiar su velocidad antes de que la siguiente muestra
procedente del sistema de visión llegue.
Por otra parte, se debe garantizar que no se excedan los rangos de movimiento
en orientación de 330 grados, ni en elevación entre 20 y 70 grados. Para esto
será necesario el uso de encoder absolutos que guarden la posición de la
plataforma aún después de haber quedado sin energía y finales de carrera para
el caso de que se llegase a las posiciones extremas.
199
Figura 3.1. Diagrama de bloques del sistema.
El punto suma representa los módulos responsables de adquirir las señales del
sistema de visión artificial, de los encoders absolutos, el joystick y los
dispositivos de emergencia para posteriormente procesarlas en el controlador.
Los servo drives llevan a cabo la tarea de mover los servomotores de acuerdo a
las instrucciones enviadas por el controlador.
200
Luego de llevarse a cabo cada movimiento, la entrada de la cámara térmica, los
encoders y dispositivos de emergencia cambian y se inicia nuevamente el lazo
de control.
Partiremos del cálculo del tiempo que tarda una aeronave a la velocidad
máxima de seguimiento de 360[m/s] en salir del campo de visión de una cámara
convencional de acuerdo a la figura 3.2.
Ê >
tan - . 3.1
2 2∙
> © ∙ ∆£ 3.2
201
2∙ Ê 2 ∙ 1000 50°
∆£ ∙ tan - . ∙ tan - . = 2.591 3.3
© 2
360 ¾ À 2
Â
2.591
d = = 51.820
50
El encoder del servomotor le informa al servo drive sobre el estado real del
servomotor, y de esta manera el servo drive se encarga de realizar el control en
base a la orden de mando.
202
en base al catálogo de configuración del sistema para los servomotores
SureServo. En el cuál se especifica que servo drive es el ideal para cada
servomotor. Existen servo drives que son compatibles para diferentes marcas
de servomotores, pero es recomendable utilizar el servo drive de la misma
marca que el servomotor.
Para el servomotor SVM 210B el catálogo muestra que el servo drive para este
servomotor es el SVA 2100, y para el servomotor SVM 220B el servo drive es el
SVA 2300. Como se muestra en la tabla 3.1.
Estos servo drives son AC; pueden ser de simple fase como es el caso del SVA
2100 y trifásico como es el caso del SVA 2300; poseen un módulo de
comunicación Modbus RTU, permiten el ingreso de datos de forma externa ya
sea por computador o un controlador, o de forma interna ya sea ingresando los
203
parámetros en el display. Las tablas con todas las especificaciones de los servo
drives se muestran en los apéndices A.3 y A.4.
17
Ibíd., P. e16-25, Vol. 13
204
Figura 3.4. Conexión del servo drive SVA 2100 / SVA 230018.
Los cables para conectar el servomotor al servo drive (cable de potencia del
servomotor) y el cable del encoder del servomotor se adquieren en base a la
tabla 3.1, y se los muestra en la figura 3.5.
18
Ibíd., P. e16-35, Vol. 13
205
Figura 3.5. Cable de potencia del motor y cable del encoder19.
19
Ibíd., P. e16-25, Vol. 13
206
el año de [Link] un marco de referencia para la definición de arquitecturas de
interconexión de sistemas de comunicaciones.
1) Es público.
2) Su implementación es fácil y requiere poco desarrollo.
3) Maneja bloques de datos sin suponer restricciones.
20
< [Link]
207
Modbus permite el control de una red de dispositivos y además se usa para la
conexión de un ordenador de supervisión con una unidad remota RTU. Un
dispositivo de unidad remota RTU es un dispositivo basado en
microprocesadores, el cual permite obtener señales independientes de los
procesos y enviar la información a un sitio remoto donde se procese. La
implementación de este protocolo con estos dispositivos RTU, dio origen al
protocolo Modbus RTU, el cuál puede ser implementado sobre una red RS-485
o RS-232.
RS232 es una interfaz que designa una norma para el intercambio de una serie
de datos binarios entre un equipo terminal de datos (DTE) y un equipo de
comunicación de datos (DCE), para distancias cortas, de hasta 15 metros
según la norma , y para velocidades de comunicación bajas, de no más de
20 kbps. A pesar de esto, muchas veces se utiliza a mayores velocidades con
un resultado aceptable. La interfaz puede trabajar en comunicación
asíncrona o síncrona y tipos de canal simplex, half duplex o full dúplex, estos
últimos simplemente indican la dirección de cómo viajan los datos.
208
puerto del PLC se controle los dos servodrives, y estos estén colocados a
largas distancias, se utilizaría la red de comunicación RS-485.
Como se trata de un protocolo que trabaja sobre una red serial RS-232, es
necesario conocer el valor de los parámetros para la comunicación, y estos son:
la velocidad de transmisión, los bits de datos, el bit de paridad y los bits de
parada.
Es así como para poder tener la comunicación es necesario que todos los
dispositivos que envíen y recepten información entre sí, se comuniquen con los
mismos valores de los parámetros expuestos anteriormente.
209
El resto de dispositivos conectados al PLC, como los servo drives trabajarán
como esclavos.
[Link] Antecedentes.
210
En base a estos parámetros se seleccionó los encoders de Allen Bradley de la
serie 842A, los cuáles son encoders absolutos magnéticos multivuelta de 25
bits de dos direcciones, que poseen una precisión de 8192 pulsos por
revolución, el formato de código es gray o binario y su interfaz de comunicación
es SSI. En la figura 3.8 se muestra este tipo de encoders.
21
Catálogo “Encoders Allen Bradley”, P. 6-2
211
Figura 3.9. Convertidor Motrona IV25122.
• Sensor capacitivo
• Sensor inductivo
• Sensor fin de carrera
• Sensor infrarrojo
• Sensor ultrasónico
• Sensor magnético
22
Manual “Motrona IV251”, P. 1
212
El sensor fin de carrera es un dispositivo que puede ser eléctrico, neumático o
mecánico, y es situado por lo general al final del recorrido de un elemento móvil,
con el objetivo de enviar señales que puedan modificar el estado de un circuito.
Internamente puede contener interruptores normalmente abiertos (NA),
cerrados (NC) o conmutadores dependiendo de la operación que cumpla al ser
accionado.
Las ventajas que este sensor muestra para esta aplicación con relación al resto
de sensores son:
Existen diversas marcas de finales de carrera, entre las que se puede nombrar
Allen Bradley y OMRON, que han sido las marcas analizadas para seleccionar
el final de carrera que mejor se adapte para esta aplicación.
213
Figura 3.10. Sensor de final de carrera OMRON D4B5A71-N23.
Este sensor trabaja hasta con 250 [VAC] y 60 [A], valor ideal ya que los
servomotores trabajan con 220 [VAC], para su conexión es necesario un fusible
de 10 [A]. Los cables de conexión pueden ser de tamaño de 20 a 14 AWG, es
necesario utilizar terminales de crimpar redondos M3.5 como se observa en la
figura 3.11.
23
Catálogo “OMRON Finales de Carrera”, P. G-249
24
Catálogo “OMRON Finales de Carrera”, P. G-263
214
La conexión es en serie con la alimentación de los servomotores. Este sensor
posee dos contactos internos cerrados, y al ser pulsado el émbolo superior
estos se abren. En la tabla 3.2 se puede ver la configuración de los contactos.
Este sensor tiene una larga vida útil con más de 30.000.000 operaciones y es
un sensor bastante rápido, que permite accionar inmediatamente los contactos
apenas el émbolo es presionado.
[Link] Antecedentes.
215
• Ambiente de operación, el controlador seleccionado debe poseer un rango
de trabajo dentro de las condiciones ambientales en las cuales
desempeñará las acciones de control.
• Dispositivos digitales, se debe determinar el número de dispositivos digitales
que requiere la aplicación.
• Dispositivos analógicos, de igual manera es necesario conocer el número y
la resolución de los dispositivos analógicos que formarán parte de la
aplicación.
• Características especiales, de existir alguna característica especial de la
aplicación, debe tomarse en cuenta en la selección; como por ejemplo el
entorno de programación a seleccionar, contar con contadores de alta
velocidad, etc.
• Requerimientos del CPU, entre los requerimientos fundamentales tenemos
el tamaño de la memoria, la velocidad de procesamiento y las capacidades
de cálculo con las que cuenta.
• Ubicación de E/S, se debe conocer la ubicación de las entradas y salidas, si
serán remotas o estarán en el mismo chasis del CPU.
• Comunicaciones, los protocolos a implementar en la aplicación.
216
Tabla 3.3. Requerimientos del controlador.
Fuente: Propia.
217
Figura 3.12. PLC Direct Logic 205 Koyo25.
Fuente: Propia.
218
En la figura 3.13 se muestra un ejemplo de módulo del PLC.
El PLC DirectLogic 205 cuenta con soporte para Windows sin la necesidad de
un programador adicional. Su programación se realiza vía Ethernet a través del
entorno de programación Think&Do.
Encontramos los tipos de dato siguientes: Array, ASCII, Axis, Comm, Counter,
Flag, Float, Input, Number, Output, Register, String, System y Timer. Podemos
visualizar fácilmente las variables de cualquier tipo de dato en la tabla de
variables del programa. Seleccionamos el tipo de dato y se observa la lista de
variables con su tagname, es importante notar que no se puede desplazar los
tagnames borrando las variables, estos quedan asignados en base al ID dentro
de la tabla.
26
<[Link]
ic_Series_PLCs_(Micro_to_Small,_Brick_-a-
_Modular)/DirectLogic_205_(Micro_Modular_PLC)/Motion_-a-_Specialty_Modules/H2-SERIO>
219
Figura 3.14. Tabla de variables de Think&Do.
220
Figura 3.15. Pantalla de configuración del PLC en Think&Do.
222
Figura 3.18. Flowcharts para la programación.
223
flechas hacia los bloques siguientes. Los bloques de bifurcación son rombos
verdes y los bloques que ejecutan operaciones son cuadrados, el color
dependerá de la operación que realicen.
224
El bloque de comparación devuelve verdadero o falso frente a una comparación
numérica de mayor que, menor que o igual a, entre variables internas.
225
El bloque de cálculo nos permite asignar valores correspondientes a
operaciones matemáticas con cualquier tipo de dato.
226
El bloque de llamada nos permite saltar hacia una subrutina estándar o
programada por el usuario, luego de ejecutada la subrutina el programa
continuará con los bloques siguientes. Se puede utilizar bloques de llamada en
las rutinas y subrutinas de la misma manera.
Se observan algunas de las subrutinas a las cuales se pueden llamar dentro del
programa en cuestión.
227
Figura 3.27. Bloque de comentarios con su ventana de edición.
[Link] Antecedentes.
228
operaciones de entrada y salida de datos. Para esta aplicación es necesario
una pantalla táctil que nos permita setear y monitorear el comportamiento de la
plataforma, y un joystick que nos permitirá mover a la plataforma en orientación
y en elevación en modo manual, y permita setear mediante un botón el modo
automático del movimiento.
Para asegurar que exista un correcto funcionamiento del sistema se realizó una
búsqueda para ver si existe una pantalla táctil de la misma marca que el PLC.
229
Figura 3.29. Pantalla táctil EA7-T12C de Automation Direct27.
27
Catálogo “C-MORE HMI AUTOMATION DIRECT”, P. e10-31, Vol. 14
28
Catálogo “C-MORE HMI AUTOMATION DIRECT”, P. e10-41, Vol. 14
230
[Link] Selección del joystick.
Lo ideal para la aplicación es contar con una fuente de poder que entregue la
carga necesaria a todos los componentes de control como lo son: el PLC, la
pantalla HMI, el joystick y los encoders. Por lo tanto se requiere conocer acerca
del consumo de estos elementos para tomar la mejor decisión, ya que una
fuente mal dimensionada no proporcionará la energía necesaria al sistema de
231
control para que funcione correctamente o puede que su vida útil se vea
reducida.
Joystick 12 100
Fuente: Propia.
232
Seleccionamos por estar dentro del rango requerido a la fuente de poder PSC-
24-090. Esta fuente al ser de la misma procedencia del fabricante garantiza un
buen acoplamiento de impedancias y además notamos que entrega 24 [V] y
hasta 3.75 [A], por lo cual se considera que está apta para suplir la energía
necesaria al sistema de control. Ver figura 3.32.
29
Catálogo “ POWER SUPPLIES AUTOMATION DIRECT”, P. e30-33, Vol. 14
30
Catálogo “ POWER SUPPLIES AUTOMATION DIRECT”, P. e30-33, Vol. 14
233
la eficiencia energética de una fuente adicional es mejor que la de un regulador,
que en bajar de 24 a 12 [V] disipa una buena parte de la potencia en calor.
234
Los parámetros indispensables a configurar mediante el archivo
“nombre_del_archivo.ssc” o manualmente son:
235
5) Entradas Digitales (Digital Input), el servodrive cuenta con 8 registros que
pueden funcionar a manera de entradas digitales. Para la operación correcta
de la plataforma las entradas digitales deben configurarse de la siguiente
manera:
Nota: Para más detalles de los parámetro del servo drive referirse al Apéndice
A.20.
236
Figura 3.35. Algoritmo de escritura.
[Link] Subrutina de lectura del servo drive.
237
Figura 3.36. Algoritmo de lectura.
[Link] Rutina principal.
238
Figura 3.37. Algoritmo de la rutina principal.
Cada vez que el programa necesita comunicarse con el servo drive para
configurar algún registro, recurre a las subrutinas de escritura y lectura antes
descritas.
239
[Link] Rutina de movimiento del servomotor.
240
acuerdo al sentido de giro respectivo, mientras que en los estados 2 y 4 la
velocidad de salida trabaja de manera proporcional al voltaje de entrada. En
cuestión al estado 3 la velocidad de salida es 0. La presencia de los 5 estados
de funcionamiento se da debido a que el valor analógico obtenido del joystick
varía aun sin mover el joystick, lo cual generaría variaciones indeseables en la
velocidad de movimiento de la plataforma. Se realizó la lectura de los valores
analógicos mínimo y máximo para el rango en donde el estado 3 será
considerado, obteniendo que la velocidad será 0 entre los valores 2100 y 2200.
ª = 2«±m=! 3.4
ª 8 = 0 3.6
ª = «±m=! 3.8
241
[Link] Lista de registros del servo drive utilizados.
Fuente: Propia.
242
3.10 DESARROLLO DE LA INTERFAZ HUMANO-MÁQUINA
(HMI).
243
que la pantalla táctil trabajará. Es importante tener el conocimiento acerca de
que variables se van a controlar y que variables se van a monitorear.
También se configuran algunos parámetros más que para esta aplicación se los
puede dejar con el valor que vienen por default. En las figuras 3.41 y 3.42 se
muestra la ventana del Panel Manager para la configuración del tipo de
pantalla y para la configuración de la comunicación de la pantalla con el PLC.
Figura 3.41. Ventana del Panel Manager para la configuración del tipo de
pantalla.
244
Figura 3.42. Ventana del Panel Manager para la configuración de la
comunicación de la pantalla con el PLC.
245
Figura 3.43. Ventana del Panel Manager mostrando el buscador del
archivo del PLC.
En la zona derecha del programa existen dos listas, la lista de objetos y la lista
de partes. En la lista de partes y de objetos, es en donde se selecciona el objeto
a colocar en la pantalla, aquí se encuentran botones, displays, switchs,
cualquier tipo de imagen que se quiera colocar en la pantalla, etc. En estas
listas encontramos todos los objetos necesarios para programar la pantalla. La
figura 3.44 muestra las listas de partes y de objetos.
246
Para colocar los objetos en la pantalla simplemente es necesario arrastrarlos
hacia la pantalla.
Una vez colocados los elementos necesarios para esta aplicación, se tiene una
pantalla inicial del programa como se muestra en la figura 3.45.
247
de operación, y en esta se puede monitorear el movimiento del helicóptero y el
seguimiento automático. Para esta aplicación está desarrollada una simulación
del movimiento del helicóptero, la cual permite comprobar el correcto
funcionamiento del modo automático de la plataforma. En la figura 3.46 se
muestra la pantalla de modo automático en la versión de simulación.
248
Figura 3.47. Pantalla de configuración de los objetos de la pantalla.
Una vez configurados todos los objetos de las pantallas, es necesario descargar
el programa a la pantalla táctil, para esto se debe conectar la pantalla a la
computadora mediante un cable Ethernet. Dentro del programa de C-MORE es
necesario dar click sobre la pestaña “SEND PROJECT TO PANEL”, de esta
manera se abre una nueva ventana en donde nos pide colocar la dirección IP
de la pantalla táctil, la cuál es [Link]; la dirección viene especificada en la
pantalla táctil. Una vez comprobada la conexión, es posible enviar el programa
desarrollado a la pantalla táctil, dando click en el botón “TRANSFER”. La
ventana de configuración de la dirección IP se muestra en la figura 3.48.
249
Figura 3.48. Ventana de configuración de la dirección IP para la
comunicación de la pantalla táctil con la computadora.
31
Catálogo “C-MORE HMI AUTOMATION DIRECT”, P. e10-42, Vol. 14
250
3.11 SIMULACIÓN DEL SISTEMA DE CONTROL SIN CARGA.
251
Figura 3.50. Algoritmo de simulación del controlador.
252
Figura 3.51. Algoritmo del botón de enganche encendido y apagado.
253
CAPITULO 4
PRUEBAS DE FUNCIONAMIENTO
254
Inicialmente el PLC no se comunicaba con la computadora, se probó varias
hipótesis de error, entre las cuáles se probó como primera fuente de error si el
problema era el cable Ethernet, el mismo que al conectarlo a internet funcionó
sin ningún inconveniente. Por tal motivo se buscó nuevas hipótesis del error.
255
Figura 4.3. Intento de comunicación con el PLC mediante la ventada de
comandos de Windows.
256
Figura 4.5. Ventana de comandos de Windows, haciendo ping con el PLC.
Una vez lograda la comunicación con el PLC, la ventana del paquete ESP de
Think&Do se muestra en la figura 4.6 con la comunicación establecida.
257
Figura 4.7. Comunicación con el PLC en Think&Do.
El programa que permite acceder a los registros del servo drive desde el
computador es el SureServo Pro, software que se lo puede descargar desde
internet de la página de Automation Direct.
258
Figura 4.8. Programa SureServo Pro.
259
Y es este el problema que se generó inicialmente, por el desconocimiento de la
configuración inicial de comunicación del servo drive, se intentó realizar la
comunicación con el computador y esta no se daba. Se probó los cables, y
estos no eran el error, se probó la entrada USB del computador, y esta
funcionaba sin ningún inconveniente. Hasta que finalmente leyendo el manual
de comunicación del servo drive, se pudo conocer cuál era la configuración
inicial para la comunicación del servo drive, y como configurarla de diferente
manera. De esta forma configurando el computador con los mismos valores de
los parámetros de comunicación del servo drive, se logró la comunicación y se
pudo leer y escribir los registros del servo drive desde el computador.
Comprobando así que la comunicación del servo drive funciona sin ningún
inconveniente y de que es posible controlar y monitorear el servo drive desde
controladores externos.
260
Como era conocida la dirección IP de la pantalla, y además su dirección se
encuentra establecida dentro del dominio de comunicación del PLC y del
computador, esta se comunicó sin problemas, comprobado esto, se realizó una
prueba de comunicación dentro del software C-MORE, la cual se logró con éxito
sin ningún inconveniente (ver figura 4.11) y se procedió a transferir un programa
básico en C-MORE a la pantalla (explicación capítulo 3.7), para comprobar que
esta funciona sin problemas. El programa constaba de un botón pulsador y de
un led. Y sin ningún problema los elementos aparecieron en la pantalla táctil.
Posteriormente se realizó el programa final que servirá para esta aplicación y se
lo cargó en la pantalla táctil.
261
Figura 4.12. Pantalla táctil con el programa desarrollado.
262
DB9 en ambos terminales. De manera semejante a la prueba anterior, se utilizó
un cable conversor de DB9 a USB para la entrada en el computador.
Para poder hallar una solución al problema, se leyó varias veces el manual de
configuración del Motrona, y se buscó en internet algún documento que pueda
de ser ayuda para este caso, y es aquí en donde se encontró un documento de
comunicación del Motrona, en el cuál se especifica que el cable de
comunicación del Motrona debe tener la señal invertida, esto es que comunica
la señal de recepción de entrada con la de transmisión de salida y viceversa. De
esta manera se procedió a realizar este tipo de cable (ver apéndice E.9 figura
E.9.3). Con el cable realizado se procedió a realizar la prueba de comunicación,
y finalmente se logró leer el encoder desde el computador.
264
Durante el desarrollo del programa que realmente irá en esta aplicación, se
presentaron varias complicaciones, principalmente cuando se configuró mal el
programa del PLC y se intento escribir y leer al mismo tiempo registros del
servo drive. Cuando se generó esta situación el PLC se colgaba, y el servo
drive entraba en modo de error “ALE 19”, el cuál es un error de comunicación.
Para solucionar este problema fue necesario comprender que el PLC trabaja en
serie y no en paralelo, por tal motivo no es posible escribir y leer registros al
mismo tiempo, entonces fue importante modificar el programa del PLC con
banderas, para evitar que varias funciones corran al mismo tiempo, de esta
manera se controlo que el PLC se encuentre o solo leyendo o solo escribiendo
registros en el servo drive. Y de esta manera el programa funcionó
correctamente.
Una vez cargado el programa en la pantalla táctil, para poder controlar el PLC
desde la pantalla táctil es necesario comunicarla al PLC mediante un cable
serial que posea un terminal DB15 y el otro terminal RJ12. Anterior a esto se
debe configurar el puerto serial del PLC como se explica en el capítulo 3.4.2.
Una vez conectado se encienden los dispositivos y se procedió a probar que la
pantalla se comunique correctamente con el PLC.
265
4.1.7. PRUEBAS DE COMUNICACIÓN ENCODER - PLC.
Para comunicar el encoder con el PLC se utilizó un cable serial con un terminal
DB9 y el otro extremo RJ12. Lo primordial para lograr leer los datos del encoder
en el PLC, fue desarrollar correctamente el programa de lectura en el PLC.
Desarrollado el programa se obtuvo la lectura de datos sin ningún
inconveniente.
266
envían el valor de velocidad con el que se va a mover cada servomotor.
Posteriormente se probó lo mismo con la diferencia de que el valor de la
velocidad que se envía a los servomotores será configurada en base a la
posición del joystick. Para lo cual se tiene un programa con control proporcional
dentro del PLC, que en base al valor analógico recibido por el movimiento del
joystick, lo transforma en un valor de velocidad y lo envía a los servo drives.
267
De manera similar los sensores de finales de carrera trabajan como paros de
emergencia y cumplen la misma acción que el botón de paro de emergencia.
268
Se llevó a cabo el montaje del servomotor SVM-210 en una base de fuego de
23mm con la que cuenta el CICTE a través de un acoplamiento provisional.
269
Una vez corregida la localización de los centros se manufacturó el acoplamiento
en acero de transmisión, ya que para la etapa de pruebas este material
proporciona una resistencia suficiente a un costo bajo.
Como primera prueba fue positivo observar que la plataforma se elevó, pero era
necesario conseguir un movimiento entre 0 y 90 grados tanto de subida como
bajada, por tal motivo fue necesario diseñar un segundo sistema de
acoplamiento con el cual se procuró aumentar la relación de transmisión, para
entregar mayor torque a la cuna.
Este acoplamiento fue empleado para llevar a cabo las diferentes pruebas que
se detallan en los ítems siguientes.
271
Los resultados obtenidos fueron positivos ya que la nueva relación de
transmisión permitió mover al sistema de elevación de la plataforma en el rango
completo de la cuna.
Se llevaron a cabo las pruebas descritas en 4.2 mas esta vez con el servomotor
acoplado a la plataforma. Mediante el joystick se variaba la velocidad de
operación de la cuna, comprobando un correcto funcionamiento del sistema de
control en todos los aspectos: el joystick pone a la velocidad máxima al
servomotor en las posiciones extremas, es proporcional en las posiciones
intermedias y en la posición central no se manifestó ningún movimiento.
272
4.4.4. PRUEBA DE MOVIMIENTO EN MODO AUTOMÁTICO.
273
Para estas pruebas se limitó el movimiento de la plataforma para que trabaje de
20 a 70 grados que es el rango de funcionamiento de los misiles Igla, la
medición de la posición se obtiene de los monitores del servo drive, que
almacena las lecturas del encoder.
90
| = H ∙
1690
Se realizaron pruebas para los valores de velocidad siguientes: 1, 5, 20, 40, 60,
80, 100, 120 [RPM], en un rango entre 20° y 70°.
274
[Link]. Descripción del código en Labview.
275
2) Se escriben los registros de aceleración, desaceleración y curva S para
garantizar que el servomotor operará de acuerdo a la misma configuración
cada vez y se habilita la EEPROM como se hizo en el programa del PLC
Koyo.
276
4) Habilitar el servodrive.
277
Figura 4.25. Deshabilitación del servomotor y cierre del puerto de
comunicaciones.
Se observan los resultados en las 6 pantallas con su respectivo título. Los valores de
los resultados que se presentan en el lado derecho serán proporcionales a los de la
izquierda en 18.8:1, que es la relación de transmisión empleada para acoplar el
servomotor a la base de fuego.
278
[Link]. Resultados de la prueba con Labview.
a) RESULTADOS A 1 [RPM]
279
debido a que la carga crece conforme mayor es su inclinación, el servomotor se
sobrecarga al pasar los 60° y se deshabilita.
b) RESULTADOS A 5 [RPM]
280
c) RESULTADOS A 20 [RPM]
281
d) RESULTADOS A 40 [RPM]
282
e) RESULTADOS A 60 [RPM]
283
f) RESULTADOS A 80 [RPM]
284
g) RESULTADOS A 100 [RPM]
285
h) RESULTADOS A 120 [RPM]
286
[Link]. Resumen de resultados.
Se han tomado las gráficas en una subida y una bajada de la cuna para cada
una de las velocidades del servomotor sujetas a prueba. Se consideran las
gráficas de posición para evaluar los resultados debido a que los resultados de
velocidad y torque en el tiempo resultan irrelevantes luego de conocer sus
valores máximos. Se observa que el resultado más fino con el servomotor a 1
[RPM] y a 120 [RPM], sin embargo no se puede confiar de los movimientos en
rango intermedios. Esto se debe a la falta de torque entregado a la cuna,
debiendo incrementarse el torque a través de una caja de transmisión con una
relación bastante superior a 18.8:1 o cambiando de servomotor a uno con
287
mayor torque, sin embargo esto no sería lo óptimo debido a que un servomotor
con mayor torque ocupa mayor espacio y consume una mayor potencia
eléctrica.
Con los datos obtenidos de las pruebas se puede calcular la inercia promedio
de la plataforma, dividiendo el torque transmitido para la aceleración durante el
movimiento de la cuna.
O"8
P"8 = 4.1
<"8
Para llevar a cabo este cálculo se requiere despreciar la carga estática que
aparece producto de los cañones a diferente altura. Realizando la división punto
a punto para cada una de las pruebas realizadas se obtienen las siguientes
inercias máximas:
PRUEBA ÔÔÔÔÔÔÔÔ
ÑuÒÓÓ
1 [RPM] 108.24
5 [RPM] 34.72
20 [RPM] 34.46
40 [RPM] 34.74
60 [RPM] 34.58
80 [RPM] 34.58
ÔÔÔÔÔÔÔÔ
ÑuÒÓÓ 46.17
Fuente: Propia.
288
Se observa que la inercia máxima calculada en las pruebas realizadas fue de
46.17 [kg∙m2], con lo cual resulta la relación de transmisión óptima siguiente:
Debido al limitado espacio que posee la plataforma para colocar una caja de
transmisión resulta complejo implementar un reductor con i=277, ya que no se
puede utilizar un reductor de tornillo sin fin debido a la necesidad de conservar
un sistema de respaldo para movimiento manual de la base de fuego, es
necesario bajar al máximo posible esta relación. Además una relación de 277:1
requiere que el servomotor trabaje a velocidades muy elevadas.
289
Notamos que la relación mínima para mover la plataforma es de
aproximadamente 0.05 veces la relación óptima 13.89:1, al haber utilizado una
relación de 18.8:1 en las pruebas se logró el movimiento de la cuna, más no se
obtuvieron los resultados esperados.
Además la velocidad máxima del servomotor para alcanzar las 10 [RPM] con la
relación óptima de 277:1 es de 2770 [RPM], valor que excede el rango de
trabajo del servomotor SVM-210 a emplear.
290
thauM1=2*sqrt(JM*(JL1*alphaL1^2+thauL1*alphaL1));
291
%%Cálculo de la energía cinética requerida por el servomotor.
KM1k=(1/2)*k.^2*(JL1+(thauL1/alphaL1))*ohmegaL1^2;
i1=iOPT1.*k1
thau1=thauL1MIN./i1
ohmega1=(ohmegaL1*0.5*60/pi).*i1
W1=thau1.*ohmega1*2*pi/60
292
CAPITULO 5
CONCLUSIONES Y RECOMENDACIONES
5.1 CONCLUSIONES.
293
4) Se seleccionaron los componentes electrónicos para la automatización
de la plataforma en base a la disponibilidad del CICTE.
294
10) El servomotor SVM-210 resultó bastante confiable tanto en su robustez
como en su protocolo de comunicaciones, no obstante para utilizarlo en
la automatización de la base de fuego ZU-23 es necesario el diseño de
una nueva caja de transmisión con una relación de transmisión entre
138.5 y 200, resultado obtenido de la prueba descrita en 4.3.6.
5.2 RECOMENDACIONES.
295
REFERENCIAS
BIBLIOGRAFÍA.
LINKOGRAFÍA.
1. Rodamientos SKF.
[Link]
2. Dispositivos Automation Direct.
[Link]
3. PLC Koyo.
[Link]
4. La pila OSI.
[Link]
CATÁLOGOS.
1. Servomotores SureServo.
2. Cajas reductoras Motovario.
3. Keyless Bushings Fenner Drives.
4. Encoders Allen Bradley.
5. Manual Motrona IV251.
6. OMRON finales de carrera.
7. C-MORE HMI Automation Direct.
8. Power supplies Automation Direct.
296