Hi,
I installed the binary for Linux, also installed vtk and numpy from my package manager but get the following ModuleNotFoundError

Traceback (most recent call last): File "/media/ihor/WORK/Calculix/ccx2paraview/ccx2paraview/__init__.py", line 21, in <module> File "/media/ihor/WORK/Calculix/ccx2paraview/ccx2paraview/writer.py", line 14, in <module writer> File "/home/ihor/.local/lib/python3.8/site-packages/vtk.py", line 30, in <module vtk> File "/usr/lib/python3.8/importlib/__init__.py", line 127, in import_module return _bootstrap._gcd_import(name[level:], package, level) ModuleNotFoundError: No module named 'vtkmodules'

I am on Linux Mint Cinnamon 20.2 and cannot work out what is missing.
Would appreciate any hint.

Hello. Thank you for your nice work. I notice that in the converted .vtu files, there is no real timesteps, only seriers like [0,1,2,3...].
However this can be solved by by apendding one .pvd file like this:
<VTKFile type="Collection" version="0.1" byte_order="LittleEndian">
<Collection>
<DataSet timestep="0.01" file="tank.01.vtu"/>
<DataSet timestep="0.02" file="tank.02.vtu"/>
<DataSet timestep="0.03" file="tank.03.vtu"/>
</Collection>
</VTKFile>
which 0.01,0.02,0.03 are real timesteps.
How about inplenting auto generation of this .pvd file into code?

Thanks a lot. I had already tried it before. But it doesn't work for me. On the other hand, version 2.3.1 contained in the (bin) folder of the zip file does work correctly. I guess it has to do with the fact that I don't have Python installed under Windows. Could it be that this version requires some component that I don't have?

Originally posted by @emilio-martin in #18 (comment)

Hello Ihor,

I have been trying to post-process a large .frd file (Almost ~7GB).
This is mainly due to very small time-steps (dt=1e-05) and a relatively fine mesh (t=0.1s).
I would also like to mention that this is only when I write

I tried to convert this file into .vtu.
But my laptop runs of RAM and it soon hangs or crashes.
Even with exiting from all other unnecessary applications, this happens.

INFO: Writing Nozzle_hex.166.vtu
INFO: Step 166, time 1.66e-03, U, 3 components, 89054 values
INFO: Writing Nozzle_hex.167.vtu
Killed

I think this is mainly because first the data is converted for all the time-steps and then it starts writing.

I tried understanding the code, but I am not entirely sure where to make the modifications.
I think this issue will be solved when the code reads and writes the file immediately.
Then I think I can do clean.cache() after each time step and avoid the issue.

Please let me know if this is possible.

Thank you very much.

Hi, do you know how can I get a input deck to visualize in Paraview? is it possible with your (awesome) ccx2paraview? I would be interested in observe CLOAD nodal forces as vectors just like the CGX -c mode, I think it could be useful to check correct setup in models.

Thank you for your work!

Could this be changed to use a generic python via PYTHONHOME?

Can't use:
pip3 install VTK
to install VTK on Win 10 Ver 1909. It says:

  ERROR: Could not find a version that satisfies the requirement VTK (from versions: none)
ERROR: No matching distribution found for VTK

