Re: [Libmesh-users] Problems in simple program

2014-08-23 Thread Rodrigo Broggi
Ok guys, I think I've solved it: it was as predicted the matrix assembling. If anyone want the code let me know... Thanks, Rodrigo Broggi 2014-08-20 23:45 GMT+02:00 John Peterson : > > > > On Wed, Aug 20, 2014 at 3:36 PM, Rodrigo Broggi > wrote: > >> Yeah that's write but the pressure in both

Re: [Libmesh-users] Problems in simple program

2014-08-21 Thread Rodrigo Broggi
Hey guys, I really can't figure it out what's going on with that program. I've printed the matrix but actually this doesn't helped me at all. I still think the problem is most likely to be in the assembling cause I've not completely understand how the shape functions are handled and I've used the

Re: [Libmesh-users] Problems in simple program

2014-08-20 Thread John Peterson
On Wed, Aug 20, 2014 at 2:39 PM, Rodrigo Broggi wrote: > ​ > Articolo.pdf > > ​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::cou

Re: [Libmesh-users] Problems in simple program

2014-08-20 Thread Rodrigo Broggi
​ Articolo.pdf ​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 od

Re: [Libmesh-users] Problems in simple program

2014-08-20 Thread John Peterson
On Wed, Aug 20, 2014 at 8:42 AM, Rodrigo Broggi 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

Re: [Libmesh-users] Problems in simple program

2014-08-20 Thread Roy Stogner
On Wed, 20 Aug 2014, Rodrigo Broggi wrote: > 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? Something like: // libMesh::out.precision(16); // if you need libMesh::out << "J = [" << *(syst

[Libmesh-users] Problems in simple program

2014-08-20 Thread Rodrigo Broggi
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 th