#1 2021-07-29 10:49:19

ravi6195

Member

Registered: 2021-07-29

Posts: 12

Glued contact between shell and solid

Hello, I'm new to Salome Meca and I'm trying to make glued contact between one plate(shell element) and two supports(solid element). I'm getting these errors and I do not understand what does it means or how to solve this.
!------------------------------------------------------------------!
   ! <EXCEPTION> <MODELISA_30>                                        !
   !                                                                  !
   !  you cannot assign values of the SHELL" type" to the model model !
   !  which does not contain only one element shell                   !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------!
   ! <EXCEPTION> <MODELISA5_57>                                 !
   !                                                            !
   !  error (S) met (S) during the checking of the assignments. !
   !------------------------------------------------------------!
I'm attaching the folder where the files have been saved for reference and here's the script from text mode.

pmesh = LIRE_MAILLAGE(FORMAT='MED',
                      UNITE=20)

s1mesh = LIRE_MAILLAGE(FORMAT='MED',
                       UNITE=2)

s2mesh = LIRE_MAILLAGE(FORMAT='MED',
                       UNITE=3)

p_s1mesh = ASSE_MAILLAGE(MAILLAGE_1=pmesh,
                         MAILLAGE_2=s1mesh,
                         OPERATION='SUPERPOSE')

MESH = ASSE_MAILLAGE(MAILLAGE_1=p_s1mesh,
                     MAILLAGE_2=s2mesh,
                     OPERATION='SUPERPOSE')

model = AFFE_MODELE(AFFE=(_F(GROUP_MA=('s_2', ),
                             MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('s_1', ),
                             MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('plate', ),
                             MODELISATION=('DKTG', ),
                             PHENOMENE='MECANIQUE')),
                    MAILLAGE=MESH)

elemprop = AFFE_CARA_ELEM(COQUE=_F(EPAIS=3.0,
                                   GROUP_MA=('plate', )),
                          MODELE=model)

mater = DEFI_MATERIAU(ELAS=_F(E=70000.0,
                              NU=0.18))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=MESH)

load = AFFE_CHAR_MECA(DDL_IMPO=(_F(DX=0.0,
                                   DY=0.0,
                                   DZ=0.0,
                                   GROUP_NO=('fix1_2', )),
                                _F(DX=0.0,
                                   DY=0.0,
                                   DZ=0.0,
                                   GROUP_NO=('fix1_1', ))),
                      MODELE=model)

load0 = AFFE_CHAR_MECA(FORCE_FACE=_F(FZ=0.0054,
                                     GROUP_MA=('plate', )),
                       MODELE=model)

CHNOR = CREA_CHAMP(GROUP_MA=('plate', ),
                   MODELE=model,
                   OPERATION='NORMALE',
                   TYPE_CHAM='NOEU_GEOM_R')

cont_2 = AFFE_CHAR_MECA(LIAISON_MAIL=(_F(CHAM_NORMALE=CHNOR,
                                         EPAIS=3.0,
                                         GROUP_MA_ESCL=('contact2_2', ),
                                         GROUP_MA_MAIT=('s_2', ),
                                         TYPE_RACCORD='COQUE_MASSIF'),
                                      _F(CHAM_NORMALE=CHNOR,
                                         EPAIS=3.0,
                                         GROUP_MA_ESCL=('contact2_1', ),
                                         GROUP_MA_MAIT=('s_1', ),
                                         TYPE_RACCORD='COQUE_MASSIF')),
                        MODELE=model)

reslin = MECA_STATIQUE(CARA_ELEM=elemprop,
                       CHAM_MATER=fieldmat,
                       EXCIT=(_F(CHARGE=load),
                              _F(CHARGE=load0),
                              _F(CHARGE=cont_2)),
                       MODELE=model)

IMPR_RESU(RESU=_F(NOM_CHAM=('DEPL', ),
                  RESULTAT=reslin),
          UNITE=80)

jmrsalonga

Member

Registered: 2022-01-20

Posts: 2

Re: Glued contact between shell and solid

I encounter quite a similar error: you cannot assign values of the SHELL" type" to the model model !
   !  which does not contain only one element shell"

but in my case I was assigning dimension to a beam "POUTRE"

I wonder if you already figured out a way to correct your code?

Last edited by jmrsalonga (2022-01-21 02:19:50)

jmrsalonga

Member

Registered: 2022-01-20

Posts: 2

Re: Glued contact between shell and solid

Hello again, so it turns out our mistake  lies within the model assignment part of our elements. In my case, the order of model assignment ruined it. I assigned a certain element as beam "POU_D_E", and then made an assignment of 3D for everything else...I realized that the latter removed my "POU_D_E" assignment, so when I assigned beam properties to the element that is supposed to be a beam through "POUTRE", I got an error.

In your case, I believe you should have assigned a shell model to your shell elements s_1 and s_2 instead of 3D. In Code Aster, that would be "COQUE_3D"