Error in my simulation between a 1D Beam and a 2D shell

alitoah · 2023-05-26 12:02:00

hey I have a 1D circular beam and a 2D shell and they have the same diameter I have this error when I do liaison_Mail
I tried to do liaison_elem between poutre and coque but they ask for axe poutre direction but I have a beam circular .
i have this error :
EXCEPTION> <PROJECTION4_55> ║
║ ║
║ Il n'y a pas de mailles à projeter ou en correspondance. ║
║ Dans le cas de l'utilisation de AFFE_CHAR_MECA / LIAISON_MAIL, les mailles maîtres ║
║ doivent avoir la même dimension que l'espace de modélisation : ║
║ - mailles volumiques pour un modèle 3D ║
║ - mailles surfaciques pour un modèle 2D
my com file :
EBUT(LANG='FR')

mesh0_grapht = LIRE_MAILLAGE(identifier='0:1',
UNITE=20)
mesh_st = LIRE_MAILLAGE(identifier='1:1',
UNITE=2)
mesh2 = ASSE_MAILLAGE(identifier='2:1',
MAILLAGE_1=mesh0_grapht,
MAILLAGE_2=mesh_st,
OPERATION='SUPERPOSE')
model0 = AFFE_MODELE(identifier='3:1',
AFFE=(_F(GROUP_MA=('Grapht_face', ),
MODELISATION='DKT',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('stent_face', ),
MODELISATION='POU_D_E',
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh2)
elemprop = AFFE_CARA_ELEM(identifier='4:1',
COQUE=_F(EPAIS=0.02,
GROUP_MA=('Grapht_face', )),
MODELE=model0,
POUTRE=_F(CARA=('R', ),
GROUP_MA=('stent_face', ),
SECTION='CERCLE',
VALE=(0.245, )))
mater_st = DEFI_MATERIAU(identifier='5:1',
ELAS=_F(E=40000.0,
NU=0.3))
fieldmat_st = AFFE_MATERIAU(identifier='6:1',
AFFE=_F(GROUP_MA=('stent_face', 'Grapht_face'),
MATER=(mater_st, )),
MODELE=model0)
mater_st = DEFI_MATERIAU(identifier='7:1',
ELAS=_F(E=1125.0,
NU=0.2))
load = AFFE_CHAR_MECA(identifier='8:1',
Any suggestion please ?

Attachments:
Capture d’écran 2023-05-26 à 12.59.03.png, Size: 109.26 KiB, Downloads: 10

jeanpierreaubry · 2023-05-26 12:57:11

hello

this may be one the case explained in my book (in the signature below)
LIAISON_ELEM is suited to that

jean pierre aubry

alitoah · 2023-05-27 19:20:39

hello jeanpierreaubry , ok thanks .
Can I use a liaison_Mail for that please ? or juste the liaison mail is between 2 shells ?
because I have used a liaison mail and the code run

alitoah · 2023-05-27 20:15:49

when I make a liaison_elem between my shell and my beam , I have this error ; <EXCEPTION> <MODELISA6_43> ║
║ ║
║ il ne faut donner qu'un seul noeud dans le GROUP_NO : stent_face .

i don't understand how to define just one noeud I dont find the way to do that ( because for my group of noeud I define a group of noeud from a line geometry ) . any suggestions ?

alitoah · 2023-05-27 21:29:00

i create one node in my st_face , but now I have this error :
EXCEPTION> <DVP_1> ║
║ ║
║ Erreur de programmation. ║
║ ║
║ Condition non respectée: ║
║ longt.le.long

this is my comm.file:
DEBUT(LANG='FR')

mesh0_gr = LIRE_MAILLAGE(identifier='0:1',
UNITE=20)

mesh_st = LIRE_MAILLAGE(identifier='1:1',
UNITE=2)

mesh2 = ASSE_MAILLAGE(identifier='2:1',
MAILLAGE_1=mesh0_gr,
MAILLAGE_2=mesh_st,
OPERATION='SUPERPOSE')

model0 = AFFE_MODELE(identifier='3:1',
AFFE=(_F(GROUP_MA=('Gr_face', ),
MODELISATION='Q4GG',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('st_face', ),
MODELISATION='POU_D_E',
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh2)

elemprop = AFFE_CARA_ELEM(identifier='4:1',
COQUE=_F(EPAIS=0.02,
GROUP_MA=('Gr_face', )),
MODELE=model0,
POUTRE=_F(CARA=('R', ),
GROUP_MA=('st_face', ),
SECTION='CERCLE',
VALE=(0.245, )))

mater_st = DEFI_MATERIAU(identifier='5:1',
ELAS=_F(E=40000.0,
NU=0.3))

fieldmat_st = AFFE_MATERIAU(identifier='6:1',
AFFE=_F(GROUP_MA=('st_face', 'Gr_face'),
MATER=(mater_st, )),
MODELE=model0)

mater_st = DEFI_MATERIAU(identifier='7:1',
ELAS=_F(E=1125.0,
NU=0.2))

load = AFFE_CHAR_MECA(identifier='8:1',
DDL_IMPO=_F(DRX=0.0,
DRY=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('inf_edge', )),
MODELE=model0)

load0 = AFFE_CHAR_MECA(identifier='9:1',
FORCE_NODALE=_F(FZ=5.0,
GROUP_NO=('sup_edge', )),
MODELE=model0)

load1 = AFFE_CHAR_MECA(identifier='10:1',
LIAISON_ELEM=_F(AXE_POUTRE=(1.0, ),
CARA_ELEM=elemprop,
GROUP_MA_1=('Gr_face', ),
GROUP_NO_2=('st_noeud', ),
OPTION='COQ_POU'),
MODELE=model0)

reslin = MECA_STATIQUE(identifier='11:1',
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat_st,
EXCIT=(_F(CHARGE=load),
_F(CHARGE=load0),
_F(CHARGE=load1)),
MODELE=model0)

IMPR_RESU(identifier='12:1',
FORMAT='GMSH',
MODELE=model0,
RESU=_F(CARA_ELEM=elemprop,
NOM_CHAM=('DEPL', ),
RESULTAT=reslin),
UNITE=3,
VERSION=1.0)

FIN()

#1 2023-05-26 12:02:00

Error in my simulation between a 1D Beam and a 2D shell

#2 2023-05-26 12:57:11

Re: Error in my simulation between a 1D Beam and a 2D shell

#3 2023-05-27 19:20:39

Re: Error in my simulation between a 1D Beam and a 2D shell

#4 2023-05-27 20:15:49

Re: Error in my simulation between a 1D Beam and a 2D shell

#5 2023-05-27 21:29:00

Re: Error in my simulation between a 1D Beam and a 2D shell

Board footer