Hi Hasan, I loaded up your mesh in gmsh and it looks like something went wrong with the extrusion: I only see quadrilaterals on two opposing faces and no volume mesh at all. I think that if you fix the mesh you will get past this error.
I attached a picture of what I see in gmsh to this message. Thanks, David On Thu, Jan 2, 2020 at 2:03 PM Mohammed Hassan <[email protected]> wrote: > > Hi Dealii Users; > > I have a similar problem, I am running a 3D problem for a beam with the > following dimensions x=10 mm and y= 1 mm and z =1 mm I am trying to adapt > step-8 to my case the I am facing the following challenges : > > 1- importing msh file to the dealii, I copied a print_mesh_info function to > see the vtu file but the following error pop up. > > An error occurred in line <2018> of file > </home/mohammed/FEM/dealii/dealii/source/grid/grid_in.cc> in function > void dealii::GridIn<dim, spacedim>::read_msh(std::istream&) [with int dim > = 3; int spacedim = 3; std::istream = std::basic_istream<char>] > The violated condition was: > cells.size() > 0 > > 2- the purpose of this work is to learn how to implement the boundary > conditions. the .msh file has been meshed as hexahedral where it has 100 > cells along x axis and 10 cells along y and 10 cells along z axis. I would > like to apply 100 N force on the right side of the beam pointing to ( + x > axis and the second case + y axis " I mean the direction of the force") and a > fixed wall on the left side knowing that I named the cells for the right side > of the beam (in .msh file) force while for the left side of the beam wall. > > Thank you in advance ! > > Kind Regards. > > > > > > On Tuesday, December 10, 2019 at 4:35:20 PM UTC+1, Hasan Tasneem wrote: >> >> I am a beginner on Deal.ii. I have installed and run some of the tutorials. >> I am now trying to run one specific problem. This problem is 2D and is a >> beam problem. My dimension of the beam is 3*100mm. The beam is fixed on one >> end and there is 1000N of transverse loading on the other end. I am >> following Step-8 which solves an elastic equation. Everything I have >> inputted correctly. However, I am stuck on applying the fixed support. I >> want the fixed support at x equal to 0. I don't know how to proceed further. >> >> VectorTools::interpolate_boundary_values(dof_handler, >> 0, >> Functions::ZeroFunction<dim>(dim), >> boundary_values); >> MatrixTools::apply_boundary_values(boundary_values, >> system_matrix, >> solution, >> system_rhs); >> >> How should I alter this to get fixed support at x equal to 0. >> >> > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/db40fc39-9631-4eed-90e8-0def8efb1441%40googlegroups.com. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CABrTbYSXTN4Ewi5VcryBAU%3DR5r9s2TAqZyURpJXghsDvJnkzsg%40mail.gmail.com.
