Hello Professor Bangerth, Thanks for the reply. After doing multiple review on the code, I found I had some errors in the assembly. some brackets were misaligned and off so the global contributions were added at each quadrature point!
Now the linear, quadratic and cubic shape functions all approach the same solution as Euler-Bernoulli ode. Many thanks, Mohammed On Friday, March 29, 2024 at 2:55:07 AM UTC+9 Wolfgang Bangerth wrote: > > > I noticed that when I use linear shape functions with an increasing > mesh x > > resolution, the solution seems to approach euler-bernoulli ode solution. > > however, when I use higher order shape functions, fem solution do not > converge > > as the linear shape function case. > > This would be concerning. The finite element approximation should converge > to > the same solution whether you use linear or quadratic elements. It should > just > converge faster with quadratic elements, but the limit point of the > convergence needs to be the same. > > If that's not the case, it's worth exploring why that would be so, perhaps > starting with articulating how exactly the two limits differ. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/fe511216-11cb-4380-bce5-67c75e9e3d51n%40googlegroups.com.
