Homec4science
Diffusion MSM (command_line) Browse

msm/ShoulderModel_HHT/GeneralizedForces_6DOF.m
73754adfee53command_line

Last Change
Raw File

GeneralizedForces_6DOF.m
View Options

Q(1,1) = cos(q2)*(F_final(1)*Dir_inT(1, 1)*Rho_inT(3, 1) - F_final(1)*Dir_inT(3, 1)*Rho_inT(1, 1) + F_final(2)*Dir_inT(1, 2)*Rho_inT(3, 2) - F_final(2)*Dir_inT(3, 2)*Rho_inT(1, 2) + F_final(3)*Dir_inT(1, 3)*Rho_inT(3, 3) - F_final(3)*Dir_inT(3, 3)*Rho_inT(1, 3) + F_final(4)*Dir_inT(1, 4)*Rho_inT(3, 4) - F_final(4)*Dir_inT(3, 4)*Rho_inT(1, 4) + F_final(5)*Dir_inT(1, 5)*Rho_inT(3, 5) - F_final(5)*Dir_inT(3, 5)*Rho_inT(1, 5) + F_final(6)*Dir_inT(1, 6)*Rho_inT(3, 6) - F_final(6)*Dir_inT(3, 6)*Rho_inT(1, 6) + F_contact(1)*Dir_contact(1, 1)*Rho_contact(3, 1) - F_contact(1)*Dir_contact(3, 1)*Rho_contact(1, 1) + F_contact(2)*Dir_contact(1, 2)*Rho_contact(3, 2) - F_contact(2)*Dir_contact(3, 2)*Rho_contact(1, 2) + F_contact(3)*Dir_contact(1, 3)*Rho_contact(3, 3) - F_contact(3)*Dir_contact(3, 3)*Rho_contact(1, 3) + F_contact(4)*Dir_contact(1, 4)*Rho_contact(3, 4) - F_contact(4)*Dir_contact(3, 4)*Rho_contact(1, 4) + F_contact(5)*Dir_contact(1, 5)*Rho_contact(3, 5) - F_contact(5)*Dir_contact(3, 5)*Rho_contact(1, 5) + F_contact(6)*Dir_contact(1, 6)*Rho_contact(3, 6) - F_contact(6)*Dir_contact(3, 6)*Rho_contact(1, 6) + F_contact(7)*Dir_contact(1, 7)*Rho_contact(3, 7) - F_contact(7)*Dir_contact(3, 7)*Rho_contact(1, 7) + F_contact(8)*Dir_contact(1, 8)*Rho_contact(3, 8) - F_contact(8)*Dir_contact(3, 8)*Rho_contact(1, 8) + F_contact(9)*Dir_contact(1, 9)*Rho_contact(3, 9) - F_contact(9)*Dir_contact(3, 9)*Rho_contact(1, 9) + F_contact(10)*Dir_contact(1, 10)*Rho_contact(3, 10) - F_contact(10)*Dir_contact(3, 10)*Rho_contact(1, 10) + F_contact(11)*Dir_contact(1, 11)*Rho_contact(3, 11) - F_contact(11)*Dir_contact(3, 11)*Rho_contact(1, 11) + F_contact(12)*Dir_contact(1, 12)*Rho_contact(3, 12) - F_contact(12)*Dir_contact(3, 12)*Rho_contact(1, 12) + F_contact(13)*Dir_contact(1, 13)*Rho_contact(3, 13) - F_contact(13)*Dir_contact(3, 13)*Rho_contact(1, 13) + F_contact(14)*Dir_contact(1, 14)*Rho_contact(3, 14) - F_contact(14)*Dir_contact(3, 14)*Rho_contact(1, 14) + F_contact(15)*Dir_contact(1, 15)*Rho_contact(3, 15) - F_contact(15)*Dir_contact(3, 15)*Rho_contact(1, 15) + F_contact(16)*Dir_contact(1, 16)*Rho_contact(3, 16) - F_contact(16)*Dir_contact(3, 16)*Rho_contact(1, 16) + F_contact(17)*Dir_contact(1, 17)*Rho_contact(3, 17) - F_contact(17)*Dir_contact(3, 17)*Rho_contact(1, 17) + F_contact(18)*Dir_contact(1, 18)*Rho_contact(3, 18) - F_contact(18)*Dir_contact(3, 18)*Rho_contact(1, 18) + F_contact(19)*Dir_contact(1, 19)*Rho_contact(3, 19) - F_contact(19)*Dir_contact(3, 19)*Rho_contact(1, 19) + F_contact(20)*Dir_contact(1, 20)*Rho_contact(3, 20) - F_contact(20)*Dir_contact(3, 20)*Rho_contact(1, 20)) - cos(q3)*sin(q2)*(F_final(1)*Dir_inT(1, 1)*Rho_inT(2, 1) - F_final(1)*Dir_inT(2, 1)*Rho_inT(1, 1) + F_final(2)*Dir_inT(1, 2)*Rho_inT(2, 2) - F_final(2)*Dir_inT(2, 2)*Rho_inT(1, 2) + F_final(3)*Dir_inT(1, 3)*Rho_inT(2, 3) - F_final(3)*Dir_inT(2, 3)*Rho_inT(1, 3) + F_final(4)*Dir_inT(1, 4)*Rho_inT(2, 4) - F_final(4)*Dir_inT(2, 4)*Rho_inT(1, 4) + F_final(5)*Dir_inT(1, 5)*Rho_inT(2, 5) - F_final(5)*Dir_inT(2, 5)*Rho_inT(1, 5) + F_final(6)*Dir_inT(1, 6)*Rho_inT(2, 6) - F_final(6)*Dir_inT(2, 6)*Rho_inT(1, 6) + F_contact(1)*Dir_contact(1, 1)*Rho_contact(2, 1) - F_contact(1)*Dir_contact(2, 1)*Rho_contact(1, 1) + F_contact(2)*Dir_contact(1, 2)*Rho_contact(2, 2) - F_contact(2)*Dir_contact(2, 2)*Rho_contact(1, 2) + F_contact(3)*Dir_contact(1, 3)*Rho_contact(2, 3) - F_contact(3)*Dir_contact(2, 3)*Rho_contact(1, 3) + F_contact(4)*Dir_contact(1, 4)*Rho_contact(2, 4) - F_contact(4)*Dir_contact(2, 4)*Rho_contact(1, 4) + F_contact(5)*Dir_contact(1, 5)*Rho_contact(2, 5) - F_contact(5)*Dir_contact(2, 5)*Rho_contact(1, 5) + F_contact(6)*Dir_contact(1, 6)*Rho_contact(2, 6) - F_contact(6)*Dir_contact(2, 6)*Rho_contact(1, 6) + F_contact(7)*Dir_contact(1, 7)*Rho_contact(2, 7) - F_contact(7)*Dir_contact(2, 7)*Rho_contact(1, 7) + F_contact(8)*Dir_contact(1, 8)*Rho_contact(2, 8) - F_contact(8)*Dir_contact(2, 8)*Rho_contact(1, 8) + F_contact(9)*Dir_contact(1, 9)*Rho_contact(2, 9) - F_contact(9)*Dir_contact(2, 9)*Rho_contact(1, 9) + F_contact(10)*Dir_contact(1, 10)*Rho_contact(2, 10) - F_contact(10)*Dir_contact(2, 10)*Rho_contact(1, 10) + F_contact(11)*Dir_contact(1, 11)*Rho_contact(2, 11) - F_contact(11)*Dir_contact(2, 11)*Rho_contact(1, 11) + F_contact(12)*Dir_contact(1, 12)*Rho_contact(2, 12) - F_contact(12)*Dir_contact(2, 12)*Rho_contact(1, 12) + F_contact(13)*Dir_contact(1, 13)*Rho_contact(2, 13) - F_contact(13)*Dir_contact(2, 13)*Rho_contact(1, 13) + F_contact(14)*Dir_contact(1, 14)*Rho_contact(2, 14) - F_contact(14)*Dir_contact(2, 14)*Rho_contact(1, 14) + F_contact(15)*Dir_contact(1, 15)*Rho_contact(2, 15) - F_contact(15)*Dir_contact(2, 15)*Rho_contact(1, 15) + F_contact(16)*Dir_contact(1, 16)*Rho_contact(2, 16) - F_contact(16)*Dir_contact(2, 16)*Rho_contact(1, 16) + F_contact(17)*Dir_contact(1, 17)*Rho_contact(2, 17) - F_contact(17)*Dir_contact(2, 17)*Rho_contact(1, 17) + F_contact(18)*Dir_contact(1, 18)*Rho_contact(2, 18) - F_contact(18)*Dir_contact(2, 18)*Rho_contact(1, 18) + F_contact(19)*Dir_contact(1, 19)*Rho_contact(2, 19) - F_contact(19)*Dir_contact(2, 19)*Rho_contact(1, 19) + F_contact(20)*Dir_contact(1, 20)*Rho_contact(2, 20) - F_contact(20)*Dir_contact(2, 20)*Rho_contact(1, 20)) - sin(q2)*sin(q3)*(F_final(1)*Dir_inT(2, 1)*Rho_inT(3, 1) - F_final(1)*Dir_inT(3, 1)*Rho_inT(2, 1) + F_final(2)*Dir_inT(2, 2)*Rho_inT(3, 2) - F_final(2)*Dir_inT(3, 2)*Rho_inT(2, 2) + F_final(3)*Dir_inT(2, 3)*Rho_inT(3, 3) - F_final(3)*Dir_inT(3, 3)*Rho_inT(2, 3) + F_final(4)*Dir_inT(2, 4)*Rho_inT(3, 4) - F_final(4)*Dir_inT(3, 4)*Rho_inT(2, 4) + F_final(5)*Dir_inT(2, 5)*Rho_inT(3, 5) - F_final(5)*Dir_inT(3, 5)*Rho_inT(2, 5) + F_final(6)*Dir_inT(2, 6)*Rho_inT(3, 6) - F_final(6)*Dir_inT(3, 6)*Rho_inT(2, 6) + F_contact(1)*Dir_contact(2, 1)*Rho_contact(3, 1) - F_contact(1)*Dir_contact(3, 1)*Rho_contact(2, 1) + F_contact(2)*Dir_contact(2, 2)*Rho_contact(3, 2) - F_contact(2)*Dir_contact(3, 2)*Rho_contact(2, 2) + F_contact(3)*Dir_contact(2, 3)*Rho_contact(3, 3) - F_contact(3)*Dir_contact(3, 3)*Rho_contact(2, 3) + F_contact(4)*Dir_contact(2, 4)*Rho_contact(3, 4) - F_contact(4)*Dir_contact(3, 4)*Rho_contact(2, 4) + F_contact(5)*Dir_contact(2, 5)*Rho_contact(3, 5) - F_contact(5)*Dir_contact(3, 5)*Rho_contact(2, 5) + F_contact(6)*Dir_contact(2, 6)*Rho_contact(3, 6) - F_contact(6)*Dir_contact(3, 6)*Rho_contact(2, 6) + F_contact(7)*Dir_contact(2, 7)*Rho_contact(3, 7) - F_contact(7)*Dir_contact(3, 7)*Rho_contact(2, 7) + F_contact(8)*Dir_contact(2, 8)*Rho_contact(3, 8) - F_contact(8)*Dir_contact(3, 8)*Rho_contact(2, 8) + F_contact(9)*Dir_contact(2, 9)*Rho_contact(3, 9) - F_contact(9)*Dir_contact(3, 9)*Rho_contact(2, 9) + F_contact(10)*Dir_contact(2, 10)*Rho_contact(3, 10) - F_contact(10)*Dir_contact(3, 10)*Rho_contact(2, 10) + F_contact(11)*Dir_contact(2, 11)*Rho_contact(3, 11) - F_contact(11)*Dir_contact(3, 11)*Rho_contact(2, 11) + F_contact(12)*Dir_contact(2, 12)*Rho_contact(3, 12) - F_contact(12)*Dir_contact(3, 12)*Rho_contact(2, 12) + F_contact(13)*Dir_contact(2, 13)*Rho_contact(3, 13) - F_contact(13)*Dir_contact(3, 13)*Rho_contact(2, 13) + F_contact(14)*Dir_contact(2, 14)*Rho_contact(3, 14) - F_contact(14)*Dir_contact(3, 14)*Rho_contact(2, 14) + F_contact(15)*Dir_contact(2, 15)*Rho_contact(3, 15) - F_contact(15)*Dir_contact(3, 15)*Rho_contact(2, 15) + F_contact(16)*Dir_contact(2, 16)*Rho_contact(3, 16) - F_contact(16)*Dir_contact(3, 16)*Rho_contact(2, 16) + F_contact(17)*Dir_contact(2, 17)*Rho_contact(3, 17) - F_contact(17)*Dir_contact(3, 17)*Rho_contact(2, 17) + F_contact(18)*Dir_contact(2, 18)*Rho_contact(3, 18) - F_contact(18)*Dir_contact(3, 18)*Rho_contact(2, 18) + F_contact(19)*Dir_contact(2, 19)*Rho_contact(3, 19) - F_contact(19)*Dir_contact(3, 19)*Rho_contact(2, 19) + F_contact(20)*Dir_contact(2, 20)*Rho_contact(3, 20) - F_contact(20)*Dir_contact(3, 20)*Rho_contact(2, 20));
Q(2,1) = sin(q3)*(F_final(1)*Dir_inT(1, 1)*Rho_inT(2, 1) - F_final(1)*Dir_inT(2, 1)*Rho_inT(1, 1) + F_final(2)*Dir_inT(1, 2)*Rho_inT(2, 2) - F_final(2)*Dir_inT(2, 2)*Rho_inT(1, 2) + F_final(3)*Dir_inT(1, 3)*Rho_inT(2, 3) - F_final(3)*Dir_inT(2, 3)*Rho_inT(1, 3) + F_final(4)*Dir_inT(1, 4)*Rho_inT(2, 4) - F_final(4)*Dir_inT(2, 4)*Rho_inT(1, 4) + F_final(5)*Dir_inT(1, 5)*Rho_inT(2, 5) - F_final(5)*Dir_inT(2, 5)*Rho_inT(1, 5) + F_final(6)*Dir_inT(1, 6)*Rho_inT(2, 6) - F_final(6)*Dir_inT(2, 6)*Rho_inT(1, 6) + F_contact(1)*Dir_contact(1, 1)*Rho_contact(2, 1) - F_contact(1)*Dir_contact(2, 1)*Rho_contact(1, 1) + F_contact(2)*Dir_contact(1, 2)*Rho_contact(2, 2) - F_contact(2)*Dir_contact(2, 2)*Rho_contact(1, 2) + F_contact(3)*Dir_contact(1, 3)*Rho_contact(2, 3) - F_contact(3)*Dir_contact(2, 3)*Rho_contact(1, 3) + F_contact(4)*Dir_contact(1, 4)*Rho_contact(2, 4) - F_contact(4)*Dir_contact(2, 4)*Rho_contact(1, 4) + F_contact(5)*Dir_contact(1, 5)*Rho_contact(2, 5) - F_contact(5)*Dir_contact(2, 5)*Rho_contact(1, 5) + F_contact(6)*Dir_contact(1, 6)*Rho_contact(2, 6) - F_contact(6)*Dir_contact(2, 6)*Rho_contact(1, 6) + F_contact(7)*Dir_contact(1, 7)*Rho_contact(2, 7) - F_contact(7)*Dir_contact(2, 7)*Rho_contact(1, 7) + F_contact(8)*Dir_contact(1, 8)*Rho_contact(2, 8) - F_contact(8)*Dir_contact(2, 8)*Rho_contact(1, 8) + F_contact(9)*Dir_contact(1, 9)*Rho_contact(2, 9) - F_contact(9)*Dir_contact(2, 9)*Rho_contact(1, 9) + F_contact(10)*Dir_contact(1, 10)*Rho_contact(2, 10) - F_contact(10)*Dir_contact(2, 10)*Rho_contact(1, 10) + F_contact(11)*Dir_contact(1, 11)*Rho_contact(2, 11) - F_contact(11)*Dir_contact(2, 11)*Rho_contact(1, 11) + F_contact(12)*Dir_contact(1, 12)*Rho_contact(2, 12) - F_contact(12)*Dir_contact(2, 12)*Rho_contact(1, 12) + F_contact(13)*Dir_contact(1, 13)*Rho_contact(2, 13) - F_contact(13)*Dir_contact(2, 13)*Rho_contact(1, 13) + F_contact(14)*Dir_contact(1, 14)*Rho_contact(2, 14) - F_contact(14)*Dir_contact(2, 14)*Rho_contact(1, 14) + F_contact(15)*Dir_contact(1, 15)*Rho_contact(2, 15) - F_contact(15)*Dir_contact(2, 15)*Rho_contact(1, 15) + F_contact(16)*Dir_contact(1, 16)*Rho_contact(2, 16) - F_contact(16)*Dir_contact(2, 16)*Rho_contact(1, 16) + F_contact(17)*Dir_contact(1, 17)*Rho_contact(2, 17) - F_contact(17)*Dir_contact(2, 17)*Rho_contact(1, 17) + F_contact(18)*Dir_contact(1, 18)*Rho_contact(2, 18) - F_contact(18)*Dir_contact(2, 18)*Rho_contact(1, 18) + F_contact(19)*Dir_contact(1, 19)*Rho_contact(2, 19) - F_contact(19)*Dir_contact(2, 19)*Rho_contact(1, 19) + F_contact(20)*Dir_contact(1, 20)*Rho_contact(2, 20) - F_contact(20)*Dir_contact(2, 20)*Rho_contact(1, 20)) - cos(q3)*(F_final(1)*Dir_inT(2, 1)*Rho_inT(3, 1) - F_final(1)*Dir_inT(3, 1)*Rho_inT(2, 1) + F_final(2)*Dir_inT(2, 2)*Rho_inT(3, 2) - F_final(2)*Dir_inT(3, 2)*Rho_inT(2, 2) + F_final(3)*Dir_inT(2, 3)*Rho_inT(3, 3) - F_final(3)*Dir_inT(3, 3)*Rho_inT(2, 3) + F_final(4)*Dir_inT(2, 4)*Rho_inT(3, 4) - F_final(4)*Dir_inT(3, 4)*Rho_inT(2, 4) + F_final(5)*Dir_inT(2, 5)*Rho_inT(3, 5) - F_final(5)*Dir_inT(3, 5)*Rho_inT(2, 5) + F_final(6)*Dir_inT(2, 6)*Rho_inT(3, 6) - F_final(6)*Dir_inT(3, 6)*Rho_inT(2, 6) + F_contact(1)*Dir_contact(2, 1)*Rho_contact(3, 1) - F_contact(1)*Dir_contact(3, 1)*Rho_contact(2, 1) + F_contact(2)*Dir_contact(2, 2)*Rho_contact(3, 2) - F_contact(2)*Dir_contact(3, 2)*Rho_contact(2, 2) + F_contact(3)*Dir_contact(2, 3)*Rho_contact(3, 3) - F_contact(3)*Dir_contact(3, 3)*Rho_contact(2, 3) + F_contact(4)*Dir_contact(2, 4)*Rho_contact(3, 4) - F_contact(4)*Dir_contact(3, 4)*Rho_contact(2, 4) + F_contact(5)*Dir_contact(2, 5)*Rho_contact(3, 5) - F_contact(5)*Dir_contact(3, 5)*Rho_contact(2, 5) + F_contact(6)*Dir_contact(2, 6)*Rho_contact(3, 6) - F_contact(6)*Dir_contact(3, 6)*Rho_contact(2, 6) + F_contact(7)*Dir_contact(2, 7)*Rho_contact(3, 7) - F_contact(7)*Dir_contact(3, 7)*Rho_contact(2, 7) + F_contact(8)*Dir_contact(2, 8)*Rho_contact(3, 8) - F_contact(8)*Dir_contact(3, 8)*Rho_contact(2, 8) + F_contact(9)*Dir_contact(2, 9)*Rho_contact(3, 9) - F_contact(9)*Dir_contact(3, 9)*Rho_contact(2, 9) + F_contact(10)*Dir_contact(2, 10)*Rho_contact(3, 10) - F_contact(10)*Dir_contact(3, 10)*Rho_contact(2, 10) + F_contact(11)*Dir_contact(2, 11)*Rho_contact(3, 11) - F_contact(11)*Dir_contact(3, 11)*Rho_contact(2, 11) + F_contact(12)*Dir_contact(2, 12)*Rho_contact(3, 12) - F_contact(12)*Dir_contact(3, 12)*Rho_contact(2, 12) + F_contact(13)*Dir_contact(2, 13)*Rho_contact(3, 13) - F_contact(13)*Dir_contact(3, 13)*Rho_contact(2, 13) + F_contact(14)*Dir_contact(2, 14)*Rho_contact(3, 14) - F_contact(14)*Dir_contact(3, 14)*Rho_contact(2, 14) + F_contact(15)*Dir_contact(2, 15)*Rho_contact(3, 15) - F_contact(15)*Dir_contact(3, 15)*Rho_contact(2, 15) + F_contact(16)*Dir_contact(2, 16)*Rho_contact(3, 16) - F_contact(16)*Dir_contact(3, 16)*Rho_contact(2, 16) + F_contact(17)*Dir_contact(2, 17)*Rho_contact(3, 17) - F_contact(17)*Dir_contact(3, 17)*Rho_contact(2, 17) + F_contact(18)*Dir_contact(2, 18)*Rho_contact(3, 18) - F_contact(18)*Dir_contact(3, 18)*Rho_contact(2, 18) + F_contact(19)*Dir_contact(2, 19)*Rho_contact(3, 19) - F_contact(19)*Dir_contact(3, 19)*Rho_contact(2, 19) + F_contact(20)*Dir_contact(2, 20)*Rho_contact(3, 20) - F_contact(20)*Dir_contact(3, 20)*Rho_contact(2, 20));
Q(3,1) = F_final(1)*Dir_inT(1, 1)*Rho_inT(3, 1) - F_final(1)*Dir_inT(3, 1)*Rho_inT(1, 1) + F_final(2)*Dir_inT(1, 2)*Rho_inT(3, 2) - F_final(2)*Dir_inT(3, 2)*Rho_inT(1, 2) + F_final(3)*Dir_inT(1, 3)*Rho_inT(3, 3) - F_final(3)*Dir_inT(3, 3)*Rho_inT(1, 3) + F_final(4)*Dir_inT(1, 4)*Rho_inT(3, 4) - F_final(4)*Dir_inT(3, 4)*Rho_inT(1, 4) + F_final(5)*Dir_inT(1, 5)*Rho_inT(3, 5) - F_final(5)*Dir_inT(3, 5)*Rho_inT(1, 5) + F_final(6)*Dir_inT(1, 6)*Rho_inT(3, 6) - F_final(6)*Dir_inT(3, 6)*Rho_inT(1, 6) + F_contact(1)*Dir_contact(1, 1)*Rho_contact(3, 1) - F_contact(1)*Dir_contact(3, 1)*Rho_contact(1, 1) + F_contact(2)*Dir_contact(1, 2)*Rho_contact(3, 2) - F_contact(2)*Dir_contact(3, 2)*Rho_contact(1, 2) + F_contact(3)*Dir_contact(1, 3)*Rho_contact(3, 3) - F_contact(3)*Dir_contact(3, 3)*Rho_contact(1, 3) + F_contact(4)*Dir_contact(1, 4)*Rho_contact(3, 4) - F_contact(4)*Dir_contact(3, 4)*Rho_contact(1, 4) + F_contact(5)*Dir_contact(1, 5)*Rho_contact(3, 5) - F_contact(5)*Dir_contact(3, 5)*Rho_contact(1, 5) + F_contact(6)*Dir_contact(1, 6)*Rho_contact(3, 6) - F_contact(6)*Dir_contact(3, 6)*Rho_contact(1, 6) + F_contact(7)*Dir_contact(1, 7)*Rho_contact(3, 7) - F_contact(7)*Dir_contact(3, 7)*Rho_contact(1, 7) + F_contact(8)*Dir_contact(1, 8)*Rho_contact(3, 8) - F_contact(8)*Dir_contact(3, 8)*Rho_contact(1, 8) + F_contact(9)*Dir_contact(1, 9)*Rho_contact(3, 9) - F_contact(9)*Dir_contact(3, 9)*Rho_contact(1, 9) + F_contact(10)*Dir_contact(1, 10)*Rho_contact(3, 10) - F_contact(10)*Dir_contact(3, 10)*Rho_contact(1, 10) + F_contact(11)*Dir_contact(1, 11)*Rho_contact(3, 11) - F_contact(11)*Dir_contact(3, 11)*Rho_contact(1, 11) + F_contact(12)*Dir_contact(1, 12)*Rho_contact(3, 12) - F_contact(12)*Dir_contact(3, 12)*Rho_contact(1, 12) + F_contact(13)*Dir_contact(1, 13)*Rho_contact(3, 13) - F_contact(13)*Dir_contact(3, 13)*Rho_contact(1, 13) + F_contact(14)*Dir_contact(1, 14)*Rho_contact(3, 14) - F_contact(14)*Dir_contact(3, 14)*Rho_contact(1, 14) + F_contact(15)*Dir_contact(1, 15)*Rho_contact(3, 15) - F_contact(15)*Dir_contact(3, 15)*Rho_contact(1, 15) + F_contact(16)*Dir_contact(1, 16)*Rho_contact(3, 16) - F_contact(16)*Dir_contact(3, 16)*Rho_contact(1, 16) + F_contact(17)*Dir_contact(1, 17)*Rho_contact(3, 17) - F_contact(17)*Dir_contact(3, 17)*Rho_contact(1, 17) + F_contact(18)*Dir_contact(1, 18)*Rho_contact(3, 18) - F_contact(18)*Dir_contact(3, 18)*Rho_contact(1, 18) + F_contact(19)*Dir_contact(1, 19)*Rho_contact(3, 19) - F_contact(19)*Dir_contact(3, 19)*Rho_contact(1, 19) + F_contact(20)*Dir_contact(1, 20)*Rho_contact(3, 20) - F_contact(20)*Dir_contact(3, 20)*Rho_contact(1, 20);
Q(4,1) = F_contact(1)*Dir_contact(1, 1) + F_contact(2)*Dir_contact(1, 2) + F_contact(3)*Dir_contact(1, 3) + F_contact(4)*Dir_contact(1, 4) + F_contact(5)*Dir_contact(1, 5) + F_contact(6)*Dir_contact(1, 6) + F_contact(7)*Dir_contact(1, 7) + F_contact(8)*Dir_contact(1, 8) + F_contact(9)*Dir_contact(1, 9) + F_contact(10)*Dir_contact(1, 10) + F_contact(11)*Dir_contact(1, 11) + F_contact(12)*Dir_contact(1, 12) + F_contact(13)*Dir_contact(1, 13) + F_contact(14)*Dir_contact(1, 14) + F_contact(15)*Dir_contact(1, 15) + F_contact(16)*Dir_contact(1, 16) + F_contact(17)*Dir_contact(1, 17) + F_contact(18)*Dir_contact(1, 18) + F_contact(19)*Dir_contact(1, 19) + F_contact(20)*Dir_contact(1, 20) + F_final(1)*Dir_inT(1, 1) + F_final(2)*Dir_inT(1, 2) + F_final(3)*Dir_inT(1, 3) + F_final(4)*Dir_inT(1, 4) + F_final(5)*Dir_inT(1, 5) + F_final(6)*Dir_inT(1, 6);
Q(5,1) = F_contact(1)*Dir_contact(2, 1) + F_contact(2)*Dir_contact(2, 2) + F_contact(3)*Dir_contact(2, 3) + F_contact(4)*Dir_contact(2, 4) + F_contact(5)*Dir_contact(2, 5) + F_contact(6)*Dir_contact(2, 6) + F_contact(7)*Dir_contact(2, 7) + F_contact(8)*Dir_contact(2, 8) + F_contact(9)*Dir_contact(2, 9) + F_contact(10)*Dir_contact(2, 10) + F_contact(11)*Dir_contact(2, 11) + F_contact(12)*Dir_contact(2, 12) + F_contact(13)*Dir_contact(2, 13) + F_contact(14)*Dir_contact(2, 14) + F_contact(15)*Dir_contact(2, 15) + F_contact(16)*Dir_contact(2, 16) + F_contact(17)*Dir_contact(2, 17) + F_contact(18)*Dir_contact(2, 18) + F_contact(19)*Dir_contact(2, 19) + F_contact(20)*Dir_contact(2, 20) + F_final(1)*Dir_inT(2, 1) + F_final(2)*Dir_inT(2, 2) + F_final(3)*Dir_inT(2, 3) + F_final(4)*Dir_inT(2, 4) + F_final(5)*Dir_inT(2, 5) + F_final(6)*Dir_inT(2, 6);
Q(6,1) = F_contact(1)*Dir_contact(3, 1) + F_contact(2)*Dir_contact(3, 2) + F_contact(3)*Dir_contact(3, 3) + F_contact(4)*Dir_contact(3, 4) + F_contact(5)*Dir_contact(3, 5) + F_contact(6)*Dir_contact(3, 6) + F_contact(7)*Dir_contact(3, 7) + F_contact(8)*Dir_contact(3, 8) + F_contact(9)*Dir_contact(3, 9) + F_contact(10)*Dir_contact(3, 10) + F_contact(11)*Dir_contact(3, 11) + F_contact(12)*Dir_contact(3, 12) + F_contact(13)*Dir_contact(3, 13) + F_contact(14)*Dir_contact(3, 14) + F_contact(15)*Dir_contact(3, 15) + F_contact(16)*Dir_contact(3, 16) + F_contact(17)*Dir_contact(3, 17) + F_contact(18)*Dir_contact(3, 18) + F_contact(19)*Dir_contact(3, 19) + F_contact(20)*Dir_contact(3, 20) + F_final(1)*Dir_inT(3, 1) + F_final(2)*Dir_inT(3, 2) + F_final(3)*Dir_inT(3, 3) + F_final(4)*Dir_inT(3, 4) + F_final(5)*Dir_inT(3, 5) + F_final(6)*Dir_inT(3, 6);
c4science · Help