Add continuum mechanics' P2 lab session
Change-Id: Ia6c229f2297a1ff24c27bbfb48212ba7812198ae
diff --git a/quad8/continuummechanics/lab/p2/data/gel.dat b/quad8/continuummechanics/lab/p2/data/gel.dat
new file mode 100644
index 0000000..c78324e
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/data/gel.dat
@@ -0,0 +1,16 @@
+# RPM v (mPa s) % Temperatura Shear rate (s^-1) Tau (Pa) Torque (microN)
+8 12075 16.1 21.1 1.672 20.608 10.847
+8 12096 16.1 21.1 1.672 20.608 10.847
+8 12073 16.1 21.1 1.672 20.608 10.847
+15 10606 26.5 21.1 3.135 33.920 17.853
+35 5250.7 30.7 21.1 7.315 39.296 20.683
+60 3033.5 30.3 21.2 12.540 38.784 20.413
+80 2272.3 30.4 21.2 16.720 38.912 20.480
+105 1737.9 30.4 21.2 21.945 38.912 20.480
+135 1350.5 30.5 21.2 28.215 39.040 20.548
+160 1153.4 30.8 21.2 33.440 39.424 20.750
+180 1055.2 31.2 21.3 37.620 39.936 21.019
+200 929 31 21.3 41.800 39.680 20.885
+250 743.4 31.1 21 52.250 39.808 20.952
+250 742.7 31 21 52.250 39.680 20.885
+250 745.7 31.1 21 52.250 39.808 20.952
diff --git a/quad8/continuummechanics/lab/p2/data/glicerina.dat b/quad8/continuummechanics/lab/p2/data/glicerina.dat
new file mode 100644
index 0000000..81a58d4
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/data/glicerina.dat
@@ -0,0 +1,16 @@
+# RPM v (mPa s) % Temperatura Shear rate (s^-1) Tau (Pa) Torque (microN)
+20 1013 16.9 21.1 4.280 4.343 11.386
+20 1014.8 16.9 21.1 4.280 4.343 11.386
+20 1014.9 16.9 21.1 4.280 4.343 11.386
+30 992.8 24.8 20.9 6.420 6.374 16.708
+40 992 33.1 20.9 8.560 8.507 22.299
+50 991.3 41.4 20.9 10.700 10.640 27.891
+60 990.4 49.7 20.9 12.840 12.773 33.483
+70 988.6 57.7 20.9 14.980 14.829 38.872
+80 986.9 65.8 20.9 17.120 16.911 44.329
+90 986.2 73.9 20.9 19.260 18.992 49.786
+100 980.7 81.7 21 21.400 20.997 55.041
+105 972.7 85.1 21 22.470 21.871 57.332
+120 987.1 98.7 21.1 25.680 25.366 66.494
+120 971.8 97.3 21 25.680 25.006 65.551
+120 972.4 97.3 21 25.680 25.006 65.551
diff --git a/quad8/continuummechanics/lab/p2/data/glicerina_temperatura.dat b/quad8/continuummechanics/lab/p2/data/glicerina_temperatura.dat
new file mode 100644
index 0000000..74e0160
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/data/glicerina_temperatura.dat
@@ -0,0 +1,31 @@
+# RPM v (mPa s) % Temperatura Shear rate (s^-1) Tau (Pa) Torque (microN)
+200 98.6 16.3 51.90 42.800 4.189 10.981
+200 103.3 17.2 51.00 42.800 4.420 11.588
+200 107.3 17.9 50.00 42.800 4.600 12.059
+200 110.9 18.5 49.00 42.800 4.755 12.463
+200 113.9 19.0 48.00 42.800 4.883 12.800
+200 117.7 19.6 47.00 42.800 5.037 13.205
+200 122.3 20.3 46.00 42.800 5.217 13.676
+200 126.7 21.1 45.00 42.800 5.423 14.215
+200 131.8 22.0 44.00 42.800 5.654 14.821
+200 137.9 23.0 43.00 42.800 5.911 15.495
+200 144.4 24.0 42.00 42.800 6.168 16.169
+200 151.9 25.3 41.00 42.800 6.502 17.045
+200 160.8 26.7 40.00 42.800 6.862 17.988
+200 169.9 28.3 39.00 42.800 7.273 19.066
+200 182.5 30.3 38.00 42.800 7.787 20.413
+200 199.8 33.3 36.60 42.800 8.558 22.434
+200 207.8 34.6 36.00 42.800 8.892 23.310
+200 228.1 38.0 35.00 42.800 9.766 25.601
+200 246.6 41.1 34.00 42.800 10.563 27.689
+200 270.5 45.1 33.00 42.800 11.591 30.384
+200 296.6 49.4 32.00 42.800 12.696 33.281
+200 324.7 54.1 31.20 42.800 13.904 36.447
+200 352.7 58.8 30.20 42.800 15.112 39.614
+200 366.6 61.1 29.80 42.800 15.703 41.163
+200 381.2 63.5 29.40 42.800 16.320 42.780
+200 397.1 66.2 29.00 42.800 17.013 44.599
+200 409.9 68.3 28.70 42.800 17.553 46.014
+200 426.0 70.9 28.30 42.800 18.221 47.765
+200 441.5 73.6 28.00 42.800 18.915 49.584
+200 453.9 75.6 27.70 42.800 19.429 50.932
diff --git a/quad8/continuummechanics/lab/p2/graphs/generate.bash b/quad8/continuummechanics/lab/p2/graphs/generate.bash
new file mode 100755
index 0000000..d15e9b8
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/graphs/generate.bash
@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+mkdir -p ../output
+
+./viscosity_vs_rpm.gnu '../data/gel.dat' '../output/gel_viscosity_vs_rpm' 25
+./torque_vs_rpm.gnu '../data/gel.dat' '../output/gel_torque_vs_rpm' 25
+./viscosity_vs_rpm.gnu '../data/glicerina.dat' '../output/glicerina_viscosity_vs_rpm' 2
+./torque_vs_rpm.gnu '../data/glicerina.dat' '../output/glicerina_torque_vs_rpm' 2
+./viscosity_vs_t.gnu '../data/glicerina_temperatura.dat' '../output/glicerina_viscosity_vs_t' 30
diff --git a/quad8/continuummechanics/lab/p2/graphs/torque_vs_rpm.gnu b/quad8/continuummechanics/lab/p2/graphs/torque_vs_rpm.gnu
new file mode 100755
index 0000000..4999493
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/graphs/torque_vs_rpm.gnu
@@ -0,0 +1,29 @@
+#!/usr/bin/env gnuplot -c
+# == DEFINICIONS ==
+outputfile = ARG2 # Nom de la imatge resultant (sense extensió)
+datafile = ARG1 # Nom del fitxer de dades que es vol usar
+
+# == CONFIGURACIÓ DE L'OUTPUT PEL LATEX ==
+set terminal cairolatex size 7.5cm, 5.5cm font ",10"
+set output outputfile.'.tex'
+
+# == CONFIGURACIÓ DEL PLOT ==
+set xlabel '$\omega \, (\si{\rpm})$'
+set ylabel '$T \, (\si{\micro\newton})$'
+
+set yrange [0:]
+
+# Opcions per la llegenda:
+set key above
+set key spacing 1.5
+set key font ",8.5"
+
+plot datafile u 1:7 w p pt 2 lc ARG3 t "Dades experimentals"
+
+# == CONFIGURACIÓ DE L'OUTPUT PER SVG ==
+# Això ho uso per generar també una imatge de previsualització que puc carregar
+# a l'ordinador per veure més o menys com a sortit el plot sense haver
+# d'inserir-ho al LaTeX per veure-ho.
+set terminal svg dashed size 600, 600 font "Computer Modern,Tinos,Helvetica,15"
+set output outputfile.'.svg'
+replot
diff --git a/quad8/continuummechanics/lab/p2/graphs/viscosity_vs_rpm.gnu b/quad8/continuummechanics/lab/p2/graphs/viscosity_vs_rpm.gnu
new file mode 100755
index 0000000..a78836f
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/graphs/viscosity_vs_rpm.gnu
@@ -0,0 +1,29 @@
+#!/usr/bin/env gnuplot -c
+# == DEFINICIONS ==
+outputfile = ARG2 # Nom de la imatge resultant (sense extensió)
+datafile = ARG1 # Nom del fitxer de dades que es vol usar
+
+# == CONFIGURACIÓ DE L'OUTPUT PEL LATEX ==
+set terminal cairolatex size 7.5cm, 5.5cm font ",10"
+set output outputfile.'.tex'
+
+# == CONFIGURACIÓ DEL PLOT ==
+set xlabel '$\omega \, (\si{\rpm})$'
+set ylabel '$\eta \, (\si{\milli\pascal\per\second})$'
+
+set yrange [0:]
+
+# Opcions per la llegenda:
+set key above
+set key spacing 1.5
+set key font ",8.5"
+
+plot datafile u 1:2 w p pt 2 lc ARG3 t "Dades experimentals"
+
+# == CONFIGURACIÓ DE L'OUTPUT PER SVG ==
+# Això ho uso per generar també una imatge de previsualització que puc carregar
+# a l'ordinador per veure més o menys com a sortit el plot sense haver
+# d'inserir-ho al LaTeX per veure-ho.
+set terminal svg dashed size 600, 600 font "Computer Modern,Tinos,Helvetica,15"
+set output outputfile.'.svg'
+replot
diff --git a/quad8/continuummechanics/lab/p2/graphs/viscosity_vs_t.gnu b/quad8/continuummechanics/lab/p2/graphs/viscosity_vs_t.gnu
new file mode 100755
index 0000000..4233f9d
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/graphs/viscosity_vs_t.gnu
@@ -0,0 +1,29 @@
+#!/usr/bin/env gnuplot -c
+# == DEFINICIONS ==
+outputfile = ARG2 # Nom de la imatge resultant (sense extensió)
+datafile = ARG1 # Nom del fitxer de dades que es vol usar
+
+# == CONFIGURACIÓ DE L'OUTPUT PEL LATEX ==
+set terminal cairolatex size 12cm, 8cm font ",10"
+set output outputfile.'.tex'
+
+# == CONFIGURACIÓ DEL PLOT ==
+set xlabel '$T \, (\si{\celsius})$'
+set ylabel '$\eta \, (\si{\milli\pascal\per\second})$'
+
+set yrange [0:]
+
+# Opcions per la llegenda:
+set key above
+set key spacing 1.5
+set key font ",8.5"
+
+plot datafile u 4:2 w p lc ARG3 t "Dades experimentals"
+
+# == CONFIGURACIÓ DE L'OUTPUT PER SVG ==
+# Això ho uso per generar també una imatge de previsualització que puc carregar
+# a l'ordinador per veure més o menys com a sortit el plot sense haver
+# d'inserir-ho al LaTeX per veure-ho.
+set terminal svg dashed size 600, 600 font "Computer Modern,Tinos,Helvetica,15"
+set output outputfile.'.svg'
+replot
diff --git a/quad8/continuummechanics/lab/p2/informecurt/main.tex b/quad8/continuummechanics/lab/p2/informecurt/main.tex
new file mode 100644
index 0000000..7ce09f1
--- /dev/null
+++ b/quad8/continuummechanics/lab/p2/informecurt/main.tex
@@ -0,0 +1,145 @@
+\documentclass[11pt,a4paper]{article}
+\usepackage[utf8]{inputenc}
+\usepackage[spanish]{babel}
+\usepackage{fancyhdr}
+\usepackage{graphicx}
+\usepackage[labelfont=bf]{caption}
+\usepackage{siunitx}
+\usepackage{geometry}
+\geometry{top=25mm}
+\usepackage{amsmath}
+\usepackage{amsfonts}
+\usepackage{booktabs}
+\usepackage{chemformula}
+\usepackage{multicol}
+\usepackage{hyperref}
+\usepackage{adjustbox}
+\usepackage{float}
+
+\usepackage{pgfplotstable}
+\pgfplotsset{compat=1.17}
+\pgfplotstableset{
+empty cells with={--}, % replace empty cells with ’--’
+every head row/.style={before row=\toprule,after row=\midrule},
+every last row/.style={after row=\bottomrule}%,
+%every even row/.style={
+%before row={\rowcolor[gray]{0.9}}}, % Add this for stylish tables ;)
+%begin table=\begin{longtable},
+%end table=\end{longtable}
+}
+\usepackage{chngcntr}
+\counterwithin{figure}{section}
+\counterwithin{table}{section}
+
+% Who the heck calls tables "cuadros"? We should call RAE to add an entry to thepage
+% word "tabla", so it can be used like "cuadro" ;)
+\addto\captionsspanish{%
+\def\tablename{Tabla}%
+}
+
+% Yep, it's not a SI unit, but screw it! ;)
+\DeclareSIUnit\rpm{RPM}
+
+\sisetup{separate-uncertainty=true}
+\setlength{\parskip}{1em}
+\pagestyle{fancy}
+\fancyhf{}
+\rhead{David Díaz, Adrià Vilanova}
+\lhead{Práctica 2}
+\rfoot{\thepage}
+
+%%%% Title %%%%
+\title{\vspace{-4ex}Práctica 2. Determinación de la viscosidad}
+\author{David Díaz, Adrià Vilanova}
+\date{Primavera 2020-21}
+\begin{document}
+ {\parskip=0pt
+ \maketitle
+ }
+
+ \section{Gel de baño}
+
+ \begin{table}[H]
+ \centering
+ \pgfplotstabletypeset[
+ columns={0, 1, 2, 3},
+ columns/0/.style={column name=$\omega \, (\si{\rpm})$, fixed, precision=0},
+ columns/1/.style={column name=$\eta \, (\si{\milli\pascal\per\second})$, fixed, fixed zerofill, precision=1},
+ columns/2/.style={column name=\% \textit{torque}, fixed, fixed zerofill, precision=1},
+ columns/3/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}
+ ]{../data/gel.dat}
+ \caption{Datos experimentales obtenidos para el gel de baño mediante el viscosímetro de rotación usando la varilla L4.}
+ \end{table}
+
+ \begin{figure}[ht]
+ \centering
+ \begin{minipage}{0.46\textwidth}
+ \centering
+ \input{../output/gel_viscosity_vs_rpm.tex}
+ \captionof{figure}{Viscosidad en función de la velocidad de rotación para el gel de baño.}
+ \end{minipage}\hfill
+ \begin{minipage}{0.46\textwidth}
+ \centering
+ \input{../output/gel_torque_vs_rpm.tex}
+ \captionof{figure}{\textit{Torque} en función de la velocidad de rotación para el gel de baño.}
+ \end{minipage}
+ \end{figure}
+
+ \textbf{Tipo de fluido según sus características reológicas}: pseudoplástico o \textit{shear thinning}.
+
+ \section{Glicerina}
+
+ \begin{table}[H]
+ \centering
+ \pgfplotstabletypeset[
+ columns={0, 1, 2, 3},
+ columns/0/.style={column name=$\omega \, (\si{\rpm})$, fixed, precision=0},
+ columns/1/.style={column name=$\eta \, (\si{\milli\pascal\per\second})$, fixed, fixed zerofill, precision=1},
+ columns/2/.style={column name=\% \textit{torque}, fixed, fixed zerofill, precision=1},
+ columns/3/.style={column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}
+ ]{../data/glicerina.dat}
+ \caption{Datos experimentales obtenidos para la glicerina mediante el viscosímetro de rotación usando la varilla L3.}
+ \end{table}
+
+ \vspace{4ex}
+
+ \begin{figure}[ht]
+ \centering
+ \begin{minipage}{0.45\textwidth}
+ \centering
+ \input{../output/glicerina_viscosity_vs_rpm.tex}
+ \captionof{figure}{Viscosidad en función de la velocidad de rotación para la glicerina.}
+ \end{minipage}\hfill
+ \begin{minipage}{0.45\textwidth}
+ \centering
+ \input{../output/glicerina_torque_vs_rpm.tex}
+ \captionof{figure}{\textit{Torque} en función de la velocidad de rotación para la glicerina.}
+ \end{minipage}
+ \end{figure}
+
+ \vspace{4ex}
+
+ \textbf{Tipo de fluido según sus características reológicas}: newtoniano.
+
+ \newpage
+
+ \section{Glicerina. Viscosidad en función de la temperatura}
+ \begin{table}[H]
+ \centering
+ \pgfplotstabletypeset[
+ columns={3, 1, 3, 1},
+ display columns/0/.style={select equal part entry of={0}{2}, column name=$\eta \, (\si{\milli\pascal\per\second})$, fixed, fixed zerofill, precision=1},
+ display columns/1/.style={select equal part entry of={0}{2}, column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1, column type/.add={}{||}},
+ display columns/2/.style={select equal part entry of={1}{2}, column name=$\eta \, (\si{\milli\pascal\per\second})$, fixed, fixed zerofill, precision=1},
+ display columns/3/.style={select equal part entry of={1}{2}, column name=$T \, (\si{\celsius})$, fixed, fixed zerofill, precision=1}
+ ]{../data/glicerina_temperatura.dat}
+ \caption{Datos experimentales obtenidos para la glicerina mientras se enfriaba mediante el viscosímetro de rotación usando la varilla L3, y una rotación fija de $\SI{200}{\rpm}$.}
+ \end{table}
+
+ \begin{figure}[H]
+ \centering
+ \input{../output/glicerina_viscosity_vs_t.tex}
+ \caption{Viscosidad en función de la temperatura para la glicerina.}
+ \end{figure}
+
+\end{document}