2.4.3.A.sagews

from sympy.physics.units import *
from sympy import *
from sympy import N as Num

# pie is better than pi
pie = 2*pi

print("\n--- User input -----------------------")
symbolic = True
# symbolic = False

if symbolic:
    # quantity = symbol:
    alpha = var("alpha")
    a, b = var("a, b")
    F = var("F")
else:
    # quantity =  factor times unit:
    # pie/12 = 30 deg
    # pie = 2 pi
    alpha = pie/8

    a = 1 *m
    b = 2 *mm
    F = 1 *newton

A = 12*b*b
pprint("A / mm*mm")
pprint(A/mm/mm)

Iyy = 53*b**4/4

print("\n--- a: -------------------------------")

print("\n--- b: -------------------------------")

N_M = -F
M_M = 3*a*F

print("\n--- c: -------------------------------")

print("\n Stress:")
for z in [- 9*b/4, 7*b/4]:
    sigN = N_M / A
    sigB = M_M / Iyy * z
    sig = sigB + sigN
    pprint(sig)

    if not (symbolic):
        print("sig / MPa:")
        pprint(Num(sig/10**6/pascal, 4))

--- User input ----------------------- A / mm*mm 2 12000000*b ----------- 2 m --- a: ------------------------------- --- b: ------------------------------- --- c: ------------------------------- Stress: 27*F*a F - ------ - ----- 3 2 53*b 12*b 21*F*a F ------ - ----- 3 2 53*b 12*b