diff options
author | Michael Abed <michaelabed@gmail.com> | 2012-04-07 23:14:54 -0400 |
---|---|---|
committer | Michael Abed <michaelabed@gmail.com> | 2012-04-07 23:14:54 -0400 |
commit | 81117fb87ae3edb8b7bb31962f52b8c85af431fa (patch) | |
tree | e5d824ffec9bc86ba94532b3fb7ff3e4225ca09a /truss_solver.m | |
parent | 2ad700f049fc8b62c33eab73c641cf67e81effbc (diff) | |
download | ek301-prelim-81117fb87ae3edb8b7bb31962f52b8c85af431fa.tar.gz ek301-prelim-81117fb87ae3edb8b7bb31962f52b8c85af431fa.tar.bz2 ek301-prelim-81117fb87ae3edb8b7bb31962f52b8c85af431fa.zip |
update
Diffstat (limited to 'truss_solver.m')
-rwxr-xr-x | truss_solver.m | 13 |
1 files changed, 10 insertions, 3 deletions
diff --git a/truss_solver.m b/truss_solver.m index 4a7a1aa..50086f3 100755 --- a/truss_solver.m +++ b/truss_solver.m @@ -8,6 +8,13 @@ forces(abs(forces) < 0.00001) = 0; [j, m] = size(C); Ln = L(L~=0); + +D = memberDistances(C, X, Y); + +cost = computeTrussCost(C,D); +weakest = weakestMember(D,forces); +ml = maxLoad(forces,L,D,weakest); + for i = 1:length(Ln) fprintf('Load: %.3f N\n', abs(Ln(i))); @@ -21,7 +28,7 @@ for i = 1:length(forces)-m fprintf('r%d: %.3f\n', i, forces(m+i)); end -cost = computeTrussCost(C, X, Y); fprintf('Cost of truss: $%d\n', round(cost)); - -fprintf('Theoretical max load/cost ratio in N/$: %.4f\n', abs(sum(L))/cost);
\ No newline at end of file +fprintf('Theoretical max load is %.4f\n', ml); +fprintf('Theoretical max load/cost ratio in N/$: %.4f\n', ml/cost); +fprintf('First member to break is member %d\n', weakest);
\ No newline at end of file |