shear2d

PURPOSE ^

SHEAR2D - computes 2d strain rate

SYNOPSIS ^

function [sx,sy,sxy,s]=shear2d(md)

DESCRIPTION ^

SHEAR2D - computes 2d strain rate

   This routine computes the strain rate of 2d models

   Usage:
      [sx,sy,sxy,s]=shear2d(md);

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [sx,sy,sxy,s]=shear2d(md)
0002 %SHEAR2D - computes 2d strain rate
0003 %
0004 %   This routine computes the strain rate of 2d models
0005 %
0006 %   Usage:
0007 %      [sx,sy,sxy,s]=shear2d(md);
0008 
0009 alpha=zeros(md.numberofelements,3); beta=zeros(md.numberofelements,3); 
0010 gamma=zeros(md.numberofelements,3); area=zeros(md.numberofelements,1);
0011 
0012 for n=1:md.numberofelements
0013     X=inv([md.x(md.elements(n,:)) md.y(md.elements(n,:)) ones(3,1)]);
0014     alpha(n,:)=X(1,:);
0015     beta(n,:)=X(2,:);
0016 end
0017 clear X;
0018 
0019 
0020 summation=[1;1;1];
0021 sx=(md.vx(md.elements).*alpha)*summation;
0022 uy=(md.vx(md.elements).*beta)*summation;
0023 vx=(md.vy(md.elements).*alpha)*summation;
0024 sy=(md.vy(md.elements).*beta)*summation;                        
0025 sxy=(uy+vx)/2;
0026 s=sqrt(sx.^2+sy.^2+sxy.^2+sx.*sy);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003