williamjohnstone/uni-raymers
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Changed numbers in stress calcs (also testing commit signing)

Browse source
This commit is contained in:
William Johnstone 2024-11-09 20:35:39 +00:00
parent 7985dfa1d3
commit 10d0ff5f60
1 changed files with 21 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

25
Raymers Code/raymers.py
View file

@ -115,16 +115,33 @@ print(f"Bending Moment at of 1 Wing from root: {bending_moment_wing:.4f} Nm\n")
# Max Stress Calcs # Max Stress Calcs
B = 4 * 0.001 # Width in m
D = 4 * 0.001 # Height or depth in m
# 4mm Spars # 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") 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") 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 # 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") 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") 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")