Shell element problem

Cyl-noz.comm.tar.gz
I am trying to model part of a pressure vessel with shell elements, but I get stuck very early in my comm file.

In AFFE_CARA_ELEM it accepts the first mesh group, but for the second one the error message is that is not part of the mesh. This message is irrespective of which mesh group is first or second, so my logical conclusion is that I am missing something in my comm file.

I attach a copy of the comm file. (The comment lines may not be valid for this application.) Any guidance is highly appreciated.

Thank you in advance.

I can't get the file to upload, so I am just copying the relevant section of the file below:

DEBUT(LANG='ENG',);

MA=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,),);

Inmesh=LIRE_MAILLAGE(FORMAT='MED',);

Shmesh=CREA_MAILLAGE(MAILLAGE=Inmesh,
MODI_MAILLE=_F(TOUT='OUI',
OPTION='TRIA6_7',),);

Model=AFFE_MODELE(MAILLAGE=Shmesh,
AFFE=(_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='COQUE_3D',),
_F(GROUP_NO='C_node',
PHENOMENE='MECANIQUE',
MODELISATION='DIS_TR',),),);

Carel=AFFE_CARA_ELEM(MODELE=Model,
COQUE=(_F(GROUP_MA='Bshell',
EPAIS=8.0,),
_F(MAILLE=('Mshell','Nozzle',),
EPAIS=10.0,),
_F(GROUP_MA='Pad',
EPAIS=18.0,),
_F(GROUP_MA='Cover',
EPAIS=50.0,),
_F(GROUP_MA='Flange',
EPAIS=43.0,),
_F(GROUP_MA='Cones',
EPAIS=10.0,),),
DISCRET=_F(CARA='K_TR_D_N',
GROUP_NO='C_node',
VALE=(0.0,0.0,0.0,0.0,0.0,0.0,),),);

# Model defined
#Assign material.

etc

Kees

Thanks for the reply. It was the error and I feel embarrassed that I did not see it. Vele dank.

Kees

This is my error message

LE MODELE Model CONTIENT :
1 ELEMENT(S) DISCRET
13574 ELEMENT(S) COQUE

!

---

!
! <EXCEPTION> <MODELISA6_10> !
! !
! la maille Mshell ne fait pas partie du maillage Shmesh !
!

---

!

Regards

It seems that this shell model is having the better of me.

I get an error message which I do not understand:

!

---

!
! <EXCEPTION> <ELEMENTS_40> !
! !
! -> L'axe de référence pour le calcul du repère local est normal à un !
! au moins un élément de plaque. !
! -> Risque & Conseil : !
! Il faut modifier l'axe de référence (axe X par défaut) en utilisant !
! ANGL_REP ou VECTEUR. !
!

---

!



My model briefly looks like this: A cylindrical section [Mshell], 850 mm long and 1770 mm dia, axis in the Z direction; on this shell a [Nozzle] in the X-direction with a [Flange] on it's end; On each end of [Mshell] there is a reducing cone [Cones]; The one cone is closed by a flat plate [Cover] and the other cone attached to a section of a smaller cylinder [Bshell], 1575 mm dia; [Presface] includes everything but [Flange].

I guess that the offending commands are one of these:

....................
Shmesh=MODI_MAILLAGE(reuse =Shmesh,
MAILLAGE=Shmesh,
ORIE_NORM_COQUE=_F(GROUP_MA=('Presface','Flange',),),);

...................
Carel=AFFE_CARA_ELEM(MODELE=Model,
COQUE=(_F(GROUP_MA='Bshell',
EPAIS=8.0,),
_F(GROUP_MA=('Mshell','Nozzle','Cones',),
EPAIS=10.0,),
_F(GROUP_MA='Pad',
EPAIS=18.0,),
_F(GROUP_MA='Cover',
EPAIS=50.0,),
_F(GROUP_MA='Flange',
EPAIS=43.0,),),
DISCRET=_F(CARA='K_TR_D_N',
GROUP_NO='C_node',
VALE=(0.0,0.0,0.0,0.0,0.0,0.0,),),);

...................
RESU=MECA_STATIQUE(MODELE=Model,
CHAM_MATER=Matrl,
CARA_ELEM=Carel,
EXCIT=_F(CHARGE=BC,),
OPTION='SIEF_ELGA',);

The error message follows on RESU.

I have tried a few versions of VECTEUR and ANGL_REP but to no avail (admitted, I do not really know what I am doing). In the MESH module the element normals look OK to me.

I will appreciate if somebody can point me in the right direction.

Thank you in advance.