Looking for a help with 2D stat mech model

Hello,
Can someone, please, help me to solve a simple static mechanical 2D model.
Every time I try to solve it, the same error message appears:


Matrice non factorisable:
pivot vraiment nul a la ligne: 110
pour le noeud N46 et la composante DY

Cette erreur est fetale. Le code s'arrete.


Here is a geo file of 2D model:

---

//Mesh density definitions
L=0.25;

//Geometric entities
Point(1) = {0,0,0,L};
Point(2) = {2,0,0,L};
Point(3) = {2,2,0,L};
Point(4) = {1.5,2,0,L};
Point(5) = {0,2,0,L};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,5};
Line(5) = {5,1};

Line Loop(11) = {1,2,3,4,5};
Plane Surface(1) = {11};

Physical Point (1) ={1};
Physical Point (2) = {4};
Physical Line (3) = {2};
Physical Surface (4) = {1};

---

And here is a comm file for that 2D model above:

---

DEBUT();

PRE_GMSH();

mail=LIRE_MAILLAGE();

FEMLin=AFFE_MODELE(MAILLAGE=mail,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='C_PLAN',),);

ortho=DEFI_MATERIAU(ELAS_ORTH=_F(E_L=3e5,
E_T=2e5,
E_N=1e5,
NU_LT=0,
G_LT=0,),
INFO=2,);

Mat=AFFE_MATERIAU(MAILLAGE=mail,
MODELE=FEMLin,
AFFE=_F(TOUT='OUI',
MATER=ortho,),);

thickn=AFFE_CARA_ELEM(MODELE=FEMLin,
INFO=1,
MASSIF=_F(GROUP_MA='GM4',
ANGL_REP=0,),);

BCnd=AFFE_CHAR_MECA(MODELE=FEMLin,
DDL_IMPO=(_F(GROUP_NO='GM1',
DY=0.0,),
_F(GROUP_MA='GM3',
DX=0.0,),),
FORCE_NODALE=_F(GROUP_NO='GM2',
FY=-1000,),
INFO=1,);

Solution=MECA_STATIQUE(MODELE=FEMLin,
CHAM_MATER=Mat,
CARA_ELEM=thickn,
EXCIT=_F(CHARGE=BCnd,),);

Solution=CALC_ELEM(reuse=Solution,
RESULTAT=Solution,
OPTION=('SIEF_ELNO_ELGA','SIGM_ELNO_DEPL',),);

IMPR_RESU(MODELE=FEMLin,
RESU=_F(RESULTAT=Solution,),);

DEFI_FICHIER(ACTION='ASSOCIER',
UNITE=30,);

IMPR_RESU(MODELE=FEMLin,
FORMAT='GMSH',
UNITE=30,
RESU=_F(RESULTAT=Solution,),);

DEFI_FICHIER(ACTION='LIBERER',
UNITE=30,);

FIN();

I'd be rally very thankful if someone helped me to solve this model.

Thanks in advance,
Vito

Post edited by: vito, at: 2007/09/19 18:25

Post edited by: vito, at: 2007/09/19 18:28

Post edited by: vito, at: 2007/09/19 18:32<br><br>Post edited by: vito, at: 2007/09/19 18:43