On Wed, Aug 20, 2014 at 8:42 AM, Rodrigo Broggi <rbrogg...@gmail.com> wrote: > Hello guys, > > I wrote this program for an elasticity axisimmetric problem of concentric > pipes. I took the example 4 of system equations as basis and made some > changes. The problem is that the result is completely wrong: even thought I > had imposed internal and external pressures and self weight in the right > hand side my output in displacements is constant all over the pipe. At this > point I can only think that I've messed up with the matrix assembling. > Someone knows how can I print the stiffness matrix once assembled?
I don't see any of the "usual" libmesh coding gotchas in your code. Are you sure that PressurexR is actually working? Floating point comparisons like this > if (x == R_i) > return R_i*(104e6 - 10000*y); > > else if (x == R_e) > return R_e*(1.5e7 - 10500*y); are generally frowned upon, but I see you have that error message to handle that case... > I've attached a document with the weak formulation of the problem (2D) I'm > trying to solve - it's the equation 9 without the thermal contribute. Nothing was attached for me, our mailing list typically deletes attachments so please share a google doc instead. -- John ------------------------------------------------------------------------------ Slashdot TV. Video for Nerds. Stuff that matters. http://tv.slashdot.org/ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
