Shell element problem
- Johan Stadler
- Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 22
- Thank you received: 1
10 years 4 months ago - 10 years 4 months ago #7785
by Johan Stadler
Shell element problem was created by Johan Stadler
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
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
Last edit: 10 years 4 months ago by Johan Stadler. Reason: File still not attached
- kwou
-
- Offline
- Moderator
-
10 years 4 months ago - 10 years 4 months ago #7788
by kwou
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Replied by kwou on topic Re: Shell element problem
Hoi Johan
hard to tell without the mess file, but it might be this:
the second group is maille, is that correct?
Can you give the relevant part of the mess file here?
It is still a pity that we cannot upload files here.
hard to tell without the mess file, but it might be this:
....
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,),
....
the second group is maille, is that correct?
Can you give the relevant part of the mess file here?
It is still a pity that we cannot upload files here.
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Last edit: 10 years 4 months ago by kwou.
The following user(s) said Thank You: Johan Stadler
- Johan Stadler
- Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 22
- Thank you received: 1
10 years 4 months ago #7789
by Johan Stadler
Replied by Johan Stadler on topic Re: Shell element problem
Kees
Thanks for the reply. It was the error and I feel embarrassed that I did not see it. Vele dank.
Thanks for the reply. It was the error and I feel embarrassed that I did not see it. Vele dank.
- kwou
-
- Offline
- Moderator
-
10 years 4 months ago #7790
by kwou
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Replied by kwou on topic Re: Shell element problem
hoi Johan
after a while all the code looks the same.
So donot feel embarrassed, it happens to me all the time.
I still would be interested in the relevant part of the mess file.
Alles goed in Zuid-Afrika?
fijn weekend.
groeten
kees
after a while all the code looks the same.
So donot feel embarrassed, it happens to me all the time.
I still would be interested in the relevant part of the mess file.
Alles goed in Zuid-Afrika?
fijn weekend.
groeten
kees
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
- Johan Stadler
- Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 22
- Thank you received: 1
10 years 4 months ago #7791
by Johan Stadler
Replied by Johan Stadler on topic Re: Shell element problem
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
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
- Johan Stadler
- Topic Author
- Offline
- Junior Member
-
Less
More
- Posts: 22
- Thank you received: 1
10 years 4 months ago #7792
by Johan Stadler
Replied by Johan Stadler on topic Re: Shell element problem
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.
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.
Moderators: catux
Time to create page: 0.163 seconds