search

spacether / pycalculix

Python 3 library to automate and build finite element analysis (FEA) models in Calculix. Meshing uses Calculix or GMSH.
http://justinablack.com/pycalculix/
Apache License 2.0
170 stars 56 forks source link

Example `dam.py` fails to run #16

Closed yw5aj closed 7 years ago

yw5aj commented 7 years ago

Hi! Thank you so much for the great package.

I am running Anaconda 64 bit with Python 3.6 on my Windows 10 64 bit computer. I used pip install pycalculix to install pycalculix successfully.

When I cloned the entire repo and run python dam.py, everything else runs fine until this showed up:

Traceback (most recent call last): File "dam.py", line 81, in prob.solve() File "C:\Users\Shawn\Anaconda3\lib\site-packages\pycalculix\problem.py", line 285, in solve self.rfile.load() File "C:\Users\Shawn\Anaconda3\lib\site-packages\pycalculix\results_file.py", line 92, in load self.read_frd() # read nodal results File "C:\Users\Shawn\Anaconda3\lib\site-packages\pycalculix\results_file.py", line 1126, in read_frd self.set_time(self.__steps[0]) IndexError: list index out of range

And I checked the entire log and found:

STEP 1

Static analysis was selected

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 1

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 2

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 3

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 4

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 5

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 6

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 7

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 8

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 9

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 10

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 11

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

*ERROR in e_c3d: nonpositive jacobian determinant in element 12

Could you please give me a hint on where would you think have possibly gone wrong?

Thanks, Shawn

spacether commented 7 years ago

Please try uninstalling pycalculix and resinstalling it from this repository. Calculix changed how the frd file was formatted. The new version, 0.9.4 should work.

yw5aj commented 7 years ago

Thank you Justin! It seems that a new error occurred as I was trying to re-install pycalculix (https://github.com/spacether/pycalculix/issues/25). I will report back when I find something.