Введем положительно ориентированную декартову систему координат

. Ось

смотрит на нас перпендикулярно плоскости рисунка, ось

вертикально вверх. Через

обозначим центр масс бруска,

. Через

обозначим угол наклона бруска к оси

.
Давление в жидкости распределено в соответствие с известной барометрической формулой

где

-- константа из стартового поста.
Со стороны жидкости на покоящийся брусок действует сила

, где

-- контур бруска;

-- внешняя единичная нормаль к контуру. И момент
![$\boldsymbol M_S(\phi,z)=-\int_L[\boldsymbol \rho,\boldsymbol n]pds.$ $\boldsymbol M_S(\phi,z)=-\int_L[\boldsymbol \rho,\boldsymbol n]pds.$](https://dxdy-03.korotkov.co.uk/f/2/e/0/2e0c35fc190561b1752f3085a4ed2a5682.png)
(

-- радиус вектор из центра бруска в точку на контуре)
Условия равновесия:

-- масса бруска. Из закона Архимеда известно, что

направлена вертикально. Поэтому уравнения (*) это два скалярных уравнения на

. Пусть

-- одно из решений этой системы.
Дальше можно предложить следующую наивную схему исследования устойчивости. Предположим, что силы сопротивления жикости пропорциональны квадратам скоростей соответствующих точек бруска, то мы получим следующую систему линейного приближения

Ну и для начала надо посмотреть на резуьтаты, которые даст исследование устойчивости системы состоящей только из первых двух уравнений.
Замечание. Во всех формулах используется "двумерное" давление

, его размерность H/м