Meu SciELO
Serviços Personalizados
Artigo
Indicadores
- Citado por SciELO
Links relacionados
- Similares em SciELO
Compartilhar
Revista Ciencias Técnicas Agropecuarias
versão On-line ISSN 2071-0054
Rev Cie Téc Agr vol.30 no.1 San José de las Lajas jan.-mar. 2021 Epub 01-Jan-2021
SOFTWARE
Modelo y software para el cálculo de parámetros en discos centrífugos esparcidores de fertilizante
IUniversidad Agraria de La Habana (UNAH), Facultad de Ciencias Técnicas, Centro de Mecanización Agropecuaria (CEMA), San José de las Lajas, Mayabeque, Cuba.
IIUniversidad Autónoma Chapingo, Texcoco, Edo. México, Estados Unidos Mexicanos.
A escala internacional, la distribución a voleo de fertilizantes minerales se efectúa fundamentalmente con el empleo de esparcidores centrífugos de disco con paletas, debido a su alta capacidad de distribución, acompañado de simplicidad constructiva y de operación. No obstante, estas ventajas, se requiere una alta exigencia en la correcta selección de los parámetros que intervienen en el proceso de trabajo de estos órganos, como requisito fundamental para garantizar una alta uniformidad de distribución del fertilizante en el campo, que garantice los requerimientos productivos correspondientes. En el presente trabajo se presenta un modelo que ha sido elaborado con vistas a determinar la característica del movimiento de las partículas de fertilizante sobre el disco esparcidor dotado de paletas inclinadas con relación a la dirección radial. El análisis se realiza obteniendo la ecuación diferencial de la dinámica del movimiento de una partícula mediante la aplicación de las leyes de la mecánica clásica newtoniana. Las ecuaciones obtenidas de la solución de la ecuación diferencial del movimiento son programadas empleando el software Mathcad, versión 2000 Professional, lo cual posibilita evaluar el modelo con datos concretos y obtener como salida, la velocidad y ángulo de salida del fertilizante en función de parámetros de entrada tales como: el ángulo de inclinación de las paletas, el coeficiente de fricción entre el fertilizante y el material del disco y las paletas, las dimensiones y velocidad de rotación del disco y las coordenadas de la zona de caída del fertilizante desde la boca de la tolva.
Palabras-clave: Paletas inclinadas; uniformidad de distribución; software Mathcad
INTRODUCCIÓN
A escala internacional, la distribución a voleo de fertilizantes minerales se efectúa fundamentalmente con el empleo de esparcidores centrífugos de disco con paletas (Fig. 1), debido a su alta capacidad de distribución, acompañado de simplicidad constructiva y de operación. No obstante, estas ventajas, se requiere una alta exigencia en la correcta selección de los parámetros que intervienen en el proceso de trabajo de estos órganos, como requisito fundamental para garantizar una alta uniformidad de distribución del fertilizante en el campo, que garantice los requerimientos productivos correspondientes.
En este sentido, la modelación del movimiento de las partículas de fertilizante sobre el disco esparcidor ha sido objeto de estudio en diferentes épocas por diferentes autores Turbin et al. (1967); Villette et al. (2005); Cerović et al. (2018), con vistas a determinar la interrelación entre los parámetros que intervienen en este proceso, tales como: el diámetro y velocidad angular del disco esparcidor; la posición de las paletas y el lugar de caída o alimentación del fertilizante sobre el disco, entre otros.
El funcionamiento de una esparcidora centrífuga de fertilizantes presenta tres fases o etapas. En la primera el fertilizante es dejado caer por gravedad desde una boca o compuerta situada en el fondo de la tolva hacia una zona cercana al centro del disco (Fig. 1). Durante la segunda etapa, el material a distribuir es transportado hacia el borde del disco, bajo la acción de un sistema de fuerzas, encabezadas por la fuerza centrífuga y sometido a la acción de las fuerzas de fricción con las superficies del propio disco y de las paletas. En la tercera fase, el material es lanzado al campo, estando su trayectoria determinada por un vuelo balístico bajo la acción de la densidad y velocidad del aire circundante.
Es obvio que una adecuada predicción y control del movimiento de las partículas en la segunda fase, partiendo de una definición adecuada del punto de caída de fertilizante desde la compuerta de la tolva, permitirá ajustar, tanto la velocidad de salida, como el ángulo de dispersión en el campo de las partículas esparcidas.
La modelación de la segunda fase ha sido abordada por diferentes autores. Villette et al. (2005) elaboran un modelo analítico para describir el movimiento de partículas sobre un disco cóncavo equipado con paletas planas. El modelo permite establecer la relación entre los componentes horizontales radial y tangencial de la velocidad de salida, aunque no es aplicable directamente a discos planos. Cerović et al. (2018) analizan el movimiento de una partícula de fertilizante ideal, esférica y homogénea a lo largo de una paleta recta fijada a un disco plano rotatorio utilizado en las esparcidoras centrífugas de fertilizantes minerales. El análisis se realiza sobre un sistema de referencia no inercial, aplicando las leyes de la mecánica clásica. Como resultado arriba a un sistema de ecuaciones diferenciales ordinarias homogéneas de segundo orden, cuya solución representa una aproximación del movimiento relativo real de una partícula de fertilizante a lo largo de una paleta recta fijada en posición radial al disco que rota a velocidad angular constante. El modelo es de utilidad para la optimización de los parámetros de este tipo de distribuidor de fertilizantes, aunque no es aplicable con exactitud en el caso de paletas con disposición inclinada con relación a la dirección radial. Para este análisis, Cerović, et.al. se apoyan en estudios anteriores efectuados por Aphale et al. (2003); Dintwa et al. (2004); Villette et al. (2005), aunque ya estudios al respecto realizados en el siglo XX alcanzaban un alto nivel de desarrollo, tales como los realizados por Mennel y Reece (1962); Turbin et al. (1967); Griffis et al. (1983); Olieslagers et al. (1996). Este último posibilita el cálculo del movimiento de las partículas cuando el disco está provisto de paletas inclinadas con relación a la dirección radial, habiendo servido de base para la realización del modelo y software que se expone en este trabajo, siendo completado hasta la determinación del abanico de dispersión de las partículas de fertilizante a la salida del disco esparcidor. La tercera fase correspondiente al vuelo de la partícula una vez que es impulsada por el disco, ha sido abordada por diferentes autores Walker et al. (1997); Van Liedekerke et al. (2009); Cool et al. (2014, 2016), no constituyendo objeto de estudio en este trabajo, cuyo objetivo es exponer un modelo mecánico matemático que describa la dinámica del movimiento de las partículas sobre un disco distribuidor de fertilizantes con paletas inclinadas, así como un software que, basado en la aplicación del modelo, viabilice el cálculo de los parámetros geométricos y cinemáticos que determinan la velocidad de salida y el abanico de dispersión del fertilizante a la salida de los discos.
MATERIALES Y MÉTODOS
El modelo que se presenta es elaborado con vistas a determinar la característica del movimiento de las partículas de fertilizante sobre el disco esparcidor y en contacto con las paletas, lo que permitirá determinar la magnitud y dirección de la velocidad con que éstas abandonan el disco en función de los parámetros geométricos y cinemáticos de éste.
Para la elaboración del modelo se parte de considerar el disco plano, aunque también se emplean discos cónicos. En cuanto a las paletas, las mismas pueden estar orientadas en la dirección radial, así como con inclinación contraria a la rotación del disco (paletas atrasadas) y con inclinación a favor de la rotación del disco. Se selecciona para la elaboración del modelo este último caso, por considerarse el más eficiente desde el punto de vista de lograr una mayor velocidad del lanzamiento de las partículas con la misma velocidad de rotación del disco.
El análisis se realiza obteniendo la ecuación diferencial de la dinámica del movimiento de una partícula mediante la aplicación de las leyes de la mecánica clásica (2da. Ley de Newton). El análisis se efectúa sobre un sistema de referencia no inercial que, rota conjuntamente con el disco, por lo que es necesario tener en cuenta las fuerzas “ficticias” (fuerza centrífuga y fuerza de Coriolis) que actúan sobre la partícula en este sistema móvil. La solución de la ecuación diferencial del movimiento se efectúa mediante métodos clásicos de análisis matemático.
Las ecuaciones obtenidas de la solución de la ecuación diferencial del movimiento son programadas empleando el software Mathcad, versión 2000 Professional, lo cual posibilita evaluar el modelo con datos concretos y obtener como salida, la velocidad y ángulo de salida del fertilizante en función de parámetros de entrada tales como: el ángulo de inclinación de las paletas, el coeficiente de fricción entre el fertilizante y el material del disco y las paletas, las dimensiones y velocidad de rotación del disco y las coordenadas de la zona de caída del fertilizante desde la boca de la tolva.
DESARROLLO
En la Fig. 2 se muestra las fuerzas que actúan sobre una partícula de fertilizante en su interacción con el disco y las paletas del distribuidor centrífugo.
Como se aprecia en la figura, sobre la partícula actúan las siguientes fuerzas:
mg - Peso de la partícula, que actúa perpendicularmente a la superficie del disco (Fig. 5.6a), igual al producto de la masa de la partícula (m) y la aceleración de la gravedad (g);
Fcf - Fuerza centrífuga, la cual está dirigida en la dirección y sentido del radio-vector
siendo su magnitud:
Fco - Fuerza de Coriolis, dada por:
siendo su dirección perpendicular a los vectores
Fpcf - Fuerza de fricción entre la paleta y la partícula producto de la fuerza centrífuga:
donde
Fpco - Fuerza de fricción entre la paleta y la partícula producto de la fuerza de Coriolis:
Fpg - Fuerza de fricción entre el disco y la partícula:
Por lo general el material del disco es el mismo de las paletas, en cuyo caso el coeficiente de fricción entre el disco y la partícula
Planteando la 2da. Ley de Newton en el sistema no inercial se obtiene la ecuación diferencial del movimiento de la partícula en la dirección de la paleta:
De la Fig. 2 b) se aprecia que:
donde:
ro - radio vector inicial, dirigido desde el centro del disco hasta el comienzo de la paleta;
Asimismo es posible plantear que:
Por otro lado, el coeficiente de fricción
Sustituyendo 9, 10 y 11 en 8 y efectuando algunas transformaciones se arriba a la siguiente expresión para la ecuación diferencial del movimiento relativo de la partícula sobre el disco:
que es una ecuación diferencial ordinaria de 2do orden, lineal y no homogénea con coeficientes constantes, cuya solución general puede ser determinada como la suma de la solución de la ecuación homogénea y la solución particular de la ecuación no homogénea.
La parte homogénea de esta ecuación se expresa como:
teniendo la solución general la forma:
siendo C1 y C2 las constantes de integración, mientras que λ1 y λ2 son las raíces de la ecuación característica:
Una solución particular para la ecuación no homogénea 12 puede ser determinada aplicando el método de la constante indeterminada:
Diferenciando 19 y sustituyendo en la ecuación 12 se obtiene la siguiente solución particular de la ecuación diferencial no homogénea:
Sustituyendo 15 y 20 en 13, la solución de la ecuación diferencial 12 toma la siguiente forma:
Finalmente se evalúa la ecuación 5.25 para las condiciones iniciales, correspondiente a la posición de comienzo de las paletas, donde para t = 0;
Una gráfica de esta expresión, evaluada para determinadas condiciones, se muestra en la Fig. 5.7, apreciándose que el tiempo de permanencia de las partículas sobre el disco es del orden de las centésimas de segundo.
Derivando la expresión 22 con respecto al tiempo, se obtiene la velocidad relativa de las partículas a lo largo de las paletas del distribuidor:
La magnitud del radiovector que une el centro del disco con la posición instantánea de las partículas, se determina como sigue:
Sustituyendo la expresión 22 de
En ese tiempo, el disco habrá girado un ángulo
Evaluando estas expresiones para los puntos de alimentación 1 y 2 (Fig. 4) que forman una franja bf al inicio de la paleta, es posible obtener los ángulos (
El módulo de la velocidad absoluta se determina por la expresión:
La velocidad de arrastre
El abanico de dispersión del fertilizante (
donde:
La programación en Mathcad de las ecuaciones obtenidas se muestra a través de un ejercicio que se expone a continuación, mostrándose las capturas de pantalla de la corrida del programa que ha sido denominado “CENTRIFERT”:
Ejercicio demostrativo:
Determine la zona (bf) de colocación del fertilizante sobre un disco centrífugo con paletas adelantadas, de manera de obtener un abanico de dispersión de 90o ± 2o en dirección opuesta al avance de la máquina. Se conocen los siguientes datos:
COMENTARIO SOBRE LA CORRIDA DEL PROGRAMA:
Como se aprecia de la corrida del programa, se obtienen en forma gráfica (gráficos en coordenadas polares) los ángulos correspondientes a los puntos de salida del fertilizante en las posiciones extremas 1 y 2 correspondientes a la banda de posicionamiento del fertilizante bf = 7.945 cm. Ahora bien, como se ha explicado, la banda de aplicación o lanzamiento del fertilizante (zona sombreada en la Fig. 5) estará enmarcada por las direcciones de las velocidades absolutas (Vab1 y Vab2), que han sido determinadas sobre la base de los valores de las velocidades relativas (Vr1 y Vr2) y la velocidad de arrastre Va. En la Fig. 5, el ángulo βo1 ha sido retradado 90o aprosimadamente, para lograr que el abanico de lanzamiento del fertilizante esté situado en dirección contraria a la del movimiento de la máquina.
CONCLUSIONES
Se obtiene un modelo analítico que describe el movimento de partículas sobre un disco esparcidor de fertilizantes del tipo centrífugos con paletas rectas inclinadas con relación a la dirección radial;
Como resultado de la programación en Mathcad del modelo se viabiliza el cálculo de los paràmetros de salida del fertilizante (velocidad de salida y ángulo de dispersión) en función de parámetros de entrada tales como: el ángulo de inclinación de las paletas, el coeficiente de fricción entre el fertilizante y el material del disco y las paletas, las dimensiones y velocidad de rotación del disco y las coordenadas de la zona de caída del fertilizante desde la boca de la tolva;
Como resultado de la ejecución de un ejercicio demostrativo, aplicando el programa “CENTRIFERT”, se obtiene un ángulo de dispersión de las partículas θs ≈ 90°. El ancho de la banda de alimentación del fertilizante en el disco resultó bf≈7.5 cm.
REFERENCES
APHALE, A.; BOLANDER, N.; PARK, J.; SHAW, L.; SVEC, J.; WASSGREN, C.: “Granular fertilizer particle dynamics on and off a spinner spreader”, Biosystem Engineering, 85: 319-329, 2003, ISSN: 1537-5110, e-ISSN: 1537-5129. [ Links ]
CEROVIĆ, V.B.; DRAGAN, V.P.; RADE, L.R.; SAŠA, R.B.; VUKOVIĆ, A.: “On the Fertilizer Particle Motion Along the Vane of a Centrifugal Spreader Disc Assuming Pure Sliding of the Particle”, Journal of Agricultural Sciences, 63(1): 83-97, 2018, ISSN: 0021-8596, DOI: https://dx.doi.org/10.2298/JAS1801083C. UDC: 631.33.022. [ Links ]
COOL, R.S.; PIETERS, G.J.; ACKER, V.J.; VAN DEN, B.J.; MERTENS, C.K.; NUYTTENS, R.E.D.; VAN DE GUCHT, C.T.; VANGEYTE, J.: “Determining the effect of wind on the ballistic flight of fertilizer particles”, Biosystems Engineering, 151: 425-434, 2016, ISSN: 1537-5110, e-ISSN: 1537-5129. [ Links ]
COOL, R.S.; PIETERS, J.; MERTENS, C.K.; HIJAZI, B.; VANGEYTE, J.: “A simulation of the influence of spinning on the Ballistic flight of spherical fertilizer grains”, Computers and Electronics in Agriculture, 105: 121-131, 2014, ISSN: 0168-1699. [ Links ]
DINTWA, E.; TIJSKENS, E.; OLIESLAGERS, R.; DE BAERDEMAEKER, J.; RAMON, H.: “Calibration of a spinning disc spreader simulation model for accurate site specific fertilizer application”, Biosystem Engineering, 88: 49-62, 2004, ISSN: 1537-5110, e-ISSN: 1537-5129. [ Links ]
GRIFFIS, C.L.; RITTER, D.W.; MATTHEWS, E.J.: “Simulation of rotary spreader distribution patterns”, Transactions of the ASAE, 26(1): 33-37, 1983, ISSN: 2151-0032. [ Links ]
MENNEL, F.I.; REECE, A.R.: “The theory of the centrifugal distributor: II. Motion on the disc, off-center feed”, Journal of Agricultural Engineering Reearch, 7(4): 345-353, 1962. [ Links ]
OLIESLAGERS, R.; RAMON, H.; DE BAERDEMAEKER, J.: “Calculation of fertilizer distribution patterns from a spinning discspreader by means of a simulation model”, . Journal of Agricultural Engineering Reearch, 63: 137-152., 1996. [ Links ]
TURBIN, B.G.; LURIÉ, A.B.; GRIGORIEV, C.M.; IVANOVICH, Z.M.; MELNIKOV, C.B.: Máquinas Agrícolas. Teoría y Cálculo Tecnológico, Ed. Machinostroienie, En idioma ruso ed., Leningrado, Rusia, 1967. [ Links ]
VAN LIEDEKERKE, P.; THIJSKENS, E.; DINTWA, E.; RIOUAL, F.; VANGEYTE, J.; RAMON, H.: “DEM simulations of the particle flow on a centrifugal fertilizer spreader”, Powder Technology, 190: 346-360, 2009, ISSN: 0032-5910. [ Links ]
VILLETTE, S.; COINTAULT, F.; PIRON, E.; CHOPINET, B.: “Centrifugal Spreading, an Analytical Model for the Motion of Fertiliser Particles on a Spinning Disc”, Biosystem Engineering, 92(2): 157-165, 2005, ISSN: 1537-5110, e-ISSN: 1537-5129. [ Links ]
WALKER, J.T.; GRIFT, T.E.; HOFSTEE, J.W.: “Determining effects of fertiliser particle shape on aerodynamic properties”, Transactions of the ASAE, 40: 21-27, 1997, ISSN: 2151-0032. [ Links ]
Recibido: 20 de Junio de 2020; Aprobado: 04 de Diciembre de 2020