Dear all,
I asked this question at the old mail group and I got a few replies. The most of them recommended, I should use instead of spooles the pardiso solver, but sadly this is not solving my problem.
First of all, here again is my problem description:
I want to calculate bone healing. For this it is necessary to assign to each tetrahedron an own material, thus each element has to be a element set.
To specify my problem I carried out a simple test:
-I created a mesh representing a bone
-Just for this test, every elements has the material steel
-At the top of the bone a force is defined, at the bottom the bone is fixed
-The mesh and the load case are for both cases the same
Case 1:
-I created for each element an element set, and assign the material steel to each element set
-The number of materials are the same as the number of elements
Case 2:
-I created 6 element sets which contain the same elements as in case 1, but these sets “bundle up” these elements
-There is just one material definition
The resulting equation system is in both cases the same, just the amount of element sets and materials are different, as you can see in the output files here (the input files can be found there too)
The calculation time of spooles is in both cases the same, but overall case one takes 2183 seconds and case two 774 seconds. Calculix needs at case one about 25 minutes longer to read the input file. For me it looks like calculix is not made for a lot of different element sets and materials, what means my work is over here, and I have to move to another solver, because I have to refine the mesh what means I have much more element sets. (atm I use a coarse mesh for testing)
Has someone an idea how to fix this or why calculix has that much problems with a large amount of element sets/material definitions?
Best regard,
Patrick