abaqusto.m

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

                 	if status==0
			   stress11(ns_node,cur_inc) = vector(1);
  	        	   stress22(ns_node,cur_inc) = vector(2);
			   stress33(ns_node,cur_inc) = vector(3);
 			   stress12(ns_node,cur_inc) = vector(4);
 			end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%                 %  12 is the  key number for invariants
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

	       elseif (label==12)
	       	if out_label~=label tmp_nn=0; end
			out_label=label;	
	       		ni_node=ni_node+1;
	       		tmp_nn=tmp_nn+1;
	       		if (cur_inc==1)
	       			tmp_node(tmp_nn)=node_number;
	       			invar_node(ni_node)=node_number;
	       		end
		        [vector,status]=read_float(fid);
                 	if status==0
             		   miss(ni_node,cur_inc)   = vector(1);
			   tres(ni_node,cur_inc)   = vector(2);
			   hydr(ni_node,cur_inc)   = vector(3);
			   princ1(ni_node,cur_inc) = vector(4);
			   princ2(ni_node,cur_inc) = vector(5);
			   princ3(ni_node,cur_inc) = vector(6);
			   inv3(ni_node,cur_inc)   = vector(7);
			  		
			end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%                 %  21 is the key number for Total strain components
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


		elseif (label==21)
			if out_label~=label tmp_nn=0; end
			out_label=label;	
	       		tmp_nn=tmp_nn+1;
	       		ne_node=ne_node+1;
	       		if (cur_inc==1)
	       			tmp_node(tmp_nn)=node_number;
	       			tostrain_node(ne_node)=node_number;
	       		end

			[vector,status]=read_float(fid);
                        if status==0	
                	   strain11(ne_node,cur_inc) = vector(1);
			   strain22(ne_node,cur_inc) = vector(2);
			   strain33(ne_node,cur_inc) = vector(3);
			   strain12(ne_node,cur_inc) = vector(4);
			
			end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%                 %  22 is the  key number for plastic strains
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
			
		elseif (label==22)
		 	if out_label~=label tmp_nn=0; end
			out_label=label;	
			np_node=np_node+1;	
	       		tmp_nn=tmp_nn+1;
	       		
	       		if (cur_inc==1)
	       			tmp_node(tmp_nn)=node_number;
	       			plstrain_node(np_node)=node_number;	       			
	       		end				       			
                 	[vector,status]=read_float(fid);
                 	if status==0
			   pstrain11(np_node,cur_inc)= vector(1);
		           pstrain22(np_node,cur_inc)= vector(2);
			   pstrain33(np_node,cur_inc)= vector(3);
			   pstrain12(np_node,cur_inc)= vector(4);
			   pequst(np_node,cur_inc)   = vector(5);
		         end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%				% 104 is record number for reaction forces
%  		1 - Node number
%		2 - First component of reaction forces (x-direction)
%		3 - Second component of reaction force  (y-direction)
%     		4 - Etc
% Only the 2 first component is used in this program the first component
% is in reac_set1 ... The node number is only collected in the first
% increment
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	

        elseif (label==104),
        	if out_label~=label tmp_nn=0; end
			out_label=label;	
                reac_pt=reac_pt+1;
                tmp_nn=tmp_nn+1;
                [node_re,status]=read_integer(fid);
                if status==0
                   [vector,status]=read_float(fid);
                end,
                if status==0
                   reac_set1(reac_pt,cur_inc)=vector(1);
		   reac_set2(reac_pt,cur_inc)=vector(2);
		end,
		if (cur_inc==1) & status==0
			 tmp_node(tmp_nn)=node_re;
			 reac_node(reac_pt)=node_re;
	        end
	        		
	
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%107 is the record number for coordinates, and
% give a corresponding node number  so by comparing the node number
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	

	elseif (label==107) & (cur_inc==1)
	   	if out_label~=label tmp_nn=0; end
			out_label=label;	
	        node_cnt=node_cnt+1;
		tmp_nn=tmp_nn+1;	
		[node_no_tmp,status]=read_integer(fid);
       		if status==0
			tmp_node(tmp_nn)=node_no_tmp;
   	    	        node_coord(node_cnt)=node_no_tmp;
       	        [coord_tmp,status]=read_float(fid);
       		end
       		if status==0
       		   coord(1,node_no_tmp)=coord_tmp(1);
                   coord(2,node_no_tmp)=coord_tmp(2);
                end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%				% 101 is record number for displacements
