Deformação de duas vigas perpendiculares com propriedades parabólicas

Quero calcular a deformação da seguinte estrutura:

Numerei o feixe vertical 1 e o feixe horzontal 2. A altura do feixe 2 é dada por: $h(x)=h\sqrt{1-x/L_2}$

$F$$Fa$$0$$Fa$$U=\int_0^L\frac{F(x)^2}{2AE}dx+\int_0^L\frac{M(x)^2}{2EI}dx$

$$U_1 = \int_0^{L_1}\frac{F^2}{2bhE}dx+\int_0^{L_1}\frac{F^2a^2}{2E\frac{1}{12}bh^3}dx = \frac{F^2L_1}{2bhE} + \frac{F^2a^2L_1}{2E\frac{1}{12}bh^3}$$

$$U_2 = \int_0^a\frac{F^2x^2}{2E\frac{1}{12}b(h\sqrt{1-x/L_2})^3}dx = \frac{4F^2L_2^2}{Ebh^3}\left(\frac{\sqrt{1-a/L_2}(a^2+4aL_2-8L_2^2)}{a-L_2}-8L_2\right)$$

$b$

$F$

$$\delta=\frac{\partial U}{\partial F}= \frac{F}{E}\left(\frac{L_1}{bh} + \frac{a^2L_1}{\frac{1}{12}bh^3} + \frac{8L_2^2}{bh^3}\left(\frac{\sqrt{1-a/L_2}(a^2+4aL_2-8L_2^2)}{a-L_2}-8L_2\right)\right)$$

$F=1500, E=113.76*10^9, b=0.005, h=0.007, L_1=0.05, L_2=0.04$$a=0.032$

Eu tentei isso sem a propriedade parabólica e isso gera resultados corretos. Essa equação é:

$$\frac{F}{E}\left(\frac{L_1}{A_1}+\frac{L_2^2L_1}{I_1}+\frac{L_2^3}{3I_3}\right)$$

Então, acho que há algo errado com minha equação parabólica, mas não consigo descobrir o que é.

$h(h/2)$$h$$\frac{8L_2^2}{bh^3}$

$U_2$

$$U_2 = \int_0^a\frac{F^2(x-a)^2}{2E\frac{1}{12}b(h\sqrt{1-x/L_2})^3} dx$$

$$\delta= \frac{F}{E}\left(\frac{L_1}{bh} + \frac{a^2L_1}{\frac{1}{12}bh^3} + \frac{8L_2}{bh^3}\left(-3a^2+12aL_2+8L_2 \sqrt{1-\frac{a}{L_2}}(L_2-a)-8L_2^2\right)\right)$$