
dimension 3
niveau_commentaire 3

TYPE_DE_CALCUL
dynamique_relaxation_dynam avec plus sauveMaillagesEnCours
PARA_TYPE_DE_CALCUL

typeCalRelaxation= 4 lambda= 1 type_calcul_mass= 2 option_recalcul_mass= 0
parametre_calcul_de_la_masse_ casMass_relax= 3
lambdaAvecPonderation_Globale_ f_lambda
propCinetiqueAvecPonderation_Globale_ f_propCinetique
parametre_calcul_de_la_viscosite_  type_calcul_visqu_critique= 1 opt_cal_C_critique= 1  f_= 0.9

avec_amortissement_cinetique_
   max_nb_decroit_pourRelaxDyn_   1
   taille_moyenne_glissante_      1
   coef_arret_pourRelaxDyn_       0.
   coef_redemarrage_pourRelaxDyn_ 0.
   max_deltaX_pourRelaxDyn_       0.
   nb_deb_test_amort_cinetique_   1
   nb_max_dX_OK_pourRelaxDyn_     1 # 10
   nb_deb_testfin_pourRelaxDyn_   1 # 250
fi_parametre_amortissement_cinetique_

ARRET_A_EQUILIBRE_STATIQUE_ 2
mode_debug_= 100
#------------------------------------------------------------------ 

< maillage.her
renumerotation_des_noeuds_

# def de mouvements solides
def_mouvement_solide_initiaux_
mouvement_solide_ 
 translation_= 0.1000 0.0000 0.0000
 centre_=      0.1000 0.0000 0.0000
 rotation_=     0 -1.56434480  0 en_radian_
fin_mouvement_solide_

#------------------------------------------------------------------ 
    les_courbes_1D ------------
# Fonction de pondération pour le loi_melange
f_ponderation COURBEPOLYLINEAIRE_1_D
  Debut_des_coordonnees_des_points
    Coordonnee dim= 2  0.  1.
    Coordonnee dim= 2  1.  1.
    Coordonnee dim= 2  2.  1.
    Coordonnee dim= 2  3.  1.
  Fin_des_coordonnees_des_points

# Altitude z0 correspondant à une pression nulle
f_z0 COURBEPOLYLINEAIRE_1_D
  Debut_des_coordonnees_des_points
    Coordonnee dim= 2  0.   0.
    Coordonnee dim= 2  1.   0.
    Coordonnee dim= 2  2.  10.00
    Coordonnee dim= 2  3.  10.00
  Fin_des_coordonnees_des_points 
  
 
    les_fonctions_nD ------------

# Fonction pour le pilotage du lambda
f_lambda FONCTION_EXPRESSION_LITTERALE_nD
  un_argument= amor_cinet_visqueux  
  un_argument= temps_courant  
  fct= (temps_courant <= 1.) ? \ 
        (amor_cinet_visqueux == 1) ?  0.600 :  0.900  \
      : (amor_cinet_visqueux == 1) ?  0.600 :  3.000  
fin_parametres_fonction_expression_litterale_  

# Fonction pour le pilotage du bascule cinétique -> visqueux
f_propCinetique FONCTION_EXPRESSION_LITTERALE_nD
  un_argument= norme_de_convergence
  un_argument= temps_courant
fct= (temps_courant <= 1) ? \
        (norme_de_convergence <  0.0e+00) ? 0. : 1. \
      : 1  
fin_parametres_fonction_expression_litterale_     
  
# calcul du volume total
le_volume FONCTION_EXPRESSION_LITTERALE_nD 
   deb_list_var_   X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_   
   fct=     -(X1*N_surf_1 + X2*N_surf_2 + X3*N_surf_3)/3
   fin_parametres_fonction_expression_litterale_   

f_radiale FONCTION_EXPRESSION_LITTERALE_nD 
   deb_list_var_   X1 X2 X3 fin_list_var_   
   fct=     X1/sqrt(X1^2+X2^2), X2/sqrt(X1^2+X2^2), 0
   fin_parametres_fonction_expression_litterale_      

    
#------------------------------------------------------------------ 

    choix_materiaux ---------
    
E_quad	     loi_film
E_biel      loi_critere_sangle
E_biel_pole loi_sangleBC

    materiaux ----------
        
#-- Loi sangle bielle
loi_sangleBC ISOELAS1D
 3.550e+08  0.3
 
#-- Loi elastique 1D avec critere pli
loi_critere_sangle LOI_CRITERE
TYPE_DE_CRITERE_ PLISSEMENT_BIEL
    ISOELAS1D
    3.550e+08  0.3
    fin_loi_interne
fin_loi_critere 
 
