Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2017-09-14 11:50:04

biscarri
Member
Registered: 2012-01-25
Posts: 44

[SOLVED] joint 3D

I'm trying to model a 3D joint to use it as isostatic mechanical 3D model support.

The support consists of:
- a group containing the nodes laying on the circumference of a bolt hole (group 'support_1')
- a group containing a single node located at the circumference center (group 'N1')
- a LIAISON_RB3 making a rigid link between the node in group 'N1' (master node) and the nodes in group 'support_1' (slave nodes)
- a DDL_IMPO specifying DX=DY=DZ=0 at node in group 'N1'
- a group containing a single 0D element (group 'E1')

When running the case I get the following error message:

ERREUR A L'INTERPRETATION DANS ACCAS - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR phase d'initialisation
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.          !
   ! Les bases sont fermees.                                        !
   ! Type de l'exception : error                                    !
   !                                                                !
   ! Erreur utilisateur :                                           !
   !   Vous voulez contraindre le ddl DY sur un ensemble de noeuds, !
   !   Mais ce ddl n'existe sur aucun de ces noeuds.                !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR phase d'initialisation
>> JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_5545-biscarri-ews=1
<I>_EXIT_CODE = 1

I have tried to add to the model a 0D element in order to provide the node in group 'N1' with the DDL needed, using MODELISATION='3D_JOINT', but I don't find the correct way to make the model work.

Please can you give some hints to model this 3D joint?

It follows the .comm file of the case, thank you very much:

DEBUT(PAR_LOT='NON');
#
# mallat
#
mesh_pl=LIRE_MAILLAGE(FORMAT='MED',UNITE=20,VERI_MAIL=_F(VERIF='NON'));

#mesh_pl=DEFI_GROUP(reuse=mesh_pl,MAILLAGE=mesh_pl,
                  #CREA_GROUP_MA=(
                                 #_F(NOM=('placa'),GROUP_MA=('placa'),),
                                 #),
                  #CREA_GROUP_NO=(
                                 #_F(NOM=('suport_1'),GROUP_NO=('suport_1'),),
                                 #_F(NOM=('suport_2'),GROUP_NO=('suport_2'),),
                                 #_F(NOM=('suport_3'),GROUP_NO=('suport_3'),),
                                 #_F(NOM=('suport_4'),GROUP_NO=('suport_4'),),
                                 #_F(NOM=('N1'),GROUP_NO=('N1'),),
                                 #_F(NOM=('N2'),GROUP_NO=('N2'),),
                                 #_F(NOM=('N3'),GROUP_NO=('N3'),),
                                 #_F(NOM=('N4'),GROUP_NO=('N4'),),
                                 #_F(NOM=('carrega'),GROUP_NO=('carrega'),),
                                #),
                  #);
#
# creacio model mecanic
#
modelM=AFFE_MODELE(MAILLAGE=mesh_pl,Lluis
                 #INFO=2,
                 AFFE=(
                 _F(GROUP_MA=('placa',),
                          PHENOMENE='MECANIQUE',
                          MODELISATION='3D',),                           
                 _F(GROUP_MA=('E1','E2','E3','E4',),
                          PHENOMENE='MECANIQUE',
                          MODELISATION='3D_JOINT',),                           
                          ),
                 );
#
# materials
#
Esteel=DEFI_CONSTANTE(VALE=210.0e9,);
NUsteel=DEFI_CONSTANTE(VALE=0.3,);
steel_M=DEFI_MATERIAU(ELAS_FO=_F(E=Esteel,
                                 NU=NUsteel,
                                 RHO=7.85e3,
                                 ),                                                                                         
                     );

materiaM=AFFE_MATERIAU(MAILLAGE=mesh_pl,
                      MODELE=modelM,
                      AFFE=(
                      _F(GROUP_MA='placa',MATER=steel_M,),
                            ), 
                      );
