Reologı́a, Medios Granulares e Impresión 3D

Aurea N. Sánchez Espı́n, José Fernando Valencia Hernández

Aarón Domenzain Del Castillo C., Patricia Lecuona Sánchez
Universidad Autónoma Metropolitana - Iztapalapa Av. San Rafael Atlixco 186,
A.P. 55-534, C.P. 09340, Ciudad de México

7 de abril de 2020

Resumen

El presente trabajo consta de tres partes independientes: Se estudió la deformación de la

materia (reologı́a) para determinar la viscosidad de una muestra de miel, se diseñó un modelo
tridimensional para posteriormente construirse mediante impresión 3D y se estudiaron los sis-
temas granulares para determinar las trayectorias en un experimento de dinámica de medios
granulares.

1. Reologı́a Módulo de almacenamiento y módulo de per-

didas (comportamiento viscoelástico lineal)
La reologı́a es el estudio del flujo y la deformación
de la materia sometidas a fuerzas, y que por lo ge- Funciones complejas de viscoelasticidad no li-
neral se mide utilizando un reómetro. La medición neal
de propiedades reológicas se aplica a todos los ma- Viscosidad aparente (relación entre esfuerzo de
teriales, desde fluidos como soluciones diluidas de corte y velocidad de corte)
polı́meros y surfactantes hasta fórmulas concentra-
das de proteı́nas, y desde semi sólidos como pastas Viscosidad compleja (respuesta ante esfuerzos
y cremas hasta polı́meros derretidos o sólidos. de corte oscilatorio)
Las propiedades reológicas pueden medirse a
través de la deformación de una muestra en un volu- Centraremos nuestra atención en el estudio de la
men grande usando un reómetro mecánico, o en una viscosidad de un fluido conforme aumenta la tem-
escala microscópica mediante el uso de un viscóme- peratura. Para ello, revisaremos algunos conceptos:
tro de microcapilaridad o una técnica óptica como
la microrreologı́a.
1.2. Fluidos Viscosos y concepto de
viscosidad
1.1. Propiedades mecánicas
El concepto de viscosidad nació con Newton,
Las propiedades mecánicas estudiadas por la cuando en su obra ”Philosophiae Naturalis Prin-
reologı́a se pueden medir mediante reómetros, apa- cipia Matematica.afirmó que la resistencia ejercida,
ratos que permiten someter al material a diferentes y que surge a partir de una falta en el deslizamiento
tipos de deformaciones controladas y medir los es- de un fluido, si el resto de factores se mantienen,
fuerzos o viceversa. Las más importantes son: es proporcional a la velocidad a la que las partes
de un fluido son separadas entre sı́. De este modo,
Coeficientes de esfuerzos normales se establece la proporcionalidad existente entre el

1
esfuerzo por unidad de área (F/A) necesario para
producir un gradiente de velocidades en un fluido,
siendo la constante de proporcionalidad un factor
que describe ”la capacidad de deslizamiento de un
fluido”(más tarde ésta constante de proporcionali-
dad fue llamada viscosidad).
La viscosidad de un fluido se suele representar con
la letra griega µ y se sabe que ésta es solo función
de la temperatura, es decir: µ=µ(T):

Figura 2: viscosidades para diversos materiales a

temperatura ambiente y presión atmosférica

Figura 1: Gráfico de fluidos Neuwtonianos y No 1.3. Variación de la viscosidad con la

Newtonianos temperatura
A parte de depender de la velocidad de corte y
Para fluidos no Newtonianos la viscosidad suele del tiempo de aplicación de la misma, la viscosi-
depender de otras propiedades como la rapidez dad es fuertemente dependiente de la temperatura.
de corte. La relación entre el corte y el esfuerzo La mayorı́a de los materiales disminuyen su visco-
cortante está dada por la expresión: sidad con la temperatura; la dependencia es expo-
nencial, siendo la expresión más común la ecuación
de Arrhenius:
τ = Gγ (1) E
η(T ) = η0 e kB T (4)

donde G es el módulo de corte. Por otro lado, la Para lı́quidos más viscosos esta dependencia es
rapidez de corte y la relación existente entre ésta y mayor, y ha de tomarse mayores precauciones en el
la viscosidad η están dadas por las expresiones control de la temperatura. Respecto a los polı́me-
ros, la dependencia con la temperatura que estos
v presentan es lógicamente una función de la estruc-
γ̇ = (2) tura y del tipo de polı́mero estudiado.
L

1.4. Aplicaciones
τ = η γ̇ (3)
Control de sustancias que sean transportadas
a lo largo de un recipiente cilı́ndrico.
respectivamente. A continuación se muestra una
tabla con valores de la viscosidad para distintas Estudio del magma en vulcanologı́a: cuanto
sustancias más fluido sea el magma más tendencia va a
tener el volcán a que provoque una erupción.

