resolve_truss

(./resolve_truss.m)

Function Synopsis

[u] = resolve_truss(in_sys_k,nodes,active_nodes,loads,dof);

Help text

  [u] = resolve_truss(in_sys_k,nodes,active_nodes,loads,dof)
 does the solving. Truss or frame.

Cross-Reference Information

assemble_loads ./assemble_loads.m

Listing of function ./resolve_truss.m

function [u] = resolve_truss(in_sys_k,nodes,active_nodes,loads,dof);
%  [u] = resolve_truss(in_sys_k,nodes,active_nodes,loads,dof)
% does the solving. Truss or frame.


num_act_nodes = max(size(active_nodes));  % number of active nodes
[num_loads,i] = size(loads);               % How many applied loads

red_load=zeros(dof*num_act_nodes,1);   % generates a reduced empty load
                                       % column vector
disp = zeros(dof*num_act_nodes,1);      % a displacement vector to solve for

[na,nc] = size(nodes);
u = zeros(na,dof);                     % the return displacements

for i = 1:num_loads
   red_load = assemble_loads(red_load,loads(i,:),active_nodes,dof);
end

disp = in_sys_k*red_load;

u(active_nodes,:) = reshape(disp,dof,num_act_nodes)';