Closed luchete80 closed 3 years ago
In 2 different python version, results differ.
Last eqn (DOF 16) (which is DOF 0 of the element)
Kt[0][0] [[-0.00111111 0. 0. 0. -0.00055556 0.
- -0.00027778 0. 0. 0. -0.00055556 0. 0. 0. ] [ 0. 0. 0. 0. 0. 0.
- 0.
- ] [ 0. 0. -0.00111111 0. 0. 0. -0.00055556 0. 0. 0. -0.00027778 0.
- -0.00055556 0. ] [ 0. 0. 0. 0. 0. 0.
- 0.
- ]
See above, this DOF is DOF 3 of the element, DOF 16 global. ALL ROWS ARE ZERO AND R IS 0:
Rglob [ 0. 0. 0. 0. -0.0025 0. 0. 0. 0.
- -0.0025 0. 0. 0. ]
SO THIS EQN IS NULL
(row) (col) 0 0 vnrow(row)vncol(col) 12.0 12.0 (row) (col) 0 1 vnrow(row)vncol(col) 12.0 13.0 (row) (col) 0 2 vnrow(row)vncol(col) 12.0 14.0 (row) (col) 0 3 vnrow(row)vncol(col) 12.0 15.0 (row) (col) 0 4 vnrow(row)vncol(col) 12.0 8.0 (row) (col) 0 5 vnrow(row)vncol(col) 12.0 9.0 (row) (col) 0 6 vnrow(row)vncol(col) 12.0 10.0 (row) (col) 0 7 vnrow(row)vncol(col) 12.0 11.0 (row) (col) 0 8 vnrow(row)vncol(col) 12.0 0.0 (row) (col) 0 9 vnrow(row)vncol(col) 12.0 1.0 (row) (col) 0 10 vnrow(row)vncol(col) 12.0 2.0 (row) (col) 0 11 vnrow(row)vncol(col) 12.0 3.0 (row) (col) 0 12 vnrow(row)vncol(col) 12.0 4.0 (row) (col) 0 13 vnrow(row)vncol(col) 12.0 5.0 (row) (col) 0 14 vnrow(row)vncol(col) 12.0 6.0 (row) (col) 0 15 vnrow(row)vncol(col) 12.0 7.0 (row) (col) 1 0 vnrow(row)vncol(col) 13.0 12.0 (row) (col) 1 1 vnrow(row)vncol(col) 13.0 13.0 (row) (col) 1 2 vnrow(row)vncol(col) 13.0 14.0 (row) (col) 1 3 vnrow(row)vncol(col) 13.0 15.0 (row) (col) 1 4 vnrow(row)vncol(col) 13.0 8.0 (row) (col) 1 5 vnrow(row)vncol(col) 13.0 9.0 (row) (col) 1 6 vnrow(row)vncol(col) 13.0 10.0 (row) (col) 1 7 vnrow(row)vncol(col) 13.0 11.0 (row) (col) 1 8 vnrow(row)vncol(col) 13.0 0.0 (row) (col) 1 9 vnrow(row)vncol(col) 13.0 1.0 (row) (col) 1 10 vnrow(row)vncol(col) 13.0 2.0 (row) (col) 1 11 vnrow(row)vncol(col) 13.0 3.0 (row) (col) 1 12 vnrow(row)vncol(col) 13.0 4.0 (row) (col) 1 13 vnrow(row)vncol(col) 13.0 5.0 (row) (col) 1 14 vnrow(row)vncol(col) 13.0 6.0 (row) (col) 1 15 vnrow(row)vncol(col) 13.0 7.0 (row) (col) 2 0 vnrow(row)vncol(col) 14.0 12.0 (row) (col) 2 1 vnrow(row)vncol(col) 14.0 13.0 (row) (col) 2 2 vnrow(row)vncol(col) 14.0 14.0 (row) (col) 2 3 vnrow(row)vncol(col) 14.0 15.0 (row) (col) 2 4 vnrow(row)vncol(col) 14.0 8.0 (row) (col) 2 5 vnrow(row)vncol(col) 14.0 9.0 (row) (col) 2 6 vnrow(row)vncol(col) 14.0 10.0 (row) (col) 2 7 vnrow(row)vncol(col) 14.0 11.0 (row) (col) 2 8 vnrow(row)vncol(col) 14.0 0.0 (row) (col) 2 9 vnrow(row)vncol(col) 14.0 1.0 (row) (col) 2 10 vnrow(row)vncol(col) 14.0 2.0 (row) (col) 2 11 vnrow(row)vncol(col) 14.0 3.0 (row) (col) 2 12 vnrow(row)vncol(col) 14.0 4.0 (row) (col) 2 13 vnrow(row)vncol(col) 14.0 5.0 (row) (col) 2 14 vnrow(row)vncol(col) 14.0 6.0 (row) (col) 2 15 vnrow(row)vncol(col) 14.0 7.0 (row) (col) 3 0 vnrow(row)vncol(col) 15.0 12.0 (row) (col) 3 1 vnrow(row)vncol(col) 15.0 13.0 (row) (col) 3 2 vnrow(row)vncol(col) 15.0 14.0 (row) (col) 3 3 vnrow(row)vncol(col) 15.0 15.0 (row) (col) 3 4 vnrow(row)vncol(col) 15.0 8.0 (row) (col) 3 5 vnrow(row)vncol(col) 15.0 9.0 (row) (col) 3 6 vnrow(row)vncol(col) 15.0 10.0 (row) (col) 3 7 vnrow(row)vncol(col) 15.0 11.0 (row) (col) 3 8 vnrow(row)vncol(col) 15.0 0.0 (row) (col) 3 9 vnrow(row)vncol(col) 15.0 1.0 (row) (col) 3 10 vnrow(row)vncol(col) 15.0 2.0 (row) (col) 3 11 vnrow(row)vncol(col) 15.0 3.0 (row) (col) 3 12 vnrow(row)vncol(col) 15.0 4.0 (row) (col) 3 13 vnrow(row)vncol(col) 15.0 5.0 (row) (col) 3 14 vnrow(row)vncol(col) 15.0 6.0 (row) (col) 3 15 vnrow(row)vncol(col) 15.0 7.0 (row) (col) 4 0 vnrow(row)vncol(col) 8.0 12.0 (row) (col) 4 1 vnrow(row)vncol(col) 8.0 13.0 (row) (col) 4 2 vnrow(row)vncol(col) 8.0 14.0 (row) (col) 4 3 vnrow(row)vncol(col) 8.0 15.0 (row) (col) 4 4 vnrow(row)vncol(col) 8.0 8.0 (row) (col) 4 5 vnrow(row)vncol(col) 8.0 9.0 (row) (col) 4 6 vnrow(row)vncol(col) 8.0 10.0 (row) (col) 4 7 vnrow(row)vncol(col) 8.0 11.0 (row) (col) 4 8 vnrow(row)vncol(col) 8.0 0.0 (row) (col) 4 9 vnrow(row)vncol(col) 8.0 1.0 (row) (col) 4 10 vnrow(row)vncol(col) 8.0 2.0 (row) (col) 4 11 vnrow(row)vncol(col) 8.0 3.0 (row) (col) 4 12 vnrow(row)vncol(col) 8.0 4.0 (row) (col) 4 13 vnrow(row)vncol(col) 8.0 5.0 (row) (col) 4 14 vnrow(row)vncol(col) 8.0 6.0 (row) (col) 4 15 vnrow(row)vncol(col) 8.0 7.0 (row) (col) 5 0 vnrow(row)vncol(col) 9.0 12.0 (row) (col) 5 1 vnrow(row)vncol(col) 9.0 13.0 (row) (col) 5 2 vnrow(row)vncol(col) 9.0 14.0 (row) (col) 5 3 vnrow(row)vncol(col) 9.0 15.0 (row) (col) 5 4 vnrow(row)vncol(col) 9.0 8.0 (row) (col) 5 5 vnrow(row)vncol(col) 9.0 9.0 (row) (col) 5 6 vnrow(row)vncol(col) 9.0 10.0 (row) (col) 5 7 vnrow(row)vncol(col) 9.0 11.0 (row) (col) 5 8 vnrow(row)vncol(col) 9.0 0.0 (row) (col) 5 9 vnrow(row)vncol(col) 9.0 1.0 (row) (col) 5 10 vnrow(row)vncol(col) 9.0 2.0 (row) (col) 5 11 vnrow(row)vncol(col) 9.0 3.0 (row) (col) 5 12 vnrow(row)vncol(col) 9.0 4.0 (row) (col) 5 13 vnrow(row)vncol(col) 9.0 5.0 (row) (col) 5 14 vnrow(row)vncol(col) 9.0 6.0 (row) (col) 5 15 vnrow(row)vncol(col) 9.0 7.0 (row) (col) 6 0 vnrow(row)vncol(col) 10.0 12.0 (row) (col) 6 1 vnrow(row)vncol(col) 10.0 13.0 (row) (col) 6 2 vnrow(row)vncol(col) 10.0 14.0 (row) (col) 6 3 vnrow(row)vncol(col) 10.0 15.0 (row) (col) 6 4 vnrow(row)vncol(col) 10.0 8.0 (row) (col) 6 5 vnrow(row)vncol(col) 10.0 9.0 (row) (col) 6 6 vnrow(row)vncol(col) 10.0 10.0 (row) (col) 6 7 vnrow(row)vncol(col) 10.0 11.0 (row) (col) 6 8 vnrow(row)vncol(col) 10.0 0.0 (row) (col) 6 9 vnrow(row)vncol(col) 10.0 1.0 (row) (col) 6 10 vnrow(row)vncol(col) 10.0 2.0 (row) (col) 6 11 vnrow(row)vncol(col) 10.0 3.0 (row) (col) 6 12 vnrow(row)vncol(col) 10.0 4.0 (row) (col) 6 13 vnrow(row)vncol(col) 10.0 5.0 (row) (col) 6 14 vnrow(row)vncol(col) 10.0 6.0 (row) (col) 6 15 vnrow(row)vncol(col) 10.0 7.0 (row) (col) 7 0 vnrow(row)vncol(col) 11.0 12.0 (row) (col) 7 1 vnrow(row)vncol(col) 11.0 13.0 (row) (col) 7 2 vnrow(row)vncol(col) 11.0 14.0 (row) (col) 7 3 vnrow(row)vncol(col) 11.0 15.0 (row) (col) 7 4 vnrow(row)vncol(col) 11.0 8.0 (row) (col) 7 5 vnrow(row)vncol(col) 11.0 9.0 (row) (col) 7 6 vnrow(row)vncol(col) 11.0 10.0 (row) (col) 7 7 vnrow(row)vncol(col) 11.0 11.0 (row) (col) 7 8 vnrow(row)vncol(col) 11.0 0.0 (row) (col) 7 9 vnrow(row)vncol(col) 11.0 1.0 (row) (col) 7 10 vnrow(row)vncol(col) 11.0 2.0 (row) (col) 7 11 vnrow(row)vncol(col) 11.0 3.0 (row) (col) 7 12 vnrow(row)vncol(col) 11.0 4.0 (row) (col) 7 13 vnrow(row)vncol(col) 11.0 5.0 (row) (col) 7 14 vnrow(row)vncol(col) 11.0 6.0 (row) (col) 7 15 vnrow(row)vncol(col) 11.0 7.0 (row) (col) 8 0 vnrow(row)vncol(col) 0.0 12.0 (row) (col) 8 1 vnrow(row)vncol(col) 0.0 13.0 (row) (col) 8 2 vnrow(row)vncol(col) 0.0 14.0 (row) (col) 8 3 vnrow(row)vncol(col) 0.0 15.0 (row) (col) 8 4 vnrow(row)vncol(col) 0.0 8.0 (row) (col) 8 5 vnrow(row)vncol(col) 0.0 9.0 (row) (col) 8 6 vnrow(row)vncol(col) 0.0 10.0 (row) (col) 8 7 vnrow(row)vncol(col) 0.0 11.0 (row) (col) 8 8 vnrow(row)vncol(col) 0.0 0.0 (row) (col) 8 9 vnrow(row)vncol(col) 0.0 1.0 (row) (col) 8 10 vnrow(row)vncol(col) 0.0 2.0 (row) (col) 8 11 vnrow(row)vncol(col) 0.0 3.0 (row) (col) 8 12 vnrow(row)vncol(col) 0.0 4.0 (row) (col) 8 13 vnrow(row)vncol(col) 0.0 5.0 (row) (col) 8 14 vnrow(row)vncol(col) 0.0 6.0 (row) (col) 8 15 vnrow(row)vncol(col) 0.0 7.0 (row) (col) 9 0 vnrow(row)vncol(col) 1.0 12.0 (row) (col) 9 1 vnrow(row)vncol(col) 1.0 13.0 (row) (col) 9 2 vnrow(row)vncol(col) 1.0 14.0 (row) (col) 9 3 vnrow(row)vncol(col) 1.0 15.0 (row) (col) 9 4 vnrow(row)vncol(col) 1.0 8.0 (row) (col) 9 5 vnrow(row)vncol(col) 1.0 9.0 (row) (col) 9 6 vnrow(row)vncol(col) 1.0 10.0 (row) (col) 9 7 vnrow(row)vncol(col) 1.0 11.0 (row) (col) 9 8 vnrow(row)vncol(col) 1.0 0.0 (row) (col) 9 9 vnrow(row)vncol(col) 1.0 1.0 (row) (col) 9 10 vnrow(row)vncol(col) 1.0 2.0 (row) (col) 9 11 vnrow(row)vncol(col) 1.0 3.0 (row) (col) 9 12 vnrow(row)vncol(col) 1.0 4.0 (row) (col) 9 13 vnrow(row)vncol(col) 1.0 5.0 (row) (col) 9 14 vnrow(row)vncol(col) 1.0 6.0 (row) (col) 9 15 vnrow(row)vncol(col) 1.0 7.0 (row) (col) 10 0 vnrow(row)vncol(col) 2.0 12.0 (row) (col) 10 1 vnrow(row)vncol(col) 2.0 13.0 (row) (col) 10 2 vnrow(row)vncol(col) 2.0 14.0 (row) (col) 10 3 vnrow(row)vncol(col) 2.0 15.0 (row) (col) 10 4 vnrow(row)vncol(col) 2.0 8.0 (row) (col) 10 5 vnrow(row)vncol(col) 2.0 9.0 (row) (col) 10 6 vnrow(row)vncol(col) 2.0 10.0 (row) (col) 10 7 vnrow(row)vncol(col) 2.0 11.0 (row) (col) 10 8 vnrow(row)vncol(col) 2.0 0.0 (row) (col) 10 9 vnrow(row)vncol(col) 2.0 1.0 (row) (col) 10 10 vnrow(row)vncol(col) 2.0 2.0 (row) (col) 10 11 vnrow(row)vncol(col) 2.0 3.0 (row) (col) 10 12 vnrow(row)vncol(col) 2.0 4.0 (row) (col) 10 13 vnrow(row)vncol(col) 2.0 5.0 (row) (col) 10 14 vnrow(row)vncol(col) 2.0 6.0 (row) (col) 10 15 vnrow(row)vncol(col) 2.0 7.0 (row) (col) 11 0 vnrow(row)vncol(col) 3.0 12.0 (row) (col) 11 1 vnrow(row)vncol(col) 3.0 13.0 (row) (col) 11 2 vnrow(row)vncol(col) 3.0 14.0 (row) (col) 11 3 vnrow(row)vncol(col) 3.0 15.0 (row) (col) 11 4 vnrow(row)vncol(col) 3.0 8.0 (row) (col) 11 5 vnrow(row)vncol(col) 3.0 9.0 (row) (col) 11 6 vnrow(row)vncol(col) 3.0 10.0 (row) (col) 11 7 vnrow(row)vncol(col) 3.0 11.0 (row) (col) 11 8 vnrow(row)vncol(col) 3.0 0.0 (row) (col) 11 9 vnrow(row)vncol(col) 3.0 1.0 (row) (col) 11 10 vnrow(row)vncol(col) 3.0 2.0 (row) (col) 11 11 vnrow(row)vncol(col) 3.0 3.0 (row) (col) 11 12 vnrow(row)vncol(col) 3.0 4.0 (row) (col) 11 13 vnrow(row)vncol(col) 3.0 5.0 (row) (col) 11 14 vnrow(row)vncol(col) 3.0 6.0 (row) (col) 11 15 vnrow(row)vncol(col) 3.0 7.0 (row) (col) 12 0 vnrow(row)vncol(col) 4.0 12.0 (row) (col) 12 1 vnrow(row)vncol(col) 4.0 13.0 (row) (col) 12 2 vnrow(row)vncol(col) 4.0 14.0 (row) (col) 12 3 vnrow(row)vncol(col) 4.0 15.0 (row) (col) 12 4 vnrow(row)vncol(col) 4.0 8.0 (row) (col) 12 5 vnrow(row)vncol(col) 4.0 9.0 (row) (col) 12 6 vnrow(row)vncol(col) 4.0 10.0 (row) (col) 12 7 vnrow(row)vncol(col) 4.0 11.0 (row) (col) 12 8 vnrow(row)vncol(col) 4.0 0.0 (row) (col) 12 9 vnrow(row)vncol(col) 4.0 1.0 (row) (col) 12 10 vnrow(row)vncol(col) 4.0 2.0 (row) (col) 12 11 vnrow(row)vncol(col) 4.0 3.0 (row) (col) 12 12 vnrow(row)vncol(col) 4.0 4.0 (row) (col) 12 13 vnrow(row)vncol(col) 4.0 5.0 (row) (col) 12 14 vnrow(row)vncol(col) 4.0 6.0 (row) (col) 12 15 vnrow(row)vncol(col) 4.0 7.0 (row) (col) 13 0 vnrow(row)vncol(col) 5.0 12.0 (row) (col) 13 1 vnrow(row)vncol(col) 5.0 13.0 (row) (col) 13 2 vnrow(row)vncol(col) 5.0 14.0 (row) (col) 13 3 vnrow(row)vncol(col) 5.0 15.0 (row) (col) 13 4 vnrow(row)vncol(col) 5.0 8.0 (row) (col) 13 5 vnrow(row)vncol(col) 5.0 9.0 (row) (col) 13 6 vnrow(row)vncol(col) 5.0 10.0 (row) (col) 13 7 vnrow(row)vncol(col) 5.0 11.0 (row) (col) 13 8 vnrow(row)vncol(col) 5.0 0.0 (row) (col) 13 9 vnrow(row)vncol(col) 5.0 1.0 (row) (col) 13 10 vnrow(row)vncol(col) 5.0 2.0 (row) (col) 13 11 vnrow(row)vncol(col) 5.0 3.0 (row) (col) 13 12 vnrow(row)vncol(col) 5.0 4.0 (row) (col) 13 13 vnrow(row)vncol(col) 5.0 5.0 (row) (col) 13 14 vnrow(row)vncol(col) 5.0 6.0 (row) (col) 13 15 vnrow(row)vncol(col) 5.0 7.0 (row) (col) 14 0 vnrow(row)vncol(col) 6.0 12.0 (row) (col) 14 1 vnrow(row)vncol(col) 6.0 13.0 (row) (col) 14 2 vnrow(row)vncol(col) 6.0 14.0 (row) (col) 14 3 vnrow(row)vncol(col) 6.0 15.0 (row) (col) 14 4 vnrow(row)vncol(col) 6.0 8.0 (row) (col) 14 5 vnrow(row)vncol(col) 6.0 9.0 (row) (col) 14 6 vnrow(row)vncol(col) 6.0 10.0 (row) (col) 14 7 vnrow(row)vncol(col) 6.0 11.0 (row) (col) 14 8 vnrow(row)vncol(col) 6.0 0.0 (row) (col) 14 9 vnrow(row)vncol(col) 6.0 1.0 (row) (col) 14 10 vnrow(row)vncol(col) 6.0 2.0 (row) (col) 14 11 vnrow(row)vncol(col) 6.0 3.0 (row) (col) 14 12 vnrow(row)vncol(col) 6.0 4.0 (row) (col) 14 13 vnrow(row)vncol(col) 6.0 5.0 (row) (col) 14 14 vnrow(row)vncol(col) 6.0 6.0 (row) (col) 14 15 vnrow(row)vncol(col) 6.0 7.0 (row) (col) 15 0 vnrow(row)vncol(col) 7.0 12.0 (row) (col) 15 1 vnrow(row)vncol(col) 7.0 13.0 (row) (col) 15 2 vnrow(row)vncol(col) 7.0 14.0 (row) (col) 15 3 vnrow(row)vncol(col) 7.0 15.0 (row) (col) 15 4 vnrow(row)vncol(col) 7.0 8.0 (row) (col) 15 5 vnrow(row)vncol(col) 7.0 9.0 (row) (col) 15 6 vnrow(row)vncol(col) 7.0 10.0 (row) (col) 15 7 vnrow(row)vncol(col) 7.0 11.0 (row) (col) 15 8 vnrow(row)vncol(col) 7.0 0.0 (row) (col) 15 9 vnrow(row)vncol(col) 7.0 1.0 (row) (col) 15 10 vnrow(row)vncol(col) 7.0 2.0 (row) (col) 15 11 vnrow(row)vncol(col) 7.0 3.0 (row) (col) 15 12 vnrow(row)vncol(col) 7.0 4.0 (row) (col) 15 13 vnrow(row)vncol(col) 7.0 5.0 (row) (col) 15 14 vnrow(row)vncol(col) 7.0 6.0 (row) (col) 15 15 vnrow(row)vncol(col) 7.0 7.0
Error, BsigF and B4i are null and Bs not:
B4i [[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]] i 2 B4i [[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]] i 3 B4i [[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0]]] BsigF [[[0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0]]
[[0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0] [0 0]]] Bs [[ 0.21132487 -0.21132487 -0.78867513 0.78867513] [ 0.21132487 0.78867513 -0.78867513 -0.21132487]]
Bv [[ 0.21132487 0. -0.21132487 0. -0.78867513 0. 0.78867513 0. ] [ 0.21132487 0. 0.78867513 0. -0.78867513 0. -0.21132487 0. ] [ 0. 0.21132487 0. -0.21132487 0. -0.78867513
- 0.78867513] [ 0. 0.21132487 0. 0.78867513 0. -0.78867513
- -0.21132487]]
Created for both methods, picard 965af04 and NR 422c59f