3.3.3 Valores en los puntos de integración

Para analizar los resultados, emplearemos la función gef_hist_intpoints descrita en otros tutoriales.

En este caso, conviene por ejemplo chequear la naturaleza no drenada de la respuesta, es decir, verificar que el cambio de volumen $\varepsilon_v= \varepsilon_{xx} + \varepsilon_{yy} + \varepsilon_{zz}$ es nulo durante todo el ensayo¹. Para ello, podemos emplear el código:

>> data.filename='TRXNMC2';
>> data.win=[.5 1. .5 1.];
>> data.var=1:9;
>> data.pos=0;
>> data.group=1;
>> [POS,HIST]=gef_hist_intpoints(data);
>> Ezz=squeeze(HIST(1,6,:));
>> Eyy=squeeze(HIST(1,5,:));
>> Exx=squeeze(HIST(1,8,:));

>> figure
>> plot(-Ezz,-(Ezz+Eyy+Exx),'LineWidth',2)
>> xlabel('\epsilon_{zz} [-]')
>> ylabel('\epsilon_v [-]')
>> set(gca,'YDir','Reverse')
>> grid on

que corresponde a un valor razonablemente pequeño para verificar que el comportamiento es no drenado.

Notas al pie

... ensayo¹

en este caso se debe considerar la deformación en la dirección $x$ ya que se trata de un problema de simetría de revolución y no en deformaciones planas