Fix the rotation of aero force vector to body frame (#3)

zuorenchen · web-flow · commit ed2fe8b82321 · 2026-03-05T01:26:17.000Z
diff --git a/rocketpy/rocket/aero_surface/generic_surface.py b/rocketpy/rocket/aero_surface/generic_surface.py
@@ -292,14 +292,14 @@ def compute_forces_and_moments(
         rotation_matrix = Matrix(
             [
                 [1, 0, 0],
-                [0, math.cos(alpha), -math.sin(alpha)],
-                [0, math.sin(alpha), math.cos(alpha)],
+                [0, math.cos(alpha), math.sin(alpha)],
+                [0, -math.sin(alpha), math.cos(alpha)],
             ]
         ) @ Matrix(
             [
-                [math.cos(beta), 0, -math.sin(beta)],
+                [math.cos(beta), 0, math.sin(beta)],
                 [0, 1, 0],
-                [math.sin(beta), 0, math.cos(beta)],
+                [-math.sin(beta), 0, math.cos(beta)],
             ]
         )
         R1, R2, R3 = rotation_matrix @ Vector([side, -lift, -drag])

Original file line number	Diff line number	Diff line change
`@@ -292,14 +292,14 @@ def compute_forces_and_moments(`
`292`	`292`	`rotation_matrix = Matrix(`
`293`	`293`	`[`
`294`	`294`	`[1, 0, 0],`
`295`		`- [0, math.cos(alpha), -math.sin(alpha)],`
`296`		`- [0, math.sin(alpha), math.cos(alpha)],`
	`295`	`+ [0, math.cos(alpha), math.sin(alpha)],`
	`296`	`+ [0, -math.sin(alpha), math.cos(alpha)],`
`297`	`297`	`]`
`298`	`298`	`) @ Matrix(`
`299`	`299`	`[`
`300`		`- [math.cos(beta), 0, -math.sin(beta)],`
	`300`	`+ [math.cos(beta), 0, math.sin(beta)],`
`301`	`301`	`[0, 1, 0],`
`302`		`- [math.sin(beta), 0, math.cos(beta)],`
	`302`	`+ [-math.sin(beta), 0, math.cos(beta)],`
`303`	`303`	`]`
`304`	`304`	`)`
`305`	`305`	`R1, R2, R3 = rotation_matrix @ Vector([side, -lift, -drag])`