Skip to content

Commit 6e27d40

Browse files
committed
docs: sync the rotor-lever comment with the coned hub-radius formula
A stale inline comment still described the polar lever as r = hub_r + span*cos(cone), from before the hub radius was coned. Update it to r = (hub_r + span)*cos(cone) so the comment matches the implementation and the surrounding docstrings.
1 parent e38287b commit 6e27d40

1 file changed

Lines changed: 4 additions & 3 deletions

File tree

src/pybmodes/io/windio.py

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1206,9 +1206,10 @@ def _nac_geom(key: str) -> float:
12061206

12071207
# Rotor inertia from the spanwise blade mass (issue #130). In the shaft
12081208
# frame the rotor tensor is diag([I_polar, I_diam, I_diam]): the polar
1209-
# term ``I_polar = N_bl·∫dm·r²`` (radial ``r = hub_r + span·cos(cone)``),
1210-
# and the transverse (diametral) term. For a coned rotor the blade mass
1211-
# sits off the hub plane by ``a = span·sin(cone)``, so the rotor CM shifts
1209+
# term ``I_polar = N_bl·∫dm·r²`` (radial ``r = (hub_r + span)·cos(cone)``,
1210+
# the coned pitch-axis distance from the hub centre), and the transverse
1211+
# (diametral) term. For a coned rotor the blade mass sits off the hub
1212+
# plane by ``a = (hub_r + span)·sin(cone)``, so the rotor CM shifts
12121213
# along the shaft axis by ``offset = ∫dm·a / m_blades`` and the tensor is
12131214
# formed about that CM: ``I_diam = I_polar/2 + N_bl·∫dm·a² − m_blades·
12141215
# offset²`` (the axial second moment reduced to the CM by the parallel-

0 commit comments

Comments
 (0)