#
# carregues i condicions de contorn
#
loadM=AFFE_CHAR_MECA(MODELE=modelM,
                        DDL_IMPO=( 
                                 #_F(GROUP_NO=("symXZ"),DY=0.0),                        # simetria pla XZ
                                 #_F(GROUP_NO=("symYZ"),DX=0.0),                        # simetria pla YZ
                                 _F(GROUP_NO=('N1'),DX=0.0,DY=0.0,DZ=0.0),              # suport 1
                                 _F(GROUP_NO=('N2'),DY=0.0,DZ=0.0),                     # suport 2
                                 _F(GROUP_NO=('N3'),DZ=0.0),                            # suport 3
                                 _F(GROUP_NO=('N4'),DZ=0.0),                            # suport 4
                                 #_F(GROUP_MA=('E1'),DX=0.0,DY=0.0,DZ=0.0),              # suport 1
                                 #_F(GROUP_MA=('E2'),DY=0.0,DZ=0.0),                     # suport 2
                                 #_F(GROUP_MA=('E3'),DZ=0.0),                            # suport 3
                                 #_F(GROUP_MA=('E4'),DZ=0.0),                            # suport 4
                                 ),
                        FORCE_NODALE=(
                                 _F(GROUP_NO=('carrega'),FZ=7333.33),
                                     ),
                        LIAISON_RBE3=(                                     # rigiditza circumferenca suport amb node centre
                                     _F(NOEUD_MAIT=('N1'),GROUP_NO_ESCL=('suport_1',),DDL_MAIT=('DX','DY','DZ','DRX','DRY','DRZ'),DDL_ESCL=('DX','DY','DZ')), 
                                     _F(NOEUD_MAIT=('N2'),GROUP_NO_ESCL=('suport_2',),DDL_MAIT=('DX','DY','DZ','DRX','DRY','DRZ'),DDL_ESCL=('DX','DY','DZ')), 
                                     _F(NOEUD_MAIT=('N3'),GROUP_NO_ESCL=('suport_3',),DDL_MAIT=('DX','DY','DZ','DRX','DRY','DRZ'),DDL_ESCL=('DX','DY','DZ')), 
                                     _F(NOEUD_MAIT=('N4'),GROUP_NO_ESCL=('suport_4',),DDL_MAIT=('DX','DY','DZ','DRX','DRY','DRZ'),DDL_ESCL=('DX','DY','DZ')), 
                                     ),
                                  );
#
# solucio
#
try:
                       
    Solucio=MECA_STATIQUE, reuse=Solucio (MODELE=modelM,
                 CHAM_MATER=materiaM,
#                 CARA_ELEM=
                 EXCIT=(_F(CHARGE=loadM),),
                 );
except:
        pass
#
# postproces
#
Solucio=CALC_CHAMP(reuse =Solucio,
                    RESULTAT=Solucio,
                    GROUP_MA='placa',
                    #CONTRAINTE=('SIGM_ELNO',),
                    CRITERES=('SIEQ_ELNO','SIEQ_NOEU',),
                    FORCE=('FORC_NODA','REAC_NODA',),
                     );
IMPR_RESU(FORMAT='MED',
           UNITE=81,
           RESU=(
                _F(RESULTAT=Solucio,
                   MAILLAGE=mesh_pl,
                   NOM_CHAM=('DEPL','REAC_NODA','FORC_NODA','SIEQ_ELNO','SIEQ_NOEU',),
                   ),
           #CONCEPT=(
                    ##_F(CARA_ELEM = realcsM),   # Check real constant
                   ## _F(CHAM_MATER = materiaM), # Check materials
                    #_F(CHARGE= loadM),         # Check loads NE SEMBLE PAS MARCHER
                   #),
                ),
          );
FIN();

Last edited by biscarri (2017-09-26 11:47:52)

Offline

#2 2017-09-20 10:41:13

biscarri
Member
Registered: 2012-01-25
Posts: 44

Re: [SOLVED] joint 3D

Please, can somebody help me with this question?

Thank you very much in advande,
Lluis

Offline

#3 2017-09-21 08:44:11

RichardS
Member
From: Munich, Germany
Registered: 2010-09-28
Posts: 528
Website

Re: [SOLVED] joint 3D

Hi,
I think you are using the wrong element type for your center nodes, it should be a DISCRETE element of type K_T_D_N (T means it has only translational DOFs, if it should also have rotational DOFs, use K_TR_D_N).
You should also be able to find some examples in the forum.

Best,
Richard


Richard Szoeke-Schuller
Product Management
www.simscale.com
We are hiring! https://simscale-jobs.personio.de/?language=en#all

Offline

#4 2017-09-21 13:28:42

biscarri
Member
Registered: 2012-01-25
Posts: 44

Re: [SOLVED] joint 3D

Thank you very much for your answer, Richard.

Best regards,
Lluis

Offline

Board footer