heatk

(./heatk.m)

Function Synopsis

[kt] = heatK(Pi,Pj,Area,K)

Help text

function [kt] = heatK(Pi,Pj,Area,K);
 or function [kt] = heatK(member,nodes)
 This function computes the K tensor for a truss
 element, given the component parts.
 member = [node 1,node 2, A, K], nodes = node list, [p1;p2;...]

Cross-Reference Information

dot            ./dot.m

assemble_k     ./assemble_k.m

Listing of function ./heatk.m

function [kt] = heatK(Pi,Pj,Area,K)
%function [kt] = heatK(Pi,Pj,Area,K);
% or function [kt] = heatK(member,nodes)
% This function computes the K tensor for a truss
% element, given the component parts.
% member = [node 1,node 2, A, K], nodes = node list, [p1;p2;...]

if nargin == 2 
        member = Pi;
        nodes = Pj;
        Pi = nodes(member(1),:);
        Pj = nodes(member(2),:);
        Area = member(3);
    K = member(4);
end

nL = Pj' - Pi';
L = sqrt(dot(nL,nL));

kt = -K*Area/L;