abaqustomatlab.htm

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

                       end
                       [vector,status]=read_float(fid);
                       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
     <A HREF="constraint.html">constraint</A>
     end

save outfil
clear

</PRE>
<P></P>
</BODY>
</HTML>