2
 Caracterización de metales, elastómeros y de
polı́meros.

Producción de medicamentos, pinturas , pega-

mentos y cosméticos.

Control de calidad de los alimentos y el estudio

de su textura y consistencia.
Figura 3: Rango de velocidades de corte y tempe-
raturas usadas.
1.5. Desarrollo Experimental
Objetivo: Conocer la variación de la viscosi- 1.7. Resultados
dad de una muestra de miel a diferentes tem- A continuación se presentan los gráficos obteni-
peraturas. dos de las mediciones y sus respectivas ecuaciones:

Hipótesis: Se espera que la viscosidad de

la miel presente una disminución conforme la
temperatura aumenta además de obtener gráfi-
cos caracterı́sticos de un fluido Newtoniano.

Materiales

0.5 ml de Miel de Abeja marca Promesa

Reómetro Anton Paar

Software de Análisis Rheoplus Anton Paar

Pañuelos
Figura 4: Gráfico esfuerzo cortante en función de la
1.6. Procedimiento velocidad de corte

Encendemos el equipo y procedemos a esperar a

que este alcance las especificaciones de temperatu-
ra necesarias para comenzar el procedimiento. La
muestra de miel está a presión atmosférica. Poste-
riormente se coloca el cono plato con la especifica-
ción CP35-2, es decir que el diámetro del cono sera
de 35 mm y su grado de inclinación será 2o .
Hacemos las calibraciones correspondientes y con el
control de temperatura, fijamos la temperatura a la
cual estará expuesta la muestra. Además, damos el
valor inicial y final de la velocidad de corte. Verti-
mos la muestra en el área de medición e indicamos
al equipo que puede bajar la pieza y comenzar la
medición. A continuación se presenta una tabla con
las temperaturas y rangos en la velocidad de corte Figura 5: Tabla de funciones del esfuerzo de corte
a las cuales se expuso la muestra: respecto a la velocidad de corte.

3
1.8. Análisis de Resultados 2.1. Aplicaciones
La tecnologı́a de impresión 3D puede utilizarse
En el gráfico anterior puede notarse que la
para crear todo tipo de cosas, desde prototipos y
tendencia es lineal.Esto quiere decir que, si compa-
piezas simples hasta productos finales altamente
ramos con la Fig.1 la miel es, en buena medida, un
técnicos, como piezas para aeronaves,autopartes,
fluido Newtoniano.Aplicando lı́neas de tendencia a
réplicas de objetos, implantes médicos e incluso
cada una de las rectas que se muestran, obtenemos
órganos artificiales que se producen con capas de
las ecuaciones presentadas en la Fig.5 y sus corres-
células humanas, etc.
pondientes viscosidades (dadas por la pendiente de
En el mercado existen más de 60 tipos de mate-
las rectas). notamos que, en efecto, a medida que
riales para impresión 3D, que gracias a sus carac-
la temperatura aumenta, la viscosidad disminuye
terı́sticas y propiedades fı́sico-quı́micas, posibilitan
siendo en las temperaturas de 20o C,22o C y 24o C
la creación de prototipos perfectos, de gran preci-
donde es más notorio pues hay una diferencia
sión, excelente nivel de detalle y aplicables casi a
media de aproximadamente 0.8651 Pa.s en tanto
todos los sectores industriales.
que para las temperaturas de 26o C y 28o C se
nota una disminución de apenas (0,3668 P a.s). A
medida que la temperatura aumenta, la separación 2.2. Desarrollo Experimental
entre gráficas disminuye, de modo que, a cierta Objetivo: Fabricar una pieza mediante el
temperatura, habrá una viscosidad mı́nima. método de impresión 3D.

Materiales
Computadora
1.9. Conclusiones
Ácido poliláctico
La viscosidad de la miel es de aproximadamente
η=10,000 CPs = 10 Pa.s cuyo dato es presentado Impresora 3D
en la Fig. 2. Este valor está dado para una tempe-
ratura ambiente (25o C). Con nuestros resultados,
notamos que la viscosidad entre las temperaturas
24o C y 26o C está muy alejada de la reportada en
la Fig. 2. Esto puede ser debido a que la miel que
se utilizó es de diferente elaboración y puede influir
en la viscosidad de ésta. Por otro lado, vemos que
las gráficas se apegan bastante bien a las reporta-
das para un fluido newtoniano, por lo cual podemos
decir que la miel tuvo un comportamiento de éste
tipo esto aunado a que su viscosidad disminuyó de
manera considerable con tan solo 2o C de diferencia
entre temperaturas. De modo que se obtuvieron los
resultados esperados.

