Tochnog Finite Element Analysis Files

 UPDATES TO TOCHNOG 6-26-2012

Dear Tochnogers,

Following the recent postings regarding tochnog and contact issues, I remembered having similar issues before and reviewed the options and problems. 
I realized that contact between a mesh and a point works, as it is shown in example 9. Contact between a mesh and a circle worked also (Please see attached problem)
contact with a line works perfectly well (see example 14) but contact between an ellipse and a mesh had issues and did not work.

After long hours looking at the problem, I have edited and corrected the errors in the code in file "geometry.cc" and it works very well, as it was supposed to.

I am attaching a copy of the corrected files, in the June 26 2012 folder so that those building their own versions can update their code. I am also including three small examples that show contacts between geometries (point, circle and ellipse) and meshes. I realized a small error in two of the calculations and I am including the corrected versions.

The updated files contain the  cleanup modifications to area.cc, check.cc, database.cc, force.cc, math.cc , miscel.cc, pri.cc, time.cc, tochnog.h kindly contributed by Roman Putanowicz who kindly cleaned up the files to get rid of warnings during building the binaries. THANK ROMAN!!! 

The updated files to build the binaries are in the Tochnog-Lapack-June-26-2012.zip. The Mac OSX binary is named tochno-OSX-June-26-2012.zip.

I would ask those of you who have access to the 64 bit versions of linux or windows and would like to contribute with binaries, please send them to me.
I will post them so that other users can download them. I would ask that you test your binaries before sending them to me.

If you can, for those using Windows, I know that there is a LAPACK version already built for windows. It would be nice if it can be built with LAPACK for
Windows, so that all versions would have the same capabilities.

If you do not have the LAPACK library option in your build please edit the corresponding line in the attached example files. The three small attached examples 
solve the same problem for contact with a point 0.3 units in diameters, a circle of 0.3 units radius, and an ellipse with a=0.5  and b=0.3 all having the same center.

For those trying the contact between meshes, remember to include the "node_boundary" line in the problem.

Let me have some feedback please!

Best regards,

Fernando Lorenzo