snd_dis.m

ABAQUS is a general purpose finite element analysis program which is widely used to analyses mechani

function [li2v]=Snd_dis(stress,a,x,jint,R);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% The stress at a distance R (e.g. rs0/J=2) from the crack tip diviede
% by the yield stress is calculated in this routine, if the stress-field
% gets beyound R, the results are sat to 0.0 in order to keep the rows
% in the matrice same length.
%
%    INPUT :
%       stress  =       stress matrice 
%       a	=	crack length
%	x	=	x coordinates
%	jint	=	J integral
%       R       =       for a value of r*s_/J=2;

%

%    OUTPUT :	
%	li2v   	=	linear interpolation of the stresses at a distance R. 

[n,m]=size(stress);
m=length(jint);
for i=1:m,
   for j=1:n,
      yp(j)=stress(j,i)/s0;
      xp(j)=x(j)*s0/jint(i);
   end;
   q=1;
   for t=2:n;
      if (xp(t) > R) &  q==1 
         q=2;
         li2v(i)=yp(t)-((yp(t)-yp(t-1))/(xp(t)-xp(t-1)))*(xp(t)-R);
      elseif (t==n-1) & (q==1)
         'Warning: fitting: did not get to r/(J/S_0)=2, replaced by 0.0'
         q=2;
         li2v(i)=0.0;
      end;
   end;
end;