diff options
| author | Michael Abed <michaelabed@gmail.com> | 2012-04-09 22:28:34 -0400 |
|---|---|---|
| committer | Michael Abed <michaelabed@gmail.com> | 2012-04-09 22:28:34 -0400 |
| commit | 2b61d521c51300575c981a48d256c87cbcf6e478 (patch) | |
| tree | e7eb8c6fe6e6b206bcb1285b26c63f87415a7ef8 /truss_solver.m | |
| parent | f950899cf0133d864abe4fece0fac469e28e8023 (diff) | |
| download | ek301-prelim-2b61d521c51300575c981a48d256c87cbcf6e478.tar.gz ek301-prelim-2b61d521c51300575c981a48d256c87cbcf6e478.tar.bz2 ek301-prelim-2b61d521c51300575c981a48d256c87cbcf6e478.zip | |
update thinsg
Diffstat (limited to 'truss_solver.m')
| -rwxr-xr-x | truss_solver.m | 11 |
1 files changed, 9 insertions, 2 deletions
diff --git a/truss_solver.m b/truss_solver.m index 50086f3..81e3cc8 100755 --- a/truss_solver.m +++ b/truss_solver.m @@ -13,8 +13,9 @@ D = memberDistances(C, X, Y); cost = computeTrussCost(C,D); weakest = weakestMember(D,forces); -ml = maxLoad(forces,L,D,weakest); +ml = maxLoad(forces,L,D,weakest,false); +printHeader(); for i = 1:length(Ln) fprintf('Load: %.3f N\n', abs(Ln(i))); @@ -31,4 +32,10 @@ end fprintf('Cost of truss: $%d\n', round(cost)); 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 +fprintf('First member to break is member %d\n', weakest); +fprintf('The buclking force is %.3f N\n' , bucklingStrength(D(weakest))); + +weak = maxLoad(forces,L,D,weakest,true); + +fprintf('Worst possible buckling strength is %.3f N\n', weak); +fprintf('Theoretical worst load/cost ratio: %.4f\n', weak/cost);
\ No newline at end of file |