Lagerkräfte bei vorgegebener Verschiebungsbelastung verstehen

Hallo,

habe ein STAT_NON_LINE Modell wo eine Verschiebungsbelastung vorgegeben wird (Dz=50mm).

[code:1]
..
..
Last=AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=_F(GROUP_NO='KnLast',
DZ=50.0,),);

BLOCAGE=AFFE_CHAR_MECA(MODELE=MODELE,
FACE_IMPO=_F(GROUP_MA='FlLagFac',
DX=0,
DY=0,
DZ=0,),);
..
..
..
Solution=STAT_NON_LINE(MODELE=MODELE,
CHAM_MATER=CHMAT,
EXCIT=(_F(CHARGE=BLOCAGE,),
_F(CHARGE=Last,
FONC_MULT=ramp,),),
COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
DEFORMATION='SIMO_MIEHE',
TOUT='OUI',),
INCREMENT=_F(LIST_INST=time,),
CONVERGENCE=_F(ITER_GLOB_MAXI=20,),
ARCHIVAGE=_F(PAS_ARCH=1,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),);

LinRes=MECA_STATIQUE(MODELE=MODELE,
CHAM_MATER=CHMAT,
EXCIT=(_F(CHARGE=BLOCAGE,),
_F(CHARGE=Last,),),);

Solution=CALC_ELEM(reuse =Solution,
MODELE=MODELE,
CHAM_MATER=CHMAT,
RESULTAT=Solution,
PRECISION=1e-06,
OPTION='SIEF_ELNO_ELGA',);

Solution=CALC_NO(reuse =Solution,
RESULTAT=Solution,
OPTION=('SIEF_NOEU_ELGA','FORC_NODA','REAC_NODA',),);

IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=Solution,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','DEPL','SIEF_ELNO_ELGA','SIEF_NOEU_ELGA',),),);

TB_nodf=POST_RELEVE_T(ACTION=(_F(OPERATION='EXTRACTION',
INTITULE='ReactionForces',
RESULTAT=Solution,
NOM_CHAM='REAC_NODA',
TOUT_ORDRE='OUI',
GROUP_NO='FlLagNod',
RESULTANTE=('DX','DY','DZ',),
MOYE_NOEUD='NON',),
_F(OPERATION='EXTRACTION',
INTITULE='NodalForces',
RESULTAT=Solution,
NOM_CHAM='FORC_NODA',
PRECISION=0.000001,
GROUP_NO='FlLagNod',
RESULTANTE=('DX','DY','DZ',),),),
TITRE='[Reaction] Nodal Forces',);

IMPR_TABLE(TABLE=TB_nodf,
FORMAT='TABLEAU',
UNITE=26,
SEPARATEUR=' * ',
TITRE='displacements at nodes on group Nforce',);

FIN(FORMAT_HDF='OUI',
RETASSAGE='OUI',
INFO_RESU='OUI',);
[/code:1]

Weiterhin werden die resultierenden Lagerkräfte ermittelt.



Ich hätte erwartet, dass nur Kräfte in Z-Richtung von bedeutender Größe wären, und das "Lager in x-y Ebene" wie ein Luftkissen arbeitet. Wenn ich bei:
[code:1]
AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=_F(GROUP_NO='KnLast',
DZ=50.0,),);
[/code:1]
keine DX und keine DY angebe, dann kann in dieser Ebene die Gruppe KnLast sich doch frei bewegen?

Habe ich bei den Reaktionskräften was falsch definiert? Warum bekomme ich so große F_x und F_y Kräfte?

Über eine Info würde ich mich freuen.

Gruß Markus<br /><br />Post edited by: MGolbs, at: 2011/03/07 15:04

Hallo Kees,

danke für die Rückmeldung. Eventuell kann eine Skizze helfen.



Wenn ich kein DX und DY angebe, können dann die in DZ=50 zu verschiebenden Knoten (KnLast) in der xy-Ebene sich frei bewegen?

Gruß und Dank Markus<br /><br />Post edited by: MGolbs, at: 2011/03/09 13:24

Hallo Kees,

danke für die Info!

...Und die sind nicht direkt durch KnLast diktiert....

Ich bin immer davon ausgegangen, dass durch &amp;#8721;F_x=0, &amp;#8721;F_y=0, &amp;#8721;F_z=0 die Kräfte an FlLagNod = (FlLagFac) = Kräfte an KnLast.

Habe ich da in Code Aster noch was falsch verstanden?


Gruß und vielen Dank Markus