Figura 6: Impresora 3D

2. Impresión 3D
2.3. Procedimiento
La impresión 3D, también conocida como manu- Comenzamos por pensar el uso que le vamos
factura por adición, es un proceso por el cual se a dar a la pieza que vamos a fabricar, para
crean objetos fı́sicos colocando un material por ca- ası́ poder determinar caracterı́sticas como: su
pas en base a un modelo digital. forma, la presencı́a de agujeros etc.

4
 El programa de diseño de la pieza que se quie- pero que sigan siendo lo suficientemente resisten-
re imprimir se llama TinkerCad, que es una tes como para poder ser sometidos a colisiones sin
sencilla aplicación gratuita en lı́nea, de di- ser deformados, lo cual resulta ideal para experi-
seño e impresión 3D. En este programa escoge- mentos de dinámica de muchos cuerpos pequeños,
mos la forma de un rectángulo de dimensiones o dinámica de medios granulares.
45mm × 19mm × 6mm.
De la misma manera se eligió ponerle texto, 2.5. Conclusiones
con la leyenda Equipo 4
Es posible construir prácticamente cualquier
Después de tener el diseño, este se mandará a cuerpo geométrico por medio de la impresión 3D,
el programa que controla especificaciones sobre lo cual resulta muy útil cuando se necesitan mu-
la densidad y el relleno de la pieza, llamado chas piezas de alguna forma caracterı́stica con cier-
Slic3r. Para finalizar el proceso de diseño el ta densidad y resistencia, como en el caso de los
programa transforma el archivo en un Formato experimentos de dinámica de medios granulares.
GCode para que la impresora lo pueda leer y
ejecutar. 3. Medios Granulares
Se carga el archivo en la impresora y se espe- Los medios granulares o sistemas granulares son
ra a éstas llegue a la temperatura ideal para aquellos formados por objetos macroscópicos lla-
comenzar el proceso. mados granos, que interaccionan por medio de con
tactos temporales o permanentes, ası́ como por me-
La impresora crea el modelo de capa en capa
dio de algún otro potencial. Las caracterı́sticas más
esparciendo una capa de la sección de la pieza.
importantes de los medios granulares es que las
El proceso es repetido hasta que todas las ca- partı́culas interactúan de manera disipativa; me-
pas han sido impresas. diante colisiones inelásticas.
En los medios granulares, el peso de los granos es
demasiado elevado como para que su energı́a térmi-
2.4. Resultados y Discusión ca resulte en un desplazamiento significativo y por
lo tanto no presentan movimiento browniano ob-
servable. Es habitual utilizar el concepto de “tem-
peratura granular” como una medida de la energı́a
cinética asociada a las trayectorias de las partı́culas
que conforman el medio.

3.1. Disipación de la Energı́a

En un sistema aislado la relación entre las
energı́as cinéticas antes y después del choque en-
tre dos partı́culas que colisionan está dada por:

Ecf = Eci e2 (5)

Donde Ecf es la energı́a cinética de las partı́culas

Figura 7: Pieza de ácido poliláctico impresa en 3D
después de la colisión, Eci antes de la colisión y e
con la leyenda Equipo 4, a un rellenado del 20 %.
el coeficiente de restitución de la colisión.
En un medio granular, el número de colisiones
Se obtuvo una pieza muy ligera pero bastante re- por unidad de tiempo puede llegar a ser muy eleva-
sistente, a pesar de haber tenido un rellenado del do y por eso ser fuertemente disipativo. Estos siste-
20 %, por lo que es posible construir muchos objetos mas requieren ser suministrados constantemente de
con una densidad tal que floten en algún lı́quido, energı́a, en este caso, para el sistema experimental,

5
esa energı́a se suministra mediante las vibraciones cinética del sistema disminuirá a menos que se le su-
que producen la agitación de las partı́culas [2] . ministre energı́a. Para que el sistema granular sea
ergódico (pero aún lo suficientemente simple como
3.2. Estructura y Caracterización para estudiarlo con modelos simplificados), deberá
tener un suministro contı́nuo de energı́a. Si el pro-
3.2.1. Sección de Área ceso de transferencia de energı́a es estocástico, en-
tonces el sistema cumplirá con la hipótesis de ergo-
Es el área total ocupada por las partı́culas sobre
dicidad [5] y por lo tanto, para una observable arbi-
el área total del sistema , definida como η
traria fˆ, el promedio en el ensemble microcanónico
πσρ coincidirá con el promedio temporal [3] :
η= (6)
4
ˆ τ ˆ
3.3. Función de distribución radial 1
lı́m f (t)dt ∼ pfˆ(~q, p~)ρ(~q, p~)
d~qd~
τ →∞ τ
g(r) 0

