eforce.m

function  eq=eforce(ex,fx,Eq)
 
  % INPUT:
  %   ex - współrzędne wektorowe elementu [ex1,ex2]
  %   fx - długość odcinka na którym przyłożone jest obciążenie 
  %   Eq - wektor obciążenia w układzie globalnym [Eq_x,Eq_y]
  % OUTPUT:
  %   eq - wektor obciążenia w układzie lokalnym elementu [eq_x,eq_y]
  %
  % Piotr Pluciński 29.11.2007 
  %-----------------------------------------------------------
  %
  %    fx=3
  %  |_______|
  %   _______ 
  %  | | | | | Zy=10
  %  V V V V V
  %             _       
  %          / |        
  %         /  |        
  %   y    /   |
  %  ^    /    | ex2=4       eq=eforce([3,4],3,[0,-10]) 
  %  |   /     |
  %  |  /      |
  %  | /   x   |
  %   / -->    |_
  %   ________ 
  %  | ex1=3  |
 
 Le=sqrt(ex(1)^2+ex(2)^2);
 Eqe=Eq*fx/Le;
 c=ex(1)/Le;
 s=ex(2)/Le;
 T=[c -s; s c];
 eq=Eqe*T;