Articolo.pdf <https://docs.google.com/file/d/0B49xyYU5aZowazV5TmNyeS1GYXc/edit?usp=drive_web> Thanks for the attention guys. Here it goes the file. About the program, I don't get any errors and actually there is also a "std::cout" that garantees that the PressurexR works fine. It's really an odd result, I can even postprocess on Paraview but the solution is just constant and zero. I've also changed the elastic constants to see if there was a problem with magnitudes but still got the same problem. Thanks,
2014-08-20 22:24 GMT+02:00 John Peterson <[email protected]>: > On Wed, Aug 20, 2014 at 8:42 AM, Rodrigo Broggi <[email protected]> > 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 [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
