Continuum mechanics lab: finish report for session 1
Change-Id: I0c2dfab362d37ccaa3dba1ae69c862a2cfeeca9b
diff --git a/quad8/continuummechanics/lab/p1/data/last.dat b/quad8/continuummechanics/lab/p1/data/last.dat
new file mode 100644
index 0000000..3d452ee
--- /dev/null
+++ b/quad8/continuummechanics/lab/p1/data/last.dat
@@ -0,0 +1,14 @@
+# p (N/m^2) Pendiente k
+1.00E+03 1.30E+08 1.30E-01
+5.00E+06 2.27E+08 2.27E-01
+1.00E+07 3.21E+08 3.21E-01
+1.50E+07 4.19E+08 4.19E-01
+2.00E+07 5.26E+08 5.26E-01
+2.50E+07 6.42E+08 6.42E-01
+3.00E+07 7.71E+08 7.71E-01
+3816778.91 2.05E+08 2.05E-01
+1.25E+07 3.69E+08 3.69E-01
+7.86E+03 1.30E+08 1.30E-01
+6.18E+04 1.31E+08 1.31E-01
+4.86E+05 1.40E+08 1.40E-01
+3.82E+06 2.05E+08 2.05E-01
diff --git a/quad8/continuummechanics/lab/p1/graphs/generate.bash b/quad8/continuummechanics/lab/p1/graphs/generate.bash
new file mode 100755
index 0000000..f233a3c
--- /dev/null
+++ b/quad8/continuummechanics/lab/p1/graphs/generate.bash
@@ -0,0 +1,4 @@
+#!/usr/bin/env bash
+mkdir -p ../output
+
+./last.gnu
diff --git a/quad8/continuummechanics/lab/p1/graphs/last.gnu b/quad8/continuummechanics/lab/p1/graphs/last.gnu
new file mode 100755
index 0000000..ad0b454
--- /dev/null
+++ b/quad8/continuummechanics/lab/p1/graphs/last.gnu
@@ -0,0 +1,30 @@
+#!/usr/bin/env gnuplot -c
+# == DEFINICIONS ==
+outputfile = '../output/last' # Nom de la imatge resultant (sense extensió)
+datafile = '../data/last.dat' # Nom del fitxer de dades que es vol usar
+
+# == CONFIGURACIÓ DE L'OUTPUT PEL LATEX ==
+set terminal cairolatex size 11cm, 8.5cm font ",10"
+set output outputfile.'.tex'
+
+# == CONFIGURACIÓ DEL PLOT ==
+set xlabel '$P_\text{interna} \, (\si{\pascal})$'
+set ylabel '$K \, (\si{\newton\per\meter})$'
+
+# Opcions per la llegenda:
+set key above
+set key spacing 1.5
+set key font ",8.5"
+
+set yrange [0:]
+set logscale x
+
+plot datafile u 1:3 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/p1/informe/p1.tex b/quad8/continuummechanics/lab/p1/informe/p1.tex
index c5f93a1..9c6c042 100644
--- a/quad8/continuummechanics/lab/p1/informe/p1.tex
+++ b/quad8/continuummechanics/lab/p1/informe/p1.tex
@@ -13,6 +13,7 @@
\usepackage{chemformula}
\usepackage{multicol}
\usepackage{hyperref}
+\usepackage{float}
\usepackage{pgfplotstable}
\pgfplotsset{compat=1.16}
@@ -40,17 +41,45 @@
%%%% 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}
+\date{Primavera curs 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} \]
+ \begin{itemize}
+ \item Mòdul de Young de la càpsida del virus buit (amb l'interval d'error): $Y = \SI{1.3(1)e9}{\newton\per\meter\squared}$.
+ \item Màxima pressió interna que pot suportar abans de trencar-se:
+ \begin{itemize}
+ \item Cas d'expansió no lineal: $P_{i, \text{max}} = \SI{12.6}{\mega\pascal}$.
+ \item Cas d'expansió lineal: $P_{i, \text{max}} = \SI{14.6}{\mega\pascal}$.
+ \end{itemize}
+
+ \item Taula amb els valors de la constant elàstica del virus en funció de la pressió interna:
+
+ \begin{table}[H]
+ \centering
+ \pgfplotstabletypeset[
+ columns={0, 2},
+ columns/0/.style={column name=$P \, (\si{\pascal})$},
+ columns/2/.style={column name=$k \, (\si{\newton\per\meter})$, fixed, fixed zerofill, precision=2}
+ ]{../data/last.dat}
+ \end{table}
+ \end{itemize}
+
+ \newpage
+
+ \begin{itemize}
+ \item Figura amb els valors de la constant elàstica del virus en funció de la pressió interna en escala logarítmica:
+
+ \begin{table}[H]
+ \centering
+ \input{../output/last.tex}
+ \end{table}
+
+ \item Pressió interna $P$ associada a la presència de l'ADN a l'interior del virus (amb l'interval d'error): $P = \SI{5.0(5)e6}{\pascal}$.
+ \end{itemize}
\end{document}