Hello,
I am currently trying to simulate a modal analysis with two glued bodies and two different materials.
I have managed to glue two meshes in a modal analysis, but I am not able to assign different materials to the two meshes. With ASSEMBLAGE I can only select one material field, if I double ASSEMBLAGE I get an error message. Does anyone have an idea or an example code for me?
Best regards and many thanks!
DEBUT(LANG='EN')
mesh_Traeger = LIRE_MAILLAGE(FORMAT='MED',
UNITE=3)
mesh_Spiegel = LIRE_MAILLAGE(FORMAT='MED',
UNITE=4)
Gesamtnetz = ASSE_MAILLAGE(MAILLAGE_1=mesh_Traeger,
MAILLAGE_2=mesh_Spiegel,
OPERATION='SUPERPOSE')
model = AFFE_MODELE(AFFE=_F(MODELISATION='3D',
PHENOMENE='MECANIQUE',
TOUT='OUI'),
MAILLAGE=Gesamtnetz)
Glas = DEFI_MATERIAU(ELAS=_F(E=80000.0,
NU=0.23,
RHO=2.5e-09))
zuw_Glas = AFFE_MATERIAU(AFFE=_F(MATER=(Glas, ),
TOUT='OUI'),
MAILLAGE=mesh_Spiegel)
Aluminium = DEFI_MATERIAU(ELAS=_F(E=200000.0,
NU=0.3,
RHO=2.5e-09))
zuw_Aluminium = AFFE_MATERIAU(AFFE=_F(MATER=(Glas, ),
TOUT='OUI'),
MAILLAGE=mesh_Traeger)
Fix = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('FIX', )),
MODELE=model)
Kontakt = AFFE_CHAR_MECA(LIAISON_MAIL=_F(DISTANCE_MAX=0.1,
GROUP_MA_ESCL=('11', '22'),
GROUP_MA_MAIT=('Rott', ),
TYPE_RACCORD='MASSIF'),
MODELE=model)
ASSEMBLAGE(CHAM_MATER=zuw_Aluminium,
CHARGE=(Fix, Kontakt),
MATR_ASSE=(_F(GROUP_MA=('Rott', ),
MATRICE=CO('rigidite'),
OPTION='RIGI_MECA'),
_F(GROUP_MA=('Rott', ),
MATRICE=CO('masse'),
OPTION='MASS_MECA')),
MODELE=model,
NUME_DDL=CO('numeddl'))
ASSEMBLAGE(CHAM_MATER=zuw_Glas,
CHARGE=(Kontakt, ),
MATR_ASSE=(_F(GROUP_MA=('S1', 'S2'),
MATRICE=CO('rigidite'),
OPTION='RIGI_MECA'),
_F(GROUP_MA=('S1', 'S2'),
MATRICE=CO('masse'),
OPTION='MASS_MECA')),
MODELE=model,
NUME_DDL=CO('numeddl'))
result2 = CALC_MODES(CALC_FREQ=_F(FREQ=(10.0, 10000.0)),
MATR_MASS=masse,
MATR_RIGI=rigidite,
OPTION='BANDE')
result1 = CALC_MODES(CALC_FREQ=_F(FREQ=(10.0, 10000.0)),
MATR_MASS=masse,
MATR_RIGI=rigidite,
OPTION='BANDE')
IMPR_RESU(FORMAT='MED',
RESU=(_F(RESULTAT=result2),
_F(RESULTAT=result1)),
UNITE=80)