Describe la variación de la densidad como fun- En experimentos con medios granulares es po-
ción de la distancia medida desde una partı́cula de sible inyectar energı́a mecánica a los sistemas de
referencia en un sistema de muchas partı́culas. manera estocástica de diversas maneras, tal como
La función de distribución radial se determina vibraciones mecánicas o perturbaciones a los cam-
usualmente calculando la distancia entre todos los pos eléctrico y magnético. Además, los granos pue-
pares de partı́culas y colocando los datos en un his- den presentar alguna interacción, lo que modifica
tograma. Entonces, el histograma se normaliza con la evolución de todo el sistema.
respecto a un gas ideal completamente descorrela-
cionado. 3.5. Desarrollo Experimental
Podemos definir la función de distribución ra-
dial, g(r),como una densidad local normalizada, di- Objetivos:
vidiendo la densidad local por la global:
1. Identificar los centros de las partı́culas.
ρ(r)
g(r) = (7) 2. Determinar las trayectorias que siguen las
ρ
partı́culas durante el experimento.

Hipótesis: Las trayectorias de los centros de

partı́cula son útiles para determinar la fı́sica
del sistema.

3.5.1. Procedimiento
Previamente se realizó un experimento de
dinámica de sistemas granulares con discos impre-
sos en 3D, con un diámetro de 2.1 cm y a un re-
llenado de 15 %, con imanes colocados encima, los
cuales se encontraban suspendidos en agua, en un
contenedor que vibraba a una frecuencia constante.
Figura 8: g(r) vs Distancia (en términos del diame- Obtuvimos un video de 36 segundos a 24 cuadros
tro). por segundo, con una resolución de 383×285 pı́xe-
les, en el que se muestran 9 discos moviéndose y
3.4. Dinámica de Medios Granulares formando clusters. Posteriormente, el video se ana-
lizó usando el software ImageJ, utilizando el módu-
La dinámica de medios granulares está goberna- lo Particle Tracker 2D/3D del plugin Mosaic [4] , por
da por colisiones disipativas, por lo que la energı́a lo que se obtuvieron las trayectorias de los 9 discos.

6
Figura 9: Primer cuadro del video del experimento. Figura 11: Trayectorias de las partı́culas determi-
nadas con ImageJ y Mosaic.

3.6. Resultados y Discusión

3.7. Conclusiones
Se obtuvo el archivo con las posiciones de cada Se obtuvieron las trayectorias de cada disco como
uno de los 9 discos en cada punto del tiempo (866 puntos en el tiempo, de lo cual se puede derivar
puntos para cada disco). prácticamente toda la fı́sica del sistema. Mediante
alguna herramienta de análisis de datos es posible
determinar la distribución de velocidades, lo que
proporcionarı́a una gran cantidad de información
sobre la dinámica del sistema.

Referencias
[1] A. T. Cadillo-Martı́nez and R. Sánchez. Ex-
perimental velocity distributions in a granular
submonolayer. Physica A: Statistical Mechanics
and its Applications, 465:221–228, Jan. 2017.

[2] I. del Carmen Romero Sánchez. Sistemas Gra-

Figura 10: Localización de los centros de partı́cula nulares Monodispersos Confinados Cuasi-1D y
estimados con ImageJ y Mosaic (Parámetros de de- Cuasi-2D. 2016.
tección de partı́cula: radius=25, cutoff radius= 0.0, [3] U. Lucia. Probability, ergodicity, irreversibility
percentile= 15.0, displacement=20.0, linkrange=1) and dynamical systems. Proceedings of the Ro-
yal Society A: Mathematical, Physical and En-
Teniendo las posiciones de las partı́culas para ca- gineering Sciences, 464(2093):1089–1104, Jan.
da tiempo, es posible calcular su derivada temporal 2008.
(la velocidad) mediante, por ejemplo, el método de [4] I. Sbalzarini and P. Koumoutsakos. Feature
Euler, el cual aproxima una curva por medio de pe- point tracking and trajectory analysis for video
queñas rectas recurriendo a la definición de la deri- imaging in cell biology. Journal of Structural
vada escrita para un incremento finito. Una vez de- Biology, 151(2):182–195, Aug. 2005.
terminadas las velocidades, serı́a posible construir
una distribución mediante un histograma de fre- [5] H. S. D. R. Walter Greiner, Ludwig Neise. Ther-
cuencias y comparar con distribuciones conocidas, modynamics and statistical mechanics. Classi-
como la de Maxwell [1] . cal theoretical physics. Springer-Verlag, 1995.