Dhondtguido / CalculiX

This repository contains the source files of CalculiX, a three-dimensional Finite Element Program (www.calculix.de).
GNU General Public License v2.0
72 stars 16 forks source link

Changes in *BUCKLE definition #57

Open lucas-bueno777 opened 10 months ago

lucas-bueno777 commented 10 months ago

I want to ask for two main things: the first is about the default accuracy in *Buckle cards. As the documentation states, the default accuracy is 0.01, and most of the time it should not be changed.

image

However, I did some tests that are presented in the CalculiX forum and when the default settings are used, Pardiso, PastiX and Spooles fail to find the first eigenvalue when only one is requested. However, if we change the accuracy to 0.0001 (100 times more accurate than default), all solvers now can handle the search for the first eigenvalue properly with absolutely minimal changes in time spent in the run for my test case. Therefore, my first point is to improve the default accuracy of the *Buckle card from 0.01 to 0.0001 at least.

image

The second thing I would like to ask is to add control parameters to change the buckling definition the same way they are presented in the *Frequency card. This way, one can easily avoid negative eigenvalues and search only within a defined range if is desired.

image

Best regards, Lucas

victorkemp commented 10 months ago

I agree about the accuracy parameter. This is the "tol" parameter passed to Arpack. Arpack's default value for that is machine precision which I guess is about 1e-15 or so. That's what I use in my other solver and it's fine. But anything better than 0.01 would help.