Constantes_utilisateurs_
< constantes.hz_courbe
fin_Constantes_utilisateurs_

dimension 3

niveau_commentaire 3

TYPE_DE_CALCUL
dynamique_explicite # dynamique_relaxation_dynam # avec plus visualisation

PARA_TYPE_DE_CALCUL

#note : lambda=0.554 revient a choisir 0.95 fois le pas de temps critique DFC dans un calcul explicite classique
#typeCalRelaxation= 1   lambda= 0.554   type_calcul_mass= 2 option_recalcul_mass= 1

#parametre_calcul_de_la_masse_ alpha= 1 beta= 0 gamma= 0 delta= 0  casMass_relax= 3

avec_amortissement_cinetique_
   max_nb_decroit_pourRelaxDyn_   1
   taille_moyenne_glissante_      1
   coef_arret_pourRelaxDyn_       0.
   coef_redemarrage_pourRelaxDyn_ 0.
   max_deltaX_pourRelaxDyn_       0.1
   nb_max_dX_OK_pourRelaxDyn_     1
   nb_deb_testfin_pourRelaxDyn_   1
   nb_deb_test_amort_cinetique_   1
fi_parametre_amortissement_cinetique_

ARRET_A_EQUILIBRE_STATIQUE_ 2


< flan_1D_axi.her



les_courbes_1D

courbe_pression COURBEPOLYLINEAIRE_1_D
    Debut_des_coordonnees_des_points
      Coordonnee dim= 2 0. 0.
      Coordonnee dim= 2 1. 1.
    Fin_des_coordonnees_des_points


les_fonctions_nD
< fonction_potentiel.hz_courbe


choix_materiaux
E_tout MAT


materiaux
MAT____ ISOELAS2D_C
200 0.45


MAT__  LOI_CONTRAINTES_PLANES

 NEWTON_LOCAL  avec_parametres_de_reglage_
                 nb_iteration_maxi_  30
                 nb_dichotomie_maxi_  5
                 tolerance_residu_ 1.e-6
                 tolerance_residu_rel_ 1.e-5
                 maxi_hsurh0_  4
                 mini_hsurh0_  0.25
              fin_parametres_reglage_Algo_Newton_

  LOI_ADDITIVE_EN_SIGMA

    ISOHYPERBULK3
    # 3K     
    30000.


    MAHEO_HYPER
      # ----------------------------------------------------------------------------------
      # |... loi de comportement hyper elastique 3D Maheo_hyper                       .. |
      # ----------------------------------------------------------------------------------
      Qsig_rev= 20.   mu1_rev= 200.  mu2_rev= 5. \
      avec_regularisation_  1.e-12
#      Qsig_rev= 10.   mu1_rev= 100.  mu2_rev= 5. \

  fin_liste_lois_elementaires

fin_loi_contrainte_plane
permet_affichage_ -1

#materiau : MAT
< loi_cpt.hz_loi



masse_volumique
E_tout   1.e-9

epaisseurs
E_tout 50.e-3


charges
# 0.0001 = 1 hPa
# 0.001  = 10 hPa
# 0.01   = 100 hPa
A_pression PRESSION 0.0783 COURBE_CHARGE: courbe_pression ECHELLE: 1


blocages

N_tout UZ

#flan
N_sym UX
N_rayon UX UY



controle

#DELTAtMINI 0.01
#DELTAtMAXI 0.01
TEMPSFIN 99999.
DELTAt COEF_PASCRITIQUE 0.6#0.01
ITERATIONS 999999999
PRECISION 1.e-2
SAUVEGARDE   0
MAXINCRE  500
NORME Residu/Reaction
#RESTART    ---


para_pilotage_equi_global


para_syteme_lineaire
#OPTIMISATION_POINTEURS_ASSEMBLAGE 1
#TYPE_MATRICE BANDE_SYMETRIQUE_LAPACK
#MATRICE_S_SECONDAIRE_S_   BANDE_SYMETRIQUE   FIN_TYPE_MAT_SECONDAIRES_
#TYPE_RESOLUTION_S_SECONDAIRE_S_  CHOLESKY  FIN_TYPE_RESOLUTION_SECONDAIRES_

para_affichage
FREQUENCE_SORTIE_FIL_DU_CALCUL  1


resultats pas_de_sortie_finale_
COPIE 0

_fin_point_info_
