On 12/10/19 8:35 AM, 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 > <https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FnamespaceVectorTools.html%23a187aeb575be07bc47cb3dea1a47aaf88&data=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C8c7a72a1adeb4a51a33808d77d8692da%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C1%7C637115889269529807&sdata=dBtpLgO2b%2F7VxS%2BnTkfbUklV0TZJGp9o%2FKFAoA194Rk%3D&reserved=0>(dof_handler, > 0, > Functions::ZeroFunction<dim> > <https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FclassFunctions_1_1ZeroFunction.html&data=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C8c7a72a1adeb4a51a33808d77d8692da%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C1%7C637115889269539798&sdata=Ig%2B6Dk%2FJQk2UJfaTMp210MzgdAKDl6akaT6EsGLTlQI%3D&reserved=0>(dim), > boundary_values); > MatrixTools::apply_boundary_values > <https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FnamespaceMatrixTools.html%23a9ad0eb7a8662628534586716748d62fb&data=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C8c7a72a1adeb4a51a33808d77d8692da%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C1%7C637115889269539798&sdata=VruB70ggs57EfJA5%2BuULNSVhUyDTvI%2BqwIZZwD1TtG4%3D&reserved=0>(boundary_values, > system_matrix, > solution, > system_rhs); > > How should I alter this to get fixed support at x equal to 0.
In your case, you want to apply the zero Dirichlet boundary conditions only on the left side of the beam. To do this, you need to "mark" the different parts of the boundary with different boundary indicators -- take a look a some of the tutorial programs that set the "boundary_id". Then you need to call interpolate_b_v with only that boundary_id that applies to the left side of your domain. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- 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/50e07b8b-ffa3-221c-d031-1b2788a34200%40colostate.edu.
