constrai.m

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Script  : Postcalc.m 
% 							    %
%   Copyright (c) 1993-1994 Annette D. Karstensen 	    %
%                         / University of Glasgow, SCOTLAND %
%							    %	
%  This is a script file used at the end of the superpost   %
%  it calculates various variables, for further explanation %
%  off the variables see variables.m which can be accessed  %
%  in matlab by > help variables		            %
% 
% OUTPUT:
%  	x_bot 		: x-coordinates ahead of the crack
%	s22_bot		: stress in 2 dircetion ahead of crack
%	r		: distance from crack tip 
%       pnd		: P/P_limit
%	Tnd		: T stress nondimensionlised by yield stress
%       lirv2		: Stresses at a distance 2J/So
%	reac2_bot	: Reaction forces ahead of the crack
%	K_I		: K_I calculated from reaction forces
%	J_el		: J-elastic calculated from K_I
%	J_pl		: J-plastic calculated from J-J_el
% 	Jpnd		: J-plastic nondimensionalised by (c So)
% 	Q_total		: Q_total lirv2-ssy at a distance 2J/So
%	Q_pl		: Q-plastic Q_total-S_mblf (2J/So)
%       st2nd		: Stresses in 2 direction at a distance
%			  2J/so divided by S_mblf ahead of the crack
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%							    
% Filtering the 2 direction stresses in the ligament        
% from the main stress matric				    
% It has been necessary to add this if exist statement to
% this program as well as to ABAQUStoMatlab, because sometimes
% this program has been used independent of the other one.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Some material properties is defined in the global space
% The yield stress s0 
% Youngs Modulus E
global s0=2e8 E=2e11;

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

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% here the the results for the  BOT node set is separated from
% the results of the MODEL

   [x_bot,s22_bot]=filtering(stress_node,stress22,BOT,coord,MODEL);
  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%							    
% Sorting the stresses after the coordinates, so stesses    
% closes to the crack tip is first and so fort              
	
   [s22_bot,r]=sorting(s22_bot,x_bot,a);
    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%							        
% Fit the 2 direction stresses to a distance r=dist.J/yield 

dist=2;
   
   lirv2=Snd_dis(s22_bot,a,r,J,dist);
   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%							        
% Filtering the 2 direction reaction forces in the ligament 
% from the main reaction force matric			    
  
   [x_bot,reac2_bot]=filtering(reac_node,reac_set2,BOT,coord,MODEL);
   
   [x_flanke,T_flanke]=filtering(stress_node,stress11,TSTRESS,coord,MODEL);

   T_flanke=min(T_flanke)/s0;
   Sum_reac=sum(reac2_bot);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%							        
% Calculated P/P_limit					    
   
   pnd=limit_load(reac2_bot,cp_num,x_bot,a);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%							        
%   Calculated the non dimensionalised T stress		    

   Tnd=T_stress(a,reac2_bot,cp_num,r);
   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%							        
% Calculates K, J-elastic and J-Plasitc, Jpnd is J-plastic  
% non dimensioanlised by the ligament and yield stress      
   
   [K_I,J_el,J_pl,Jpnd]=JK_calc(reac2_bot,J,cp_num,x_bot,a);
	
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%							        
% Q-total and Q-plastic is calculated from the MBLF,        
% st2nd is s22/smblf at a distance "dist"                    
	
   [Q_total,Q_pl,st2nd]=Q_stress(Tnd,lirv2,hard,dist);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%	
%							        
% for CCP's and DEC's the moment-to-force ratio and the     
% line of action is calculated, further T_mom is T from     
% a pure moment, taken at a distance m_dis, general this    
% distance is half the width as that will cancel out the    
% vertial tensile forces and T_ten is T from the tensile    
%  net force.  T_dif is Tnd-T_ten			    
% 25/7-94 (ADK) changed to be calculated for SEC and SEB    
% aswell     						    
if cp_num==1 | 4  
    m_dis=500;
 	
 	[mf_ratio_500,line_act_500]= ... 
 	                  momfor(reac2_bot,x_bot,a,m_dis);
 	                  
	[T_mom,T_ten,T_dif]=T_sym(a,reac2_bot,x_bot,Tnd);
end
							        
%%%%%%%%%%%% End of Constraint   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%