Como calcular a força de cisalhamento e os momentos de flexão e torção em uma asa?

Eu quero calcular as cargas (força de cisalhamento, flexão e momentos de torção). Supondo que eu tenha a distribuição de força ao longo de todas as superfícies da aeronave, meu ponto de partida é considerar a asa (e outras superfícies cujas cargas eu quero calcular) se comportar como um cantilever. Então, eu tenho estas fórmulas:

$$ S_z = \ int_z ^ {b / 2} F \ dy $$

$$ M = \ int_z ^ {b / 2} S_z \ dy $$

Onde

Sz = força de cisalhamento no ponto z; Mz = momento de flexão no ponto z; F = Distribuição do elevador;
b = Envergadura;

e o sistema de coordenadas é o seguinte:

Estas fórmulas estão corretas e como posso calcular o momento de torção?