trussk

(./trussk.m)

Function Synopsis

[kt,ael,n] = trussK(Pi,Pj,Area,E)

Help text

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

Cross-Reference Information

dot            ./dot.m
tensor         ./tensor.m

assemble_k     ./assemble_k.m
assemble_truss_k./assemble_truss_k.m
forces         ./forces.m

Listing of function ./trussk.m

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

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

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

kt = (Area*E/L) * tensor(n,n);
ael = Area*E/L;