3.3.3 Valores en los puntos de integración

El procedimiento de extracción de valores en los puntos de integración es idéntico al descrito en §2.3.3. Por lo tanto, si se desea construir la curva del deviador $q$ con respecto a la deformación axial se puede emplear el código:

>> data.filename='TRXDMC1';
>> 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,:));
>> Q=squeeze(HIST(1,2,:))-squeeze(HIST(1,1,:));
>> figure
>> plot(-Ezz,-1e-3*Q,'LineWidth',2)
>> xlabel('\epsilon_{zz} [-]')
>> ylabel('q [kPa]')
>> grid on

y se obtendrá:

O bien para la deformación volumétrica¹ en función de la deformación axial:

>> 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

para obtener:

Notas al pie

... volumétrica¹

en este caso se debe considerar la componente $\varepsilon_{xx}$ ya que el modelo está en simetría de revolución y no en deformaciones planas