0001 function [sx,sy,sxy,s]=shear2d(md)
0002
0003
0004
0005
0006
0007
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);