Add continuum mechanics' P1 and P5 lab sessions
Change-Id: I7dfbf1169b93c93b9b182c958694faf021c2a70c
diff --git a/quad8/continuummechanics/lab/.gitignore b/quad8/continuummechanics/lab/.gitignore
index 19d0629..3f57aa0 100644
--- a/quad8/continuummechanics/lab/.gitignore
+++ b/quad8/continuummechanics/lab/.gitignore
@@ -1,2 +1,4 @@
# Don't include compiled PDFs
*.pdf
+# Nor generated files
+output
diff --git a/quad8/continuummechanics/lab/p1/informe/p1.tex b/quad8/continuummechanics/lab/p1/informe/p1.tex
new file mode 100644
index 0000000..c5f93a1
--- /dev/null
+++ b/quad8/continuummechanics/lab/p1/informe/p1.tex
@@ -0,0 +1,56 @@
+\documentclass[11pt,a4paper]{article}
+\usepackage[utf8]{inputenc}
+\usepackage[catalan]{babel}
+\usepackage{fancyhdr}
+\usepackage{graphicx}
+\usepackage[labelfont=bf]{caption}
+\usepackage{siunitx}
+\usepackage{geometry}
+\geometry{margin=20mm}
+\usepackage{amsmath}
+\usepackage{amsfonts}
+\usepackage{booktabs}
+\usepackage{chemformula}
+\usepackage{multicol}
+\usepackage{hyperref}
+
+\usepackage{pgfplotstable}
+\pgfplotsset{compat=1.16}
+\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},
+set thousands separator={\,}%,
+%every even row/.style={
+%before row={\rowcolor[gray]{0.9}}}, % Add this for stylish tables ;)
+%begin table=\begin{longtable},
+%end table=\end{longtable}
+}
+
+\sisetup{separate-uncertainty=true}
+
+\setlength{\parskip}{1em}
+
+\pagestyle{fancy}
+\fancyhf{}
+\rhead{Adrià Vilanova Martínez}
+\lhead{Pràctica 1}
+\rfoot{\thepage}
+
+%%%% Title %%%%
+\title{\vspace{-2ex}Pràctica 1. Determinació del mòdul de Young i la pressió interna d'un virus mitjançant modelització d'experiments d'indentació}
+\author{Adrià Vilanova Martínez (D1)\vspace{-2ex} }
+\date{Primavera curso 2020-21}
+
+\begin{document}
+ \maketitle
+
+ \section{Mòdul de Young de la càpsida del virus buit}
+
+ Per trobar el mòdul de Young de la càpsida del virus buit s'ha procedit de la següent manera: s'ha suposat que la constant elàtica té una relació monòtona creixent respecte del mòdul de Young (tal com es menciona a la teoria del guió de la pràctica), i en base a això s'ha utilitzat l'algoritme de la cerca binària per trobar el valor del mòdul de Young tal que la constant elàstica és $k_0 = \SI{0.13(1)}{\newton\per\meter}$. S'ha procedit d'igual forma per trobar els extrems de l'interval de confiança pel valor del mòdul de Young.
+
+ El valor trobat és:
+ \[ Y = \SI{1.3(1)e9}{\newton\per\meter\squared} \]
+
+
+\end{document}
diff --git a/quad8/continuummechanics/lab/p5/data/noslip.dat b/quad8/continuummechanics/lab/p5/data/noslip.dat
new file mode 100644
index 0000000..212a29e
--- /dev/null
+++ b/quad8/continuummechanics/lab/p5/data/noslip.dat
@@ -0,0 +1,10 @@
+# L (m) n \lambda_x (m) K_s*d RA* (+-10)
+2 2 2 3.141592654 1720
+2.1 2 2.1 2.991993003 1720
+2.2 2 2.2 2.855993321 1740
+2.3 2 2.3 2.731819699 1760
+2.4 2 2.400000 2.617993878 1790
+2.5 3 1.666666 3.769911184 1820
+2.6 3 1.733333 3.6249146 1780
+2.7 3 1.8 3.490658504 1760
+2.8 3 1.866666667 3.365992129 1740
diff --git a/quad8/continuummechanics/lab/p5/data/slip.dat b/quad8/continuummechanics/lab/p5/data/slip.dat
new file mode 100644
index 0000000..29e89f7
--- /dev/null
+++ b/quad8/continuummechanics/lab/p5/data/slip.dat
@@ -0,0 +1,10 @@
+# L (m) n \lambda_x (m) K_s*d RA* (+-10)
+0.8 1 1.6 3.926990817 1050
+0.9 1 1.8 3.490658504 880
+1 1 2 3.141592654 780
+1.1 1 2.2 2.855993321 720
+1.2 1 2.4 2.617993878 680
+1.3 1 2.6 2.416609734 660
+1.4 1 2.8 2.243994753 660
+1.6 1 3.2 1.963495408 680
+2 1 4 1.570796327 760
diff --git a/quad8/continuummechanics/lab/p5/graph/generate.bash b/quad8/continuummechanics/lab/p5/graph/generate.bash
new file mode 100644
index 0000000..3f55e2b
--- /dev/null
+++ b/quad8/continuummechanics/lab/p5/graph/generate.bash
@@ -0,0 +1,5 @@
+#!/usr/bin/env bash
+mkdir -p ../output
+
+gnuplot -c graph.gnu "../data/noslip.dat" "../output/noslip"
+gnuplot -c graph.gnu "../data/slip.dat" "../output/slip"
diff --git a/quad8/continuummechanics/lab/p5/graph/graph.gnu b/quad8/continuummechanics/lab/p5/graph/graph.gnu
new file mode 100755
index 0000000..9da674b
--- /dev/null
+++ b/quad8/continuummechanics/lab/p5/graph/graph.gnu
@@ -0,0 +1,28 @@
+#!/usr/bin/env gnuplot -c
+outputfile=ARG2
+datafile=ARG1
+
+y_error=10
+
+# == CONFIGURACIÓ DE L'OUTPUT PEL LATEX ==
+set terminal cairolatex size 11cm, 8cm
+set output outputfile.'.tex'
+
+# == CONFIGURACIÓ DEL PLOT ==
+set xlabel '$K_s \cdot d$'
+set ylabel '$Ra^*$'
+
+# Opcions per la llegenda:
+set key above
+set key spacing 1.5
+
+plot datafile u 4:5:(0):(y_error) w xyerr t "Dades"
+
+# == 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/p5/informe/p5.tex b/quad8/continuummechanics/lab/p5/informe/p5.tex
new file mode 100644
index 0000000..c437c0c
--- /dev/null
+++ b/quad8/continuummechanics/lab/p5/informe/p5.tex
@@ -0,0 +1,83 @@
+\documentclass[11pt,a4paper]{article}
+\usepackage[utf8]{inputenc}
+\usepackage[catalan]{babel}
+\usepackage{fancyhdr}
+\usepackage{graphicx}
+\usepackage[labelfont=bf]{caption}
+\usepackage{siunitx}
+\usepackage{geometry}
+\geometry{margin=20mm}
+\usepackage{amsmath}
+\usepackage{amsfonts}
+\usepackage{booktabs}
+\usepackage{chemformula}
+\usepackage{multicol}
+\usepackage{hyperref}
+
+\usepackage{pgfplotstable}
+\pgfplotsset{compat=1.16}
+\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},
+set thousands separator={\,}%,
+%every even row/.style={
+%before row={\rowcolor[gray]{0.9}}}, % Add this for stylish tables ;)
+%begin table=\begin{longtable},
+%end table=\end{longtable}
+}
+
+\sisetup{separate-uncertainty=true}
+
+\setlength{\parskip}{1em}
+
+\pagestyle{fancy}
+\fancyhf{}
+\rhead{Adrià Vilanova Martínez}
+\lhead{Práctica 5}
+\rfoot{\thepage}
+
+%%%% Title %%%%
+\title{\vspace{-2ex}Pràctica 5. Anàlisi de la inestabilitat de Rayleigh-Bénard mitjançant dinàmica de fluids computacional}
+\author{Adrià Vilanova Martínez (D1)\vspace{-2ex} }
+\date{Primavera curs 2020-21}
+
+\begin{document}
+ \maketitle
+
+ \begin{center}
+ \centering
+ \input{../output/noslip.tex}
+ \captionof{figure}{Gràfica de la corba d'estabilitat en el cas en què les condicions de contorn són de \textit{no-slip}.}
+ \end{center}
+
+ \begin{center}
+ \centering
+ \pgfplotstabletypeset[
+ columns = {0, 2, 3, 4},
+ columns/0/.style={column name=$L \, (\si{\meter})$},
+ columns/2/.style={column name=$\lambda_x \, (\si{\meter})$, fixed zerofill, precision=2},
+ columns/3/.style={column name=$K_x \, (\si{\per\meter})$, fixed zerofill, precision=2},
+ columns/4/.style={column name=$Ra^* \, (\pm 10)$}
+ ]{../data/noslip.dat}
+ \captionof{table}{Dades obtinguts a les simulacions en el cas en què les condicions de contorn són de \textit{no-slip}.}
+ \end{center}
+
+ \begin{center}
+ \centering
+ \input{../output/slip.tex}
+ \captionof{figure}{Gràfica de la corba d'estabilitat en el cas en què les condicions de contorn són de \textit{slip}.}
+ \end{center}
+
+ \begin{center}
+ \centering
+ \pgfplotstabletypeset[
+ columns = {0, 2, 3, 4},
+ columns/0/.style={column name=$L \, (\si{\meter})$},
+ columns/2/.style={column name=$\lambda_x \, (\si{\meter})$, fixed zerofill, precision=2},
+ columns/3/.style={column name=$K_x \, (\si{\per\meter})$, fixed zerofill, precision=2},
+ columns/4/.style={column name=$Ra^* \, (\pm 10)$}
+ ]{../data/slip.dat}
+ \captionof{table}{Dades obtinguts a les simulacions en el cas en què les condicions de contorn són de \textit{slip}.}
+ \end{center}
+\end{document}