Уравнения равновесия и кривизны

[sviyav]

Здравствуйте,подскажите,пожалуйста,решаю задачу по изгибу тавра численно-заменяю после интегралы по формуле Симпсона ,для этого выписала ур.равновесия

$$\int\limits_{0}^{h_1}b_1\sigma dS+\int\limits_{h_1}^{h_1+h_2}b_2\sigma dS=0$$

,

$$\int\limits_{0}^{h_1}b_1\sigma(z+\delta) dS+\int\limits_{h_1}^{h_1+h_2}b_2\sigma(z+\delta) dS=0$$

кривизну $\kappa=M/EJ+1/J(\int\limits_{0}^{h_1} b_1\varepsilon(z+\delta)+\int\limits_{h_1}^{h_1+h_2} b_2\varepsilon(z+\delta))dz$ (из гип. плоских сечений)
и смещение нейтральной оси $\delta=1/\kappa h_1\int\limits_{0}^{h_1}\varepsilon dz+1/\kappa h_2\int\limits_{h_1}^{h_1+h_2} \varepsilon dz$.
НЕ уверена в правильности формул.Где ошибка.понять не могу. Решаю аналогично задаче балки(расчеты сошлись)