vault backup: 2025-02-27 16:40:30

补课/多体动力学/05_vectors.ipynb

@@ -75,13 +75,6 @@
    "source": [
     "C.orient_body_fixed(B, (q3, q4, 0), 'XZX')"
    ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": []
   }
  ],
  "metadata": {

补课/多体动力学/05_vectors.py

@@ -0,0 +1,30 @@
+import sympy as sm
+import sympy.physics.mechanics as me
+sm.init_printing()
+
+q1, q2, q3, q4 = sm.symbols('q1:5')
+l1, l2, l3, l4 = sm.symbols('l1:5')
+
+N, A, B, C = sm.symbols('N, A, B, C', cls=me.ReferenceFrame)
+
+A.orient_body_fixed(N, (q1, q2, 0), 'ZXZ')
+print(A.dcm(N))
+
+B.orient_axis(A, q3, A.x)
+
+C.orient_body_fixed(B, (q3, q4, 0), 'XZX')
+
+r_P1_P2 = l1 * A.z
+print(r_P1_P2)
+
+r_P2_P3 = l2 * B.z
+print(r_P2_P3)
+
+r_P3_P4 = l3 *C.z - l4 *C.y
+print(r_P3_P4)
+
+r_P1_P4 = r_P1_P2 + r_P2_P3 + r_P3_P4
+print(r_P1_P4)
+
+print(r_P1_P4.express(B))
+print(r_P1_P4.free_symbols(N))

补课/多体动力学/sympy_start.ipynb

@@ -289,7 +289,7 @@
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "Turtle",
+   "display_name": "blade",
    "language": "python",
    "name": "python3"
   },
@@ -303,7 +303,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.8"
+   "version": "3.10.16"
   }
  },
  "nbformat": 4,