```

I am trying to create a faster a binary with pyinstaller but the executable errors out with "No module found...".. Here is my full output below if there is something obvious I am doing wrong:

[feacluster@micro ccx2paraview]$ ll
total 52
-rwxrwxr-x. 1 feacluster feacluster  3796 Dec  5 19:54 __init__.py
-rwxrwxr-x. 1 feacluster feacluster  8113 Dec  5 19:54 frd2vtk.py
-rwxrwxr-x. 1 feacluster feacluster  1809 Dec  5 19:54 clean.py
-rwxrwxr-x. 1 feacluster feacluster  6305 Dec  5 19:54 writer.py
-rwxrwxr-x. 1 feacluster feacluster 20802 Dec  5 19:54 reader.py
drwxrwxr-x. 2 feacluster feacluster   153 Dec  5 19:54 __pycache__
-rw-rw-r--. 1 feacluster feacluster  1129 Dec  5 20:04 __init__.spec
[feacluster@micro ccx2paraview]$ pyinstaller __init__.py --onedir
43 INFO: PyInstaller: 4.7
43 INFO: Python: 3.6.8
43 INFO: Platform: Linux-4.18.0-348.2.1.el8_5.x86_64-x86_64-with-centos-8
44 INFO: wrote /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/__init__.spec
45 INFO: UPX is not available.
46 INFO: Extending PYTHONPATH with paths
['/home/feacluster/.local/lib/python3.6/site-packages']
176 INFO: checking Analysis
176 INFO: Building Analysis because Analysis-00.toc is non existent
176 INFO: Initializing module dependency graph...
177 INFO: Caching module graph hooks...
182 INFO: Analyzing base_library.zip ...
2227 INFO: Caching module dependency graph...
2314 INFO: running Analysis Analysis-00.toc
2336 INFO: Analyzing /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/__init__.py
3719 INFO: Processing pre-safe import module hook urllib3.packages.six.moves from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks/pre_safe_import_module/hook-urllib3.packages.six.moves.py'.
4707 INFO: Processing pre-safe import module hook six.moves from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks/pre_safe_import_module/hook-six.moves.py'.
5158 INFO: Processing module hooks...
5158 INFO: Loading module hook 'hook-lxml.py' from '/home/feacluster/.local/lib/python3.6/site-packages/_pyinstaller_hooks_contrib/hooks/stdhooks'...
5412 INFO: Loading module hook 'hook-lxml.etree.py' from '/home/feacluster/.local/lib/python3.6/site-packages/_pyinstaller_hooks_contrib/hooks/stdhooks'...
5412 INFO: Loading module hook 'hook-lxml.isoschematron.py' from '/home/feacluster/.local/lib/python3.6/site-packages/_pyinstaller_hooks_contrib/hooks/stdhooks'...
5415 INFO: Loading module hook 'hook-lxml.objectify.py' from '/home/feacluster/.local/lib/python3.6/site-packages/_pyinstaller_hooks_contrib/hooks/stdhooks'...
5415 INFO: Loading module hook 'hook-difflib.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5417 INFO: Loading module hook 'hook-encodings.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5456 INFO: Loading module hook 'hook-heapq.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5458 INFO: Loading module hook 'hook-multiprocessing.util.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5458 INFO: Import to be excluded not found: 'test'
5459 INFO: Loading module hook 'hook-pickle.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5460 INFO: Loading module hook 'hook-sqlite3.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5494 INFO: Loading module hook 'hook-xml.dom.domreg.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5494 INFO: Loading module hook 'hook-xml.etree.cElementTree.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5494 INFO: Loading module hook 'hook-xml.py' from '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks'...
5505 INFO: Looking for ctypes DLLs
5506 INFO: Analyzing run-time hooks ...
5509 INFO: Including run-time hook '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks/rthooks/pyi_rth_pkgutil.py'
5511 INFO: Including run-time hook '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks/rthooks/pyi_rth_inspect.py'
5512 INFO: Including run-time hook '/home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/hooks/rthooks/pyi_rth_multiprocessing.py'
5519 INFO: Looking for dynamic libraries
5823 INFO: Looking for eggs
5823 INFO: Using Python library /lib64/libpython3.6m.so.1.0
5828 INFO: Warnings written to /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/build/__init__/warn-__init__.txt
5856 INFO: Graph cross-reference written to /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/build/__init__/xref-__init__.html
5864 INFO: checking PYZ
5864 INFO: Building PYZ because PYZ-00.toc is non existent
5864 INFO: Building PYZ (ZlibArchive) /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/build/__init__/PYZ-00.pyz
6344 INFO: Building PYZ (ZlibArchive) /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/build/__init__/PYZ-00.pyz completed successfully.
6349 INFO: checking PKG
6350 INFO: Building PKG because PKG-00.toc is non existent
6350 INFO: Building PKG (CArchive) __init__.pkg
6364 INFO: Building PKG (CArchive) __init__.pkg completed successfully.
6365 INFO: Bootloader /home/feacluster/.local/lib/python3.6/site-packages/PyInstaller/bootloader/Linux-64bit-intel/run
6365 INFO: checking EXE
6365 INFO: Building EXE because EXE-00.toc is non existent
6365 INFO: Building EXE from EXE-00.toc
6365 INFO: Copying bootloader EXE to /home/feacluster/.local/lib/python3.6/site-packages/ccx2paraview/build/__init__/__init__
6365 INFO: Appending PKG archive to ELF section in EXE
6378 INFO: Building EXE from EXE-00.toc completed successfully.
6379 INFO: checking COLLECT
6379 INFO: Building COLLECT because COLLECT-00.toc is non existent
6379 INFO: Building COLLECT COLLECT-00.toc
6407 INFO: Building COLLECT COLLECT-00.toc completed successfully.
[feacluster@micro ccx2paraview]$ ls
build     dist        __init__.py    __pycache__  writer.py
clean.py  frd2vtk.py  __init__.spec  reader.py
[feacluster@micro ccx2paraview]$ cd dist
[feacluster@micro dist]$ ll
total 4
drwxrwxr-x. 5 feacluster feacluster 4096 Dec  5 20:06 __init__
[feacluster@micro dist]$ cd *
[feacluster@micro __init__]$ ls
base_library.zip  libexslt.so.0         libreadline.so.7  libz.so.1
__init__          libffi.so.6           libsqlite3.so.0   lxml
libbz2.so.1       libgcrypt.so.20       libssl.so.1.1     psutil
libcrypto.so.1.1  libgpg-error.so.0     libtinfo.so.6
lib-dynload       liblzma.so.5          libxml2.so.2
libexpat.so.1     libpython3.6m.so.1.0  libxslt.so.1
[feacluster@micro __init__]$ ./__init__
Traceback (most recent call last):
  File "ccx2paraview/__init__.py", line 21, in <module>
ModuleNotFoundError: No module named 'writer'

I executed following command with version 3.0.0,
ccx2paraview.exe exp_calculix.frd vtu

It gave me following error.

Fatal Python error: initfsencoding: unable to load the file system codec
ModuleNotFoundError: No module named 'encodings'

However everything was working fine in earlier version-2.3.1.

Please help.

Hi, the VTU tool also is able to animate modal analyses?

Thank you for your great work!

Problem description

When installing ccx2paraview under Ubuntu 22.04 using pip like

pip3 install --user ccx2paraview

pip will create a file ~/.local/lib/python3.10/site-packages/ccx2paraview/__init__.py that turns ccx2paraview into a Python module.

The __init__.py is empty, and therefore the module ccx2paraview -- when importing in a Python file -- has no members except again the submodule ccx2paraview (because the folder ~/.local/lib/python3.10/site-packages/ccx2paraview contains the file ccx2paraview.py).

Hence, in order to access the module elements from ccx2paraview one has to state from ccx2paraview import ccx2paraview and then all elements are again available like e.g. ccx2paraview.Converter

Desired behavior

After import ccx2paraview in a Python file have the relevant members of ccx2paraview available like ccx2paraview.Converter etc.

Potential fix

Include a file __init__.py into the ccx2paraview repository with content like

from .ccx2paraview import *

to import everything from the ccx2paraview.py or

from .ccx2paraview import Converter
from .ccx2paraview import <whatever a user needs>

I modified the empty __init__.py like this and it worked.

Good afternoon,

First I would like to thank you about the great job you are doing on this project.

I try to convert a calculix .frd file with the excutable ccx2paraview.exe. My issue is that I have an error (Numpy).

Can you give me an advice about to solve this issue. I tkink I am doing something wrong.

Thank you.

Raphaël

Here after the message :

C:\Users\RA\Documents\09_OUTILS_BE\calculix_\1\ccx2.16orig>ccx2paraview.exe fichier_parameter.frd vtu
INFO: Parsing fichier_parameter.frd
INFO: 192 nodes
INFO: 75 cells
INFO: Step 1, time 1.0, U, 3 components, 0 values
ERROR: Numpy is not installed.
Additional stresses will not be appended.
INFO: Step 1, time 1.0, S, 6 components, 0 values
ERROR: Numpy is not installed.
Additional strains will not be appended.
INFO: Step 1, time 1.0, E, 6 components, 0 values
INFO: 1 time increment
INFO: Writing fichier_parameter.vtu
INFO: Step 1, time 1.0, U, 3 components, 192 values
INFO: Step 1, time 1.0, S, 6 components, 192 values
INFO: Step 1, time 1.0, E, 6 components, 192 values

Hi,

I've noticed that when I use this converter and open the vtk file in ParaView, node and element numbers start from 0 instead of 1 like in the input file. I thought that it was a fault of ParaView but its dev told me that the converter itself does it. It can be a problem when exporting results of ParaView operations at nodes and cells (e.g. to map them to a mesh). Would it be possible to change this behavior ?

I would like to know if you are going to provide a binary release (bin) of the ccx2paraview application, as you did with 2.3.1. In any case, thank you very much for your contribution.

It seems ccx used brick elements instead of direct 2D elements for the analysis. When output the results, is that possible to convert the nodes and elements number back to the original 2D nodes and elments numbers?

Hi,
I've been using this converter for a few days now and I really like it. However, I noticed that stress components seem to be stored as vectors. This means that tensor operations in the python calculator such as eigenvector(S) are not possible. Would it possibly make sense to store tensor quantities as tensors? Or is there a specific reason that it is exported as a vector?

Currently I convert the vectors to tensors in a ProgrammableFilter:

def make_tensor(xx,yy,zz,xy,yz,xz):
t = np.array([])
for i in range(len(xx)):
t_i = np.array([[xx[i],xy[i],xz[i]],[xy[i],yy[i],yz[i]],[xz[i],yz[i],zz[i]]])
t = np.append(t, t_i)
t = dsa.VTKArray(t)
t.shape=(len(xx),3,3)
return t

sxx = inputs[0].PointData["S"][:,0]
syy = inputs[0].PointData["S"][:,1]
szz = inputs[0].PointData["S"][:,2]
sxy = inputs[0].PointData["S"][:,3]
syz = inputs[0].PointData["S"][:,4]
sxz = inputs[0].PointData["S"][:,5]

stress = make_tensor(sxx,syy,szz,sxy,syz,sxz)

But this is very slow.

Hi,
While using this tool, I have often found myself longing for the possibility to only write a single timestep.
Sometimes I simply want to check on a running simulations last time step, or have previously written out .vtus for previous timesteps and now want to write the remaining steps. Currently, it is only possible to write out every time step again. I think it would be a cool feature to have simple arguments like:
-latestTime
-newTimes
-time (specify a range of steps to write)
This is something I often use with OpenFOAM tools like reconstructPar (https://www.openfoam.com/documentation/guides/latest/man/reconstructPar.html)
I was wondering if you think this would be a good feature to implement?

Hello,
First of all, thank you for this tool!
I ran into an issue using the conversion from .frd to vtu. The conversion to to .vtu works just fine resulting in files like:
shell.01.vtu
shell.02.vtu
....
The created .pvd file then contains

		<DataSet file="shell.01vtu" timestep="0.002"/>
		<DataSet file="shell.02vtu" timestep="0.004"/>

Due to the missing dot, the correct .vtu files can't be opened.
I could resolve this issue by changing the following section:

f.write('\t\t<DataSet file="{}.vtu" timestep="{}"/>\n'\ 
         .format(os.path.basename(file_name), inc))

Can you confirm this behaviour on your side?

Hi,
when I plot a mesh using Hexa quads C3D20's, I can see that the element is built with some wrong order in the node list, making them to look somehow twisted. This is the same both in VTk and VTU scripts.

Edite: Corrected type from C3D30s to C3D20's

Hi,

I'm using "view vec" in cgx to produce so-alike stress trajectory. How can i simulate it using paraview .vtu using glyph filter? in paraview i cannot see it in orientation array (i can only see displacement). i notice that principal stress are post-produce. ty, sorry for my bad english.

and i think i cannot get ccx "worstPS" in vtu produced file?

i dont kown what's going wrong!

Your translator looks very cool! I want to post-process small models in the web browser. I think ParaViewWeb can do that. Have you ever tried setting that up? If not, do you think it will work with your vtk files?

Presently I do some basic screenshots from cgx via a batch file. See:

https://www.feacluster.com/web_calculix.php

I want to enhance that to allow people to do full 3d visualization over the web.

I suppose

p = FRDParser.Parse(self.file_name)

contains mesh nodes / elements as well as fields defined on nodes / elements. Instead of writing the legacy PVD / VTK files (one file per time step), a more modern format to use is XDMF which supports natively time-steps. Maybe meshio can be used to a generate an easy-to-use XDMF output.

with meshio.xdmf.TimeSeriesWriter(filename) as writer:
    writer.write_points_cells(points, cells)
    for t in [0.0, 0.1, 0.21]:
        writer.write_data(t, point_data={"phi": data})

https://calculix.discourse.group/t/converter-ccx2paraview-for-large-frd-files/1259/15

Hello,

Just to report a bug on the ccx2paraview converter. Not sure if is the right place. But anyway here it is.
I think the author keeps an eye into the forum, so first of all thanks for coding this very useful tool (at least for the workflow i use).

The bug occurs when placing a perturbation step for a frequency analysis after a static step.

Example:

** == ANALYSIS SECTION ==
**
*TIME POINTS, NAME=T0
0.200,0.400,0.600,0.800,1.000
*STEP
*STATIC
*BOUNDARY
RefNode,6,6,0.03
*EL FILE, TIMEPOINTS=T0
S,E
*NODE FILE, TIMEPOINTS=T0
U,RF
*END STEP
*STEP,PERTURBATION
*FREQUENCY
10
*EL FILE
S,E
*NODE FILE
U,RF
*END STEP

Error:

[...]
INFO: Writing n-twist_model.05.vtu
Traceback (most recent call last):
File "/home2/jl/workspace/software_devel/ccx2paraview/ccx2paraview.py", line 992, in
main()
File "/home2/jl/workspace/software_devel/ccx2paraview/ccx2paraview.py", line 987, in main
ccx2paraview.run()
File "/home2/jl/workspace/software_devel/ccx2paraview/ccx2paraview.py", line 904, in run
result_blocks = self.frd.parse_results(step, inc) # NOTE Could be empty list []
File "/home2/jl/workspace/software_devel/ccx2paraview/ccx2paraview.py", line 696, in parse_results
result_blocks.append(self.calculate_principal(b))
File "/home2/jl/workspace/software_devel/ccx2paraview/ccx2paraview.py", line 780, in calculate_principal
eigenvalues = np.linalg.eigvals(tensor).tolist()
File "<array_function internals>", line 5, in eigvals
File "/usr/lib/python3/dist-packages/numpy/linalg/linalg.py", line 1063, in eigvals
_assert_finite(a)
File "/usr/lib/python3/dist-packages/numpy/linalg/linalg.py", line 209, in _assert_finite
raise LinAlgError("Array must not contain infs or NaNs")
numpy.linalg.LinAlgError: Array must not contain infs or NaNs

So the error occurs when trying to convert the frequency steps.

The procedure finishes without error if we comment the code between lines 694 to 700.

               if b.name == 'S':
                    result_blocks.append(self.calculate_mises_stress(b))
                    result_blocks.append(self.calculate_principal(b))
                if b.name == 'E':
                    result_blocks.append(self.calculate_mises_strain(b))
                    result_blocks.append(self.calculate_principal(b))

but there is still a warning:

[...]
INFO: Writing n-twist_model.05.vtu
INFO: Step 6, time 2750.3, U, 3 components, 34432 values
INFO: Step 6, time 2750.3, S, 6 components, 34432 values
WARNING: 69264 NaN values are converted to 0.0
INFO: Step 6, time 2750.3, E, 6 components, 34432 values
INFO: Step 6, time 2750.3, RF, 3 components, 34432 values
WARNING: 48780 NaN values are converted to 0.0
INFO: Step 6, time 2750.3, ERROR, 1 components, 34432 values
[...]

to be honest to my workflow it would be usefull if there was a way to select if to compute these additional quantities or not (misses and principal values). This is because I use the vtk api to further compute additional quantities (including the principal vectors, which make the initial computation of principal values redundant).

Hope this info is useful,
thanks again,
regards,
Jorge

Hello,
So I wanted to run ccx2paraview inside a venv. I ran

python3 -m venv venv
source venv/bin/activate
pip install numpy vtk ccx2paraview

and then I tried to run the script

#!/usr/bin/env python3
# Converts calculix results to vtk or vtu
import ccx2paraview
c = ccx2paraview.Converter("beam.frd", 'vtu')`
c.run()
```

But this only creates an empty file named ccx2paraview

The command

pyinstaller src/ccx2paraview.py --onefile

runs successfull. But in a fresh console I get

ccx2paraview ....frd vtk
Traceback (most recent call last):
  File "site-packages/PyInstaller/loader/rthooks/pyi_rth_pkgres.py", line 13, in <module>
  File "/usr/lib/python3.8/site-packages/PyInstaller/loader/pyimod03_importers.py", line 623, in exec_module
  File "site-packages/pkg_resources/__init__.py", line 84, in <module>
ModuleNotFoundError: No module named 'packaging.requirements'
[12] Failed to execute script pyi_rth_pkgres

VTK introduced the VTKHDF file format, which supports temporal data.

This format could be a 1-file output solution. case.frd -> case.vtkhdf

Check variables CSTR and CDIS.

Hi there.
Thanks for this utility. It's great to find something which simply works without hassle. (ubuntu 16.04)
Here's a feature suggestion: Its great to have the derived stress quantities, but if we want to do things like 3D glyph plots of stress we would need the stress tensor. Of course we can reconstruct it with the programmable filter but it would be great not to need that.
One small question: does ccx2paraview assume that the node numbering and/or element numbering are contiguous? I work from CAD files and end up zapping a lot of element sets. This leaves gaps and I suspect that ccx2paraview doesnt allow that.

Best,Peter

Hi,

thank you for sharing the ccx2paraview code! I was just curious about your thoughts on creating a pip (or conda if that is preferable) package of ccx2paraview?

Having a distributable python package of ccx2paraview seems like something that could make it a bit easier to keep track of package versions when building docker images containing ccx (which is something I am working on)! if you want I can create a pull request with a pypi pipeline setup. Or if you prefer to keep the repo focused on its current tasks, I completely understand and I will instead continue with my fork of ccx2paraview.

Let me know what you think,

Kind Regards
Kristoffer

Hello,

I am really glad I found this because I am new to CCX and CGX, but I have experience with ParaView. Also, I have some issues with visualizing results using cgx -v file_name.frd.

I was able to use it on one of my results.
I opened the final time step(500th-time step) and confirm that the results I am seeing in ParaView.
I read your warning about multiple files being created for each time step.

Is it possible to open all the files at once?

I don't know if I should open a new issue for this or not since I am a new user.
Is it possible to see the deformed shape in ParaView?

In any matter, I am really glad I found this since this can at least work as a temporary fix for my problems until I figure out CGX.

Thank you @imirzov for ccx2paraview.

PS: Here are the results

Here,you can see the space which occupies the flap being deformed, but would be nice to see the first image deformed.