#-- Loi film
loi_film  LOI_DES_MELANGES_EN_SIGMA
  les_grandeurs_de_controle_= TEMPS  fin_grandeurs_
  deb_fonct_= f_ponderation fin_fonct_
  
    #-- première loi
    ISOELAS2D_C demarre_a_prop_non_nulle_puis_strictement_decroissante_
        2.580e+08  0.30
    
   #-- seconde loi
    LOI_CRITERE demarre_a_prop_non_nulle_puis_strictement_decroissante_
    TYPE_DE_CRITERE_  PLISSEMENT_MEMBRANE
        LOI_CONTRAINTES_PLANES 
            NEWTON_LOCAL  avec_parametres_de_reglage_
                nb_iteration_maxi_  20
                nb_dichotomie_maxi_ 20
                tolerance_residu_ 1e1  # tolerance residu 10 Pa
                tolerance_residu_rel_ 1.e-4
                maxi_hsurh0_  10. #4
                mini_hsurh0_  0.01 #0.25
                permet_affichage_ 0#5 
            fin_parametres_reglage_Algo_Newton_
            
           ISOELAS
           2.5800e+08  0.30
        fin_loi_contrainte_plane
    fin_loi_interne
    parametres_controle_pour_contraintes_planes_double_
      NEWTON_LOCAL avec_parametres_de_reglage_
                    nb_iteration_maxi_  20 
                    nb_dichotomie_maxi_  20 
                    tolerance_residu_ 1e1
                    tolerance_residu_rel_ 1.e-4 
                    mini_hsurh0_ 1.e-2 
                    mini_bsurb0_ 1.e-2 
                    maxi_hsurh0_ 1. 
                    maxi_bsurb0_ 1. 
                    sortie_post_ 1 
                    permet_affichage_ 0#4 
                 fin_parametres_reglage_Algo_Newton_ 
    fin_parametres_controle_pour_contraintes_planes_double_
	fin_loi_critere 
fin_liste_lois_elementaires
#-------------------------------------------------------- 

    masse_volumique ---------
E_quad       9.130000e+02 # en kg/m^3
E_biel       8.300000e+02 # en kg/m^3
E_biel_pole  8.300000e+02 # en kg/m^3

    epaisseurs ------------
nom_mail= ballon E_quad_film    2.500000e-05


    sections  ---------
nom_mail= ballon E_biel_bord    4.793000e-06
nom_mail= ballon E_biel_int  9.586000e-06
nom_mail= ballon E_biel_pole  1.000000e-03 

    integrale_sur_volume_
E_quad une_fonction_nD_ le_volume

#------------------------------------------------------------------ 

    charges ---------
 # gradient de pression
F_tout  PHYDRO  0. 0. 1. \
                       0. 0.  nom_Xi_ f_z0 \
                       7.00000000 ATTRIBUT_ sans_limitation_ TEMPS_MINI= 1 

N_f_radiales PONCTUELLE 10, 10, 0 Fonction_nD_CHARGE: f_radiale TEMPS_MAXI= 2
N_poleN PONCTUELLE 0, 0, 4000 TEMPS_MAXI= 1
                       
 # poids propre du a la gravite
nom_mail= ballon E_quad_film VOLUMIQUE  0. 0. -8.862618e+03 TEMPS_MINI= 1
nom_mail= ballon E_biel VOLUMIQUE 0. 0. -8.056925e+03 TEMPS_MINI= 1


 # poids du pole
N_poleN_int PONCTUELLE 0. 0. -0.603417 TEMPS_MINI= 1
N_poleN_ext PONCTUELLE 0. 0. -0.301708 TEMPS_MINI= 1
#------------------------------------------------------------------ 

    blocages ---------

N_poleS UX UY UZ
N_poleN UX UY
# Etape 6 : blocage
N_tout UX UY UZ TEMPS_MINI= 3.01

#<bloc_blocages_decactive>



#------------------------------------------------------------------ 
    typecharge -----------
TYPE3

    controle ------------ 
TEMPSFIN    1
DELTAt 	    1
DELTAtMINI  1
DELTAtMAXI  1
ITERATIONS 100000
NORME E_cinetique/E_statique_ET_Res/Reac_et_Fext
PRECISION 1.00e-02
SAUVEGARDE  1
RESTART 0

    para_affichage ------------ 
FREQUENCE_AFFICHAGE_INCREMENT  1
FREQUENCE_AFFICHAGE_ITERATION  1
FREQUENCE_SORTIE_FIL_DU_CALCUL 1

para_pilotage_equi_global
INIT_INCRE_AVEC_PAS_PREC 0

resultats  pas_de_sortie_finale_ 
COPIE 0   

_suite_point_info_
    controle ------------ 
TEMPSFIN    4
PRECISION   1.00e-03   
RESTART 0

  _fin_point_info_