pressure

(./pressure.m)

Help text

 dof = 3 for a truss, 6 for a frame

Cross-Reference Information

forces         ./forces.m
normal_loads   ./normal_loads.m
plot_truss     ./plot_truss.m
solve_truss    ./solve_truss.m

Listing of ./pressure.m

nodes = [0,0,0;
        10,0,0;
        10,10,0;
        20,10,0;
        20,00,0;
        30,00,0];    

A = 8;                     % in^2
E = 29000;                 % ksi
g = 11250;
j = .001;
iss = .001;
itt = .001;

prop = [A,E,0,iss,itt,j,g];
members = [
        1,2,prop;
        2,3,prop;
        3,4,prop;
        4,5,prop;
        5,6,prop];        %node 1, node 2,area, E
        
active_nodes = [2:5];

loads = zeros(4,7);        %node, fi,fj,fk,mi,mj,mk
loads(1:4,1:4) =  normal_loads(nodes,active_nodes, members, -1)

% dof = 3 for a truss, 6 for a frame

dof = 6;


displacements = solve_truss(nodes,members,active_nodes,loads,dof)
plot_truss(nodes,members,displacements,1);
member_forces = forces (displacements, members,nodes)
setup_rotate;
view(2);