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}