Теоретические основы приближенного расчета нелинейных систем подрессоривания тяговых и транспортных гусеничных машин. Савочкин В.А - 40 стр.

                                                 40
     ωI = vgi – частоты среза возмущения, с-1 (частота среза ωI зависит от
скорости v движения ГМ и дорожных частот g, м-1, которые зависят от профиля
местности).
        9.2. Полагая, что жесткость каждой подвески равна жесткости подвески
при статистическом положении катка, а коэффициент сопротивления амортиза-
тора равен среднему значению коэффициента сопротивления при прямом и
обратном ходах катка при нулевой скорости перемещения катка, вычисляется,
используя выражение (76), передаточная функция относительного хода катка,
        9.3. Найденное значение передаточной функции (76) используется для
 определения спектральной плотности относительного хода ГМ:

                          Gλj( ω ) = Gу( ω ) / Wλj (p) /2                     (99)

       9.4. По полученному выражению (99) находятся средние квадратические
значения перемещения и скорости перемещения катка с помощью соотношений:
                                                            1
                                        1   ∝             2
                             σ λj =  ∫ Gλj (ω )dω                         (100)
                                    π 0           
                                                            1
                                        1 ∝               2
                              σ λ&j   =  ∫ ω 2 Gλj (ω )dω                 (101)
                                        π 0               

     9.5. Найденные средние квадратические значения            и    используются
для определения эквивалентных параметров каждой подвески в первом прибли-
жении по формулам (90 ... 92). При этом в первом приближении математического
ожидания хода катка задается равным его статистическому перемещению.
       9.6. С учетом значений σ λ j ,σ& λ j и сj, rj решаются уравнения (57) отно-
сительно mλj.
       9.7. Используя параметры mλj(1), сj(1), rj(1) первого приближения, снова
проводятся все вычисления, указанные в пунктах 2-6 для той же скорости движе-
ния ГМ. В результате находится следующее приближение для эквивалентных
параметров подвески. Вычисления проводятся до тех пор, пока эквивалентные
параметры не будут вычислены с достаточной степенью точности.
     Найденные значения эквивалентных параметров каждой подвески подстав-
ляются в выражения для передаточных функций (65...67), после чего определя-
ются спектральные плотности вертикальных Gz(ω) продольно-угловых Gϕ(ω) и
поперечно-угловых колебаний корпуса Gϒ(ω) ГМ, используя соотношения:

                                                                2
                                        Gz( ω ) = Gу(ω) Wz ( p ) ;          (102)


                                                                2
                                        Gϕ( ω ) = Gу(ω) Wϕ ( p) ;           (103)