2.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_{yy} + \varepsilon_{zz}$ es nulo durante todo el ensayo. Para ello, podemos emplear el código:

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

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

\epsfig{file=figures/ev_ezz.eps,clip=,width=.8\textwidth}
que corresponde a un valor razonablemente pequeño para verificar que el comportamiento es no drenado.



Esteban Saez 2010-10-28