2.3.2 Valores en los nodos

Como se mencionó en otros tutoriales, los valores nodales de deplazamientos y presiones son almacenados en el archivo binario en una matriz UNOD. Esta matriz matrices posee 3 índices:

• número de nodo
• número de grado de libertad (7 campos): $u_x$, $u_y$, $u_z$, $r_{xy}$, $r_{xz}$, $r_{yz}$ y $u_w$
• paso de tiempo guardado

En este ejemplo, si se desea trazar la evolución del presión de poros del nodo número 4 se puede emplear las instrucciones siguientes:

>> UNOD=gefread('BXNMC2_SAVE UNOD');
>> Time=gefread('BXNMC2_SAVE TIME');

>> figure
>> plot(10*Time,1e-3*squeeze(UNOD(4,7,:)),'LineWidth',2)
>> xlabel('step')
>> ylabel('u_w [kPa]')
>> title('Pore pressure at node 4')
>> grid on

y se obtendrá

Cabe mencionar que se ha seleccionado el nodo número 4 en forma arbitraria, pero los resultados son idénticos para todos los nodos por tratarse de un ensayo homogéneo.

Es interesante notar que esta curva pone de manifiesto la tendencia del comportamiento volumétrico del material: inicialmente contractivo ($\dot{u_w}>0$) y luego dilatante ($\dot{u_w}<0$). Notar que por tratarse de un ensayo no drenado la variación de volumen es nula.