%  		1 - Node number
%		2 - First component of displacements (x-direction)
%		3 - Second component of displacements  (y-direction)
%     		4 - Etc
% Only the 2 first component is used in this program the first component
% is in reac_set1 ... The node number is only collected in the first
% increment
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	

        elseif (label==101),
        	if out_label~=label tmp_nn=0; end
			out_label=label;	
                disp_pt=disp_pt+1;
                tmp_nn=tmp_nn+1;
                [node_di,status]=read_integer(fid);
                if status==0
                   [vector,status]=read_float(fid);
                end,
                if status==0
                   disp1(disp_pt,cur_inc)=vector(1);
		   disp2(disp_pt,cur_inc)=vector(2);
		end,
		if (cur_inc==1) & status==0
			 tmp_node(tmp_nn)=node_di;
			 disp_node(disp_pt)=node_di;
	        end
	        		
	


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%   2000   Key number for increment start record
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

	elseif (label==2000) 		
		cur_inc=cur_inc+1;
         	[vector,status]=read_float(fid);
         	
		if status==0
		   [proc,status]=read_integer(fid);
		end
		if status==0
		   [step,status]=read_integer(fid);
		end

		if status==0
		   [inc_set_tmp,status]=read_integer(fid);		
		end
		if status==0
		   inc_set(cur_inc)=inc_set_tmp;
		end
		tmp_nn=0;
		tmp_node=0;	
		node_cnt=0;
		reac_pt=0;	
		disp_pt=0;
		np_node=0;
		ne_node=0;
		ni_node=0;
		no_node=0;
		ns_node=0;
		%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%   2001   Key number for increment end record
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

		
%	elseif (label==2001) & (cur_inc==1)				
%		if (exist('tmp_node')==1) & status==0
%		    eval([ string '= tmp_node;']);
%		  end
	
	
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
% 1991 is the key number for J-integral (2nd contour value is selected)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


	elseif (label==1991)
           [crack_num,status]=read_integer(fid);
           [string1,status]=read_text(fid);
           [dum,status]=read_integer(fid);		
	   [vector,status]=read_float(fid); 			
	   if status==0
 	      J(cur_inc)=vector(2);
 	   end
	
	
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%		% 1999 is the key number for ENERGY
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	


	elseif (label==1999) 	
	        [vector,status]=read_float(fid);
	        if status==0
		   Tese(cur_inc)=vector(2);
		   Tew(cur_inc)=vector(3);
		   Tpd(cur_inc)=vector(4);
		end
	end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%		

		[status]=skip_to_end(fid);
		[items,label,status]=read_header(fid);
   end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% If  it the ABAQUS run has not finished properly then we want to
% to clear out the last read data
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear coord_tmp cur_inc dum fid hard infil items n_node tmp_nn
clear tmp_nn node_number node_re no_node ns_node
clear status step vector  tmp_nn node_cnt reac_pt np_node ne_node ni_node
clear flag string  inc_set_tmp string1 tmp_node date_1 date_2 item
clear items read 
save outfil
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% If you want to make use of the constraint files, this is the place to notise
% that some of the nodesets has names which is known by the program. In the
% ABAQUS input file the node along the ligament, can either be called BOT,
% BOTR (for rigth hand side) BOTL (for left hand side), the whole node model
% is recognised as MODEL.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     if exist('BOTR')==1
	BOT=BOTR;
     end

      if exist('BOTL')==1
	 BOT=BOTL;
      end

     if exist('BOT')==1  & exist('stress22')==1
     constraint
     end

save outfil
clear