Code Aster Modell mit 1D Beam, 2D Plate and 3D Solid - Gruppen - Netzdefinition

Hallo,

bis jetzt habe ich immer nur Modelle gerechnet, wo nur 2D oder nur 3D Elemente enthalten waren.

Nun möchte ich gern Mixmodelle mit 1D, 2D und 3D Elementen rechnen. Bei der Gruppendefinition und der Netzzuordnung sowie den Freiheitsgraden an den Knoten benötige ich noch Unterstützung.

Ist folgendes so erst einmal richtig:
# pro 1D, 2D und 3D Elementtype (Querschnitt, Material..) muss ich im Netz erst einmal Gruppen definieren?
# für jede Gruppe von 1D und 2D Elementen muss ich Material und Geometrie in der *.comm definieren?
# das gesamte Netz spreche ich dann z.B. über ...=AFFE_MODELE(..) an?
# was muss dann aber in =AFFE_MODELE(... MODELISATION='???'.. stehen?
# was mache ich bei unterschiedlichen Materialien ..=AFFE_MATERIAU( MATER=???,_???_,?
# Wie realisiere ich dann z.B. unter ...=STAT_NON_LINE(
COMP_INCR=_F(RELATION='????',
DEFORMATION='???' für unterschiedliche Elementtypen?
# Wenn ich an einem Knoten eines Solids einen Balken anbinde werden auch im Code-Aster keine Momente übertragen?

# Gibt es ein Beispiel wo verschiedene Elementtypen (1D, 2D und 3D) zusammen ein Netz bilden? Rechnen möchte ich lineare Statik, STAT_NON_LINE sowie Modalanalyse.


Über Tipps und Infos wäre ich sehr dankbar.


Gruß und Dank Markus<br /><br />Post edited by: MGolbs, at: 2011/06/17 10:00

Hallo Markus,
hier ein Beispiel mit 1D-Balken-, 2D-Schalen und 3D-Volumen-Elementen.
Das müsste eigentlich alle deine Fragen beantworten.
MfG,
Richard



DEBUT();

MeshLin=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',);

MeshLin=MODI_MAILLAGE(reuse =MeshLin,
MAILLAGE=MeshLin,
ORIE_NORM_COQUE=_F(GROUP_MA='BEAM_2D',),);

FEMLin=AFFE_MODELE(MAILLAGE=MeshLin,
INFO=2,
AFFE=(_F(GROUP_MA=('BEAM_2D','COUPL',),
PHENOMENE='MECANIQUE',
MODELISATION='DKT',),
_F(GROUP_MA='BEAM',
PHENOMENE='MECANIQUE',
MODELISATION='POU_D_E',),
_F(GROUP_MA='BEAM_3D',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),),);

Car_Elem=AFFE_CARA_ELEM(MODELE=FEMLin,
POUTRE=_F(GROUP_MA='BEAM',
SECTION='CERCLE',
CARA=('R','EP',),
VALE=(10.0,1.0,),),
COQUE=_F(GROUP_MA=('BEAM_2D','COUPL',),
EPAIS=1.0,),);

Steel=DEFI_MATERIAU(ELAS=_F(E=201000.0,
NU=0.3,),);

Material=AFFE_MATERIAU(MAILLAGE=MeshLin,
AFFE=_F(TOUT='OUI',
MATER=Steel,),);

LinRamp=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,0.0,
1.0,1.0,
),INTERPOL='LIN',PROL_DROITE='LINEAIRE',PROL_GAUCHE='EXCLU',);

Bcn=AFFE_CHAR_MECA(MODELE=FEMLin,
DDL_IMPO=_F(GROUP_NO=('BCN_3D','BCN_SH',),
LIAISON='ENCASTRE',),
LIAISON_SOLIDE=_F(GROUP_NO='LIAS_SOL',),
LIAISON_ELEM=_F(OPTION='COQ_POU',
GROUP_MA_1='COUPL',
NOEUD_2='N10394',
CARA_ELEM=Car_Elem,
AXE_POUTRE=(0.0,0.0,-1.0,),),);

Force=AFFE_CHAR_MECA(MODELE=FEMLin,
FORCE_NODALE=_F(NOEUD='N10395',
FY=-500.0,),);

Listres=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1.0,
PAS=0.05,),);

Solution=STAT_NON_LINE(MODELE=FEMLin,
CHAM_MATER=Material,
CARA_ELEM=Car_Elem,
EXCIT=(_F(CHARGE=Bcn,),
_F(CHARGE=Force,
FONC_MULT=LinRamp,),),
COMP_INCR=(_F(RELATION='ELAS',
DEFORMATION='PETIT',
GROUP_MA=('BEAM_3D','BEAM_2D',),),
_F(RELATION='ELAS',
DEFORMATION='PETIT',
GROUP_MA='BEAM',),),
INCREMENT=_F(LIST_INST=Listres,
SUBD_METHODE='UNIFORME',
SUBD_PAS=2,
SUBD_PAS_MINI=0.005,),
NEWTON=_F(MATRICE='TANGENTE',
REAC_ITER=1,),
CONVERGENCE=_F(ITER_GLOB_MAXI=100,),
SOLVEUR=_F(METHODE='MULT_FRONT',),
ARCHIVAGE=_F(PAS_ARCH=2,),);

Solution=CALC_ELEM(reuse =Solution,
RESULTAT=Solution,
OPTION='SIEF_ELNO_ELGA',);

IMPR_RESU(MODELE=FEMLin,
RESU=_F(RESULTAT=Solution,
NOM_CHAM='DEPL',),);

FIN();

Hallo Richard,

danke für den Code. Ich werde es mal an einem einfachen Beispiel versuchen umzusetzen und anzupassen.

Gruß und Dank Markus