forces

(./forces.m)

Function Synopsis

f = forces(nodes,members,displacements);

Help text

function f = forces(nodes,members,displacements);
change for frames

Cross-Reference Information

dot            ./dot.m
trussk         ./trussk.m

AeroTruss1     ./AeroTruss1.m
plot_truss     ./plot_truss.m
pressure       ./pressure.m
truss1         ./truss1.m
wheel2         ./wheel2.m
wheellu        ./wheellu.m

Listing of function ./forces.m

function f = forces(nodes,members,displacements);
%function f = forces(nodes,members,displacements);
%change for frames

[nr,nc] = size(nodes);
truss = 1:nc;

[num_mem,i] = size(members);
f = zeros(num_mem,1);
for i = 1:num_mem
    node_i = members(i,1);
    node_j = members(i,2);
    ui = displacements(node_i,truss)';
    uj = displacements(node_j,truss)';
    [kt,ael,n] = trussk(members(i,:),nodes);
    f(i) = -ael*dot(n,(ui-uj));
end