```python
from sympy import *
from sympy.physics.units import kg, m, s

Newton = kg*m/s**2

a, F, q = var("a, F, q")

sub_list=[
    ( a, 1 *m         ),
    ( F, 1 *Newton    ),
    ( q, 1 *Newton/m  ),
    ]

Ah, Av, MA, Bv, Gh, Gv \
    = var("Ah, Av, MA, Bv, Gh, Gv")

R = 2*a*q

# equilibrium conditions:
eq1 = Eq(  Gh )
eq2 = Eq(  - Bv + F + Gv + R  )
eq3 = Eq(  2*a*F + a*(-Bv + R)  )
eq4 = Eq(  Ah + Gh)
eq5 = Eq(  Av + Gv)
eq6 = Eq(  Gv*a - MA )

unks = [Ah, Av, MA, Bv, Gh, Gv]
eqs = [eq1, eq2, eq3, eq4, eq5, eq6]

sol = solve(eqs, unks)
pprint(sol)

Ah, Av, MA, Bv, Gh, Gv \
    = sol[Ah], sol[Av], sol[MA], sol[Bv], sol[Gh], sol[Gv]

pprint("\nAh, Av, Bv, Gh, Gv, MA in SI units:")
for s in [Ah, Av, Bv, Gh, Gv]:
    tmp = s
    tmp = tmp.subs(sub_list)
    tmp /= Newton
    pprint(tmp)

tmp = MA
tmp = tmp.subs(sub_list)
tmp /= Newton*m
pprint(tmp)

pprint("\nAh, Av, Bv, MA in SI units:")
Ah, Av, MA, Bv= var("Ah, Av, MA, Bv")
eq1 = Eq(  2*F*a + a*(-Bv + R) )
eq2 = Eq(  Ah )
eq3 = Eq(  Av + Bv - R - F )
eq4 = Eq(  2*a*Bv -2*a*R -3*a*F - MA )

unks = [Ah, Av, MA, Bv]
eqs = [eq1, eq2, eq3, eq4]

sol = solve(eqs, unks)
pprint(sol)

# {Ah: 0, Av: -F, Bv: 2⋅F + 2⋅a⋅q, Gh: 0, Gv: F, MA: F⋅a}
#
# Ah, Av, Bv, Gh, Gv, MA in SI units:
# 0
# -1
# 4
# 0
# 1
# 1
#
# Ah, Av, Bv, MA in SI units:
# {Ah: 0, Av: -F, Bv: 2⋅F + 2⋅a⋅q, MA: F⋅a}

```

2a323a8b496c51c6f91b5a3bea06441a6797dd98
ipynb
Jupyter Notebook
ipynb/TM_1/5_AL/1.5.K.ipynb
kassbohm/tm-snippets
5e0621ba2470116e54643b740d1b68b9f28bff12
MIT
ipynb/TM_1/5_AL/1.5.K.ipynb
kassbohm/tm-snippets
5e0621ba2470116e54643b740d1b68b9f28bff12
MIT
ipynb/TM_1/5_AL/1.5.K.ipynb
kassbohm/tm-snippets
5e0621ba2470116e54643b740d1b68b9f28bff12
MIT
Qwen/Qwen-72B
1. YES
2. YES