Also I checked with a 4 element mesh and 1 quadrature point. However, I see a noninvertible mass matrix. I do not understand what is missing in the code.
Thank you. On Saturday, December 14, 2024 at 6:41:36 PM UTC-6 Pratyusa Kar wrote: > Hi, > > I have a simple 2D linear elasticity problem, a square domain x(0,1), > y(0,1). In that I am fixing x=0 in all directions and applying a > displacement of 0.02 along x direction at the edge x=1. > > I am facing issues with the strain L2 function convergence in my code. > Couple of checks I already did: > - The main solver and field calculation runs fine. Gives meaningful plot > of u which shows boundary conditions are fine. > - In the compute strain projection function, I printed > strain_xx_cell_rhs(i) (which prints ux,0 ) > - I also printed cell_matrix. > > The code does not converge with the CG solver. Since the problem is > simple, I am not sure if I am missing something. Any help will be > appreciated. > > I have attached the main.cc and a snapshot of the error. > > Thank you. > > > -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion visit https://groups.google.com/d/msgid/dealii/3a7d5a27-6a87-4d40-8657-86f59a66e8a1n%40googlegroups.com.
