/* UNIVERSAL - CalculiX interface program, Bernhardi, Aug 2011 */
/****************************************************************************
* THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
* IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
* AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
*
* THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
* WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
* MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
* LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
*
* YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
* ALONG WITH THIS PROGRAM; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION, INC,
* 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
*****************************************************************************/
/*
This small program converts finite elemet meshes given in the SDRL universal
file format into ABAQUS / Calculix compatible input decks. The format
of the universal data sets can be found at http://www.sdrl.uc.edu.

The element - node mapping is taken from the Python unv2abaqus script, which
can be found at /www.caelinux.org.

This program will still be far from being error - free; therefore, use
it with caution!

The program also translates surface meshes in case they are present in
the .med file: If the main mesh is planar, the line elements (if they
exist) are translated to abaqus elsets with load surface IDs P1...P4.

Correspondingly, if a volume mesh is given, the surface meshes (if they
are also present in the med file) are converted to elsets with load
IDs P1...P6.

The *elset names will be the same as the group names defined in SALOME.

I use that feature to introduce boundary conditions. However, the
calculation of the surface elsets takes quite long for very large meshes.

The program works for all CalculiX element types. It compiles under
opensuse 12.3 (but should compile on all linuxes).

To install and use:

- place the program source text into a proper directory, e.g. /path/unical

- go to that directory.

- issue a "gcc -o2 -o unical unical1.c" command.

- try it out using, for example, a universal file cube.unv:
"/path/unical cube".

- Then, configure your Abaqus/Calculix analysis run by modifying the
cube.inp file: Introduction of loads, boundary conditions, prodedures,
etc.

Remarks:

- Unical translates the following three Universal data sets:
2411, node point numbers and coordinates
2412, element type codes and connectivities
2467, groups of nodes and elements.
Everything else is ignored.

- Only tested with universal files generated by SALOME.

- Planar elements are always translated to shell elements, e.g. S8R. Thus
the *element, type= lines must be modified if e.g. CPS8 elements are to be
used.

- In case of a volume mesh, SALOME also generates surface elements;
correspondingly, in case of a planar mesh the boundary is also output
using line elements. In most cases the boundaries are not needed and
are therefore not translated to CalculiX; rather, the surface elements
are used to generate *surface input lines for ease of boundary condition
definition, as described above.

- However, there are cases where different element dimensions are desired;
for example, combination of volume elements and beam elements. If this is the
case, specify "./unical -f cube". Then, all elements are translated. If only
volumes and beams are wanted, shell elements must be manually removed from
the .inp file.