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:
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á
Es interesante notar que esta curva pone de manifiesto la tendencia del comportamiento volumétrico del material: inicialmente contractivo () y luego dilatante (). Notar que por tratarse de un ensayo no drenado la variación de volumen es nula.
Esteban Saez 2010-10-28