Function rc_alpha_y_sugano 

pub fn rc_alpha_y_sugano(
    pt: f64,
    a_over_d: f64,
    d_over_full: f64,
    n: f64,
) -> f64

RC 梁の曲げ降伏時剛性低下率 αy（菅野式）。

αy = (0.043 + 1.635·n·pt + 0.043·(a/D))·(d/D)²   (2.0 ≤ a/D ≤ 5.0)
     (−0.0836 + 0.159·(a/D))·(d/D)²              (1.0 ≤ a/D < 2.0)

• pt: 引張鉄筋比（小数）
• a_over_d: シアスパン比 a/D（a=l0/2）。適用範囲 [1.0, 5.0] にクランプする。
• d_over_full: 有効せい/全せい d/D
• n: ヤング係数比 Es/Ec

出典: 梅村魁『鉄筋コンクリート建物の動的耐震設計法』P.106-108（要・原典照合）。 トリリニア骨格の降伏点変形（θy=θe/αy）に用いる剛性低下率で、0〜1 に収まる想定。 負となる異常入力は 0 にクランプする（1 超は補正しない＝呼び出し側で扱う）。