flw2te.m

基于matlab的有限元程序源代码

function [Ke,fe]=flw2te(ex,ey,ep,D,eq)
% Ke=flw2te(ex,ey,ep,D)
% [Ke,fe]=flw2te(ex,ey,ep,D,eq)
%-------------------------------------------------------------
% PURPOSE
%  Compute element stiffness (conductivity) matrix for a 
%  triangular field element.
%
% INPUT:  ex = [x1 x2 x3]
%         ey = [y1 y2 y3]      element coordinates
%
%         ep = [t]  	       element thickness  	 
%                             
%         D = [kxx kxy;
%              kyx kyy]        constitutive matrix
%
%         eq                   heat supply per unit volume
%
% OUTPUT: Ke :  element 'stiffness' matrix (3 x 3)
%
%         fe :  element load vector (3 x 1)
%-------------------------------------------------------------

% LAST MODIFIED: K Persson   1995-08-23
% Copyright (c) 1991-94 by Division of Structural Mechanics and
%                          Department of Solid Mechanics.
%                          Lund Institute of Technology
%-------------------------------------------------------------
  t=ep(1);
  if nargin==4; eq=0; end 
%
  C=[ones(3,1) ex' ey'];   B=[0 1 0;
                              0 0 1 ]*inv(C);   A=1/2*det(C);
  Ke1=B'*D*B*t*A;
  
  fe1=eq*A*t/3*[1 1 1]'; 

  Ke=Ke1;  fe=fe1;
%--------------------------end--------------------------------