When I run your program, I have a problem with the eval methods :
mfpt.eval("c*(x[1]-r) * I",{'c':c, 'r':r_tube, 'I':I})
p_o = mfpo.eval("p(x[2]) * I",{'p':p_atm , 'I':I})
p_i = mfpi.eval("p(x[2]) * I",{'p':p_tube, 'I':I})
Normally, x[1] and x[2] are not defined. The available coordinates
are the scalar variables x, y and z.
But I have an error message ...
Yves.
Le 20/02/2012 21:38, address@hidden a écrit :
Dear all,
my getfem++ model is killed with signal 9 and I've no idea
why.
I try to use gmsh and getfem++ to calculate the deformation of
a pressurized tube. I model the tube as a thin walled cylinder
with top and bottom in gmsh. I use python as the interface to
getfem++ to describe the high inside pressure, the lower
outside pressure, the asymmetric top pressure and a diriclet
bottom. When I run this programm it starts, but it crashes
without any hint why.
I tried this with different length of the cylinder, with
different thickness of the walls, and with meshes of different
granularity. The result is always the same: the program is
killed. Only with a very small flat cylinder it succeeded, but
I'm interested in the behaviour of a long tube-like cylinder.
I've suspected the main memory, but I've added plenty of
virtual memory without an observable change.
Can you give me any hints how to overcome the problem or what
I'm doing wrong?
Many thanks in advance.
Kind regards,
Radames
---
Nutzen Sie
freenet Mail optimal angepasst für Ihr iPhone, Android oder
Nokia Handy auch von unterwegs.
Alle Infos und Download unter freenetMail Apps
_______________________________________________
Getfem-users mailing list
address@hidden
https://mail.gna.org/listinfo/getfem-users
--
Yves Renard (address@hidden) tel : (33) 04.72.43.87.08
Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
---------
|