Hi,
I recently converted an ABAQUS input file into UCD mesh using convert mesh
as follows,
nml-staff@ubuntu:~/deal.II/contrib/mesh_conversion$ ./convert_mesh 3 0
~/Desktop/input.inp ~/Desktop/output.ucd
Dimension: 3
************************************************************************************
*** FEM Mesh conversion tool
***
************************************************************************************
*** Author: Jean-Paul Pelteret
***
*** Date: June 2010
***
*** References:
***
***
http://www.dealii.org/developer/doxygen/deal.II/structGeometryInfo.html
***
*** http://people.scs.fsu.edu/~burkardt/html/ucd_format.html
***
*** http://people.scs.fsu.edu/~burkardt/data/ucd/ucd.html
***
*** http://www.cprogramming.com/tutorial/string.html
***
*** http://www.codeguru.com/forum/showthread.php?t=231054
***
************************************************************************************
*** FEATURES:
***
*** Read-in file types:
***
*** - Abaqus inp (from Cubit)
***
*** Write-out file types:
***
*** - AVS UCD
***
************************************************************************************
Reading in ABAQUS .inp FILE: /home/nml-staff/Desktop/input.inp
Writing out AVS .ucd FILE: /home/nml-staff/Desktop/output.ucd
Output successful!
Number of nodes: 10610
Number of cells: 9408
Number of boundary faces: 0
and everything seems to be working well. However when I read the UCD file in
and try to output it in any form (eps, gnuplot,...)
I get the following error,
*nml-staff@ubuntu:~/Desktop/fem/step-1$ make run*
*============================ Running step-1*
*--------------------------------------------------------*
*An error occurred in line <11617> of file
</home/nml-staff/Desktop/deal.II/source/grid/tria.cc> in function*
* typename dealii::Triangulation<dim, spacedim>::raw_hex_iterator
dealii::Triangulation<dim, spacedim>::begin_raw_hex(unsigned int) const
[with int dim = 3, int spacedim = 3]*
*The violated condition was: *
* level<levels.size()*
*The name and call sequence of the exception was:*
* ExcInvalidLevel(level)*
*Additional Information: *
*The given level 0 is not in the valid range!*
*
*
*Stacktrace:*
*-----------*
*#0 /home/nml-staff/Desktop/deal.II/lib/libdeal_II.g.so.7.0.0:
dealii::Triangulation<3, 3>::begin_raw_hex(unsigned int) const*
*#1 /home/nml-staff/Desktop/deal.II/lib/libdeal_II.g.so.7.0.0:
dealii::Triangulation<3, 3>::begin_hex(unsigned int) const*
*#2 /home/nml-staff/Desktop/deal.II/lib/libdeal_II.g.so.7.0.0:
dealii::Triangulation<3, 3>::begin_active_hex(unsigned int) const*
*#3 /home/nml-staff/Desktop/deal.II/lib/libdeal_II.g.so.7.0.0:
dealii::Triangulation<3, 3>::begin_active(unsigned int) const*
*#4 /home/nml-staff/Desktop/deal.II/lib/libdeal_II.g.so.7.0.0: *
*#5 /home/nml-staff/Desktop/deal.II/lib/libdeal_II.g.so.7.0.0: void
dealii::GridOut::write_gnuplot<3, 3>(dealii::Triangulation<3, 3> const&,
std::ostream&, dealii::Mapping<3, 3> const*) const*
*#6 ./step-1: first_grid()*
*#7 ./step-1: main*
*--------------------------------------------------------*
*make: *** [run] Aborted*
Does anyone knows what is going wrong? What is "levels.size"?
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
