From 10d0ff5f60c63d7b862613d8690ab5def449bec7 Mon Sep 17 00:00:00 2001 From: William Johnstone Date: Sat, 9 Nov 2024 20:35:39 +0000 Subject: [PATCH] Changed numbers in stress calcs (also testing commit signing) --- Raymers Code/raymers.py | 25 +++++++++++++++++++++---- 1 file changed, 21 insertions(+), 4 deletions(-) diff --git a/Raymers Code/raymers.py b/Raymers Code/raymers.py index 154e011..40fabe4 100644 --- a/Raymers Code/raymers.py +++ b/Raymers Code/raymers.py @@ -115,16 +115,33 @@ print(f"Bending Moment at of 1 Wing from root: {bending_moment_wing:.4f} Nm\n") # Max Stress Calcs +B = 4 * 0.001 # Width in m +D = 4 * 0.001 # Height or depth in m + # 4mm Spars -area_moment_inertia = ((0.004) * (0.004**3))/12 +area_moment_inertia = ((B) * (D**3))/12 print(f"Second moment of inertia for 4mm: {area_moment_inertia} m^4") -max_stress = bending_moment_wing * ((2*(10**-3)) / area_moment_inertia) +max_stress = bending_moment_wing * ((D/2) / area_moment_inertia) print(f"Max Stress for 4mm spar: {max_stress:.4f} Pa or {max_stress/1000000:.4f} MPa") +B = 6 * 0.001 # Width in m +D = 6 * 0.001 # Height or depth in m + # 6mm Spars -area_moment_inertia = ((0.006) * (0.006**3))/12 +area_moment_inertia = ((B) * (D**3))/12 print(f"Second moment of inertia for 6mm: {area_moment_inertia} m^4") -max_stress = bending_moment_wing * ((3*(10**-3)) / area_moment_inertia) +max_stress = bending_moment_wing * ((D/2) / area_moment_inertia) print(f"Max Stress for 6mm spar: {max_stress:.4f} Pa or {max_stress/1000000:.4f} MPa") + + +B = 12 * 0.001 # Width in m +D = 2 * 0.001 # Height or depth in m + +# 6mm Spars +area_moment_inertia = ((B) * (D**3))/12 +print(f"Second moment of inertia for new test: {area_moment_inertia} m^4") + +max_stress = bending_moment_wing * ((D/2) / area_moment_inertia) +print(f"Max Stress for new test spar: {max_stress:.4f} Pa or {max_stress/1000000:.